Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,171,417 G→A L850L (CTG→TTG)  mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,171,4170GA98.8% 249.5 / NA 82L850L (CTG→TTG) mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (45/35);  minor base C (1/0);  total (46/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

TCAAGACGTTTTTGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAGTTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAGGTCACGCATTCC  >  NC_000913/1171173‑1171660
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
tcatgACGTTTTTGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCTACCGCGGAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTGGATGATGGAAATCATTCATCACCCGTTCCAATtcgc                                                                                                                                                                                                                                                <  2:386380/246‑1 (MQ=255)
tCAAGACGTTTTTGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCATCCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATtcg                                                                                                                                                                                                                                                 <  2:538312/249‑1 (MQ=255)
tCAAGACGTTTTTGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCACCCGCGTAACTGGAGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATAATTCATCACCCGTTCCAATtcgcg                                                                                                                                                                                                                                               >  1:679510/1‑251 (MQ=255)
  aaGACGTTTTTGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATtcgcgct                                                                                                                                                                                                                                             <  1:619514/251‑1 (MQ=255)
   aGACGTTTGTGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGAAGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATtcgcgctc                                                                                                                                                                                                                                            <  2:264018/251‑1 (MQ=255)
            tgtgCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTg                                                                                                                                                                                                                                   >  1:162426/1‑251 (MQ=255)
            tgtgCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCACGGTTTCGATATTGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTg                                                                                                                                                                                                                                   >  2:434566/1‑251 (MQ=255)
               gCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGccc                                                                                                                                                                                                                                <  1:673583/251‑1 (MQ=255)
                      tAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGCTGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCCCGTTGAACCGTTGATGATGGAAATCATTCATAACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCg                                                                                                                                                                                                                         >  2:664143/1‑251 (MQ=255)
                           aTCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGc                                                                                                                                                                                                                    >  1:666706/1‑251 (MQ=255)
                           aTCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGc                                                                                                                                                                                                                    <  1:149933/251‑1 (MQ=255)
                              gCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGAt                                                                                                                                                                                                                 >  1:91767/1‑251 (MQ=255)
                                      actcggggtCAGCAACCATGCCTAGGCCTGATGATGCTAACGTCTGACGCGTCAGCGGAACTTGTGCAGCTGCTCCAGACCGAAGTGATCCGCGCGGTCAATGATAATAGTGTTGGCTGTCGGGATGGCGATCCCGGTTTCGATAATGTGTGTACAAACCAGCACGTTGAAACGTTGCTGATGGAAAGCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCtt                                                                                                                                                                                                         <  2:441871/244‑1 (MQ=255)
                                        gCGGTGTCAGCAACCATGCATATGCCTGCTGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGCTGTCGATCCCGGTTTCGCTAATGGTAGTACAAACCAGCAAGTTGAAACGTTGATGATGGAAATCATTCTTCCACCGTTCCAATTCACGCTCGCGCATCTGCCCGTGACCGATGACAATCCGCGCTTcg                                                                                                                                                                                                       >  2:550816/1‑250 (MQ=255)
                                        gCGGTGTCAGCAACCATGCAGATGCCTGATGATGCCAGCATCCGACGCTACCGCGTCACTGCTGCAGCTTCTCCAGACCGAAGTGATTCGAGCGTTCAATGATTATAGTTGTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAAGCATTCATCACCCGTTCCAATTCGCGATCGCTCATCTGCCCGTGACCGATGGCTATCCGCGCGTCt                                                                                                                                                                                                       <  2:262841/251‑1 (MQ=255)
                                          ggTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTgg                                                                                                                                                                                                     >  1:653881/1‑251 (MQ=255)
                                           gtgtCAGCAACCATGCATATGCCTGATGATGCGCACGTCCGACGCGACCGCGTCACTGGTGCAGCTGCGCCAGACCGCAGTCATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCAGGGTTCAATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGAAAATACTTCATCCCCCGTTCCAATTCGCGATCGCGCATCTGCCAGTGCCCGTTGGCGACCCGCGCTTCGGGc                                                                                                                                                                                                    >  2:276512/1‑251 (MQ=255)
                                                aacatccATGCATATGCCTGATGATGCGAACGTCCGACGCGACCTCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATTGAAATCATTCAGCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATTGCGATCCGCGCTTCTGGCCCCAg                                                                                                                                                                                               <  2:426710/246‑1 (MQ=255)
                                                   aaCCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGCTAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTc                                                                                                                                                                                            >  2:294264/1‑251 (MQ=255)
                                                     ccATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGCTCCGCGCTTCTGGCACCAGTTCTg                                                                                                                                                                                          >  2:56419/1‑251 (MQ=255)
                                                     ccATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGCGCAGCTGAGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGTTGGAAATCATTCATCACCCGTTCCAATTCTCGCGCGCGCATCTGCCCGTGCCCGAGGGCGATCCGCGCTTCTGGCACCAGTTCTg                                                                                                                                                                                          >  2:573806/1‑251 (MQ=255)
                                                                   gatgatGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTcggcg                                                                                                                                                                            <  1:301516/251‑1 (MQ=255)
                                                                   gatgatGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTcggc                                                                                                                                                                             <  1:468155/250‑1 (MQ=255)
                                                                    gtgatgCCAACGGCCGCGGCGACCGCGTAACTGGTGCAGCTGCGTCAGACCGAGGTGTGCCGCGTGTGCAATGATAATAGTGGTGGCCGTCGGGATGGCGATCCCGGTGTCGATAATGGTTGTACAAACCCGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAAGTCGCGCTCGCGCATCTGCCCTTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTcggcgt                                                                                                                                                                           <  2:229444/250‑2 (MQ=255)
                                                                    atggtgagaagcaccGACGCTACCGCGTAACTTGTGCAGCTGCGTCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTGGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGGTGATGATGGAAAGCATTCATCACCCGTGCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTGCTGGCACCAGTTCTGCCAGCCGTTcggcgg                                                                                                                                                                           <  2:81930/238‑1 (MQ=255)
                                                                           gAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACAGAAGTGATCCGCGCGTTCAATGATAATCGTGTTGGCTGTCGGGATGTCGCTCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCTATCCGCGCTTCTGTCACCAGTCCTGCCAGCCGTTCGGCGGCTTTCTg                                                                                                                                                                    >  2:497714/1‑251 (MQ=255)
                                                                               gTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCAACACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATg                                                                                                                                                                >  1:490334/1‑251 (MQ=255)
                                                                                  cgacgcgaccgtttatacggTGCAGGTGCGCCAGTCCGAAGTGATCCGCGCGTTCAAGGATAATAGTGTTTGCTGTCGGGATGTCGATCCCGGGTTCGATAATGGTTGTACAAACCAGCCCGTTGAAACGTTGAGGATGGAAATCATTCATCACCCGTTCCATTTCGCGCTCGCGCATCTGCCCGTGACCTATGGCGTTCCTCGTTTCTGTCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGttt                                                                                                                                                             <  2:362029/233‑1 (MQ=255)
                                                                                                   tGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGCTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGAATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCCCCAGTTCTGcc                                                                                                                                                                                        >  1:79000/1‑207 (MQ=255)
                                                                                                   tGGTGCAGCTGCGCCAGACCGAAGTGATCCGAGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGCAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCACTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTagcg                                                                                                                                            >  2:292015/1‑249 (MQ=255)
                                                                                                    ggTGCAGCTGCTCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCCGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAATCGTTGATGATGGAAATCATTCATCACCCGTTCCTATGCGTGCTCGCGCATCTGCTCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCTGTTCGGCGGCTTGCTGAATGTTTTCCACATCATTGTAgaga                                                                                                                                           <  1:49748/251‑1 (MQ=255)
                                                                                                     gTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGat                                                                                                                                          <  1:654416/251‑1 (MQ=255)
                                                                                                         aGCTGCGCCAGACCGAAGTGATCCGCGCGTGCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGataata                                                                                                                                      <  1:154479/251‑1 (MQ=255)
                                                                                                          gCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGataata                                                                                                                                      >  1:211733/1‑250 (MQ=255)
                                                                                                           cTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGCTGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGAtaataa                                                                                                                                     >  2:332122/1‑250 (MQ=255)
                                                                                                             gcgcCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCCCATCATTGTAGAGATAATAAAcc                                                                                                                                  >  2:219970/1‑251 (MQ=255)
                                                                                                                   gACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCTCTCGCGCATCTg                                                                                                                                                                                                                                   <  1:491886/148‑1 (MQ=255)
                                                                                                                   gACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTg                                                                                                                                                                                                                                   >  2:491886/1‑148 (MQ=255)
                                                                                                                    aCCGAAGTGATCCGCGTGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGtcctc                                                                                                                           <  2:344534/251‑1 (MQ=255)
                                                                                                                    aCCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGtcctc                                                                                                                           <  1:434566/251‑1 (MQ=255)
                                                                                                                           tGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGACACGTTGATGATGGAAATCATTAATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGCTGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCCTTTGAGAGCTAATAAACCTGTCCTCAGCGCa                                                                                                                     >  2:560361/1‑250 (MQ=255)
                                                                                                                             aTCCGAGCGTTCCATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGGAACGTTGATGATGGCAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGACTGTTTTACACATCATTGTAAAGATACTAAACCTGGCCTCCGCGCaca                                                                                                                   >  2:415323/1‑248 (MQ=255)
                                                                                                                                   gcgTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTc                                                                                                                                                                                            >  2:627209/1‑171 (MQ=255)
                                                                                                                                   gcgTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTc                                                                                                                                                                                            <  1:627209/171‑1 (MQ=255)
                                                                                                                                       tCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGATTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCTATGGCGATCCGCGCTTCTGGCACCCGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAGTAAACCTGTCCACCGCGCAAACTTTCACGc                                                                                                          >  2:655394/1‑249 (MQ=255)
                                                                                                                                         aaTGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGATTCGATAATGGTTGTACACACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCCCCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCCGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAACCCTGTCCTCCGCGCAAAATTTCCCGCAggg                                                                                                      >  2:461367/1‑250 (MQ=255)
                                                                                                                                           tgataacagagatGGCTGTCGGGCTGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAACCGTTGATGATGGCAATCATTCATCACCCGTTCCAATTCGCGCCCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGTATGTTTTCCACCTCATTGTAGAGATCATACACCTGTCCTCCGCGCAAAATTTCCCGCAGGacc                                                                                                    >  2:279911/13‑249 (MQ=255)
                                                                                                                                                  aGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTccc                                                                                             <  1:56419/251‑1 (MQ=255)
                                                                                                                                                   gtgtTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCg                                                                                            <  1:570947/251‑1 (MQ=255)
                                                                                                                                                    tgtTGGCTGTCGGGATGTCGATCCCGGTTTAGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCACGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAACACTGTCCTCCGCGCAAAATTTCCAGCAGGATCGCATCCCgg                                                                                           >  2:462735/1‑251 (MQ=255)
                                                                                                                                                        cgttcgcGGGATGCCGATCCCGGTTTCAATAATGGTTGTACAAACCATCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTTCGCGCAAAATTTCACGCAGGATCGCCTCCCGGacca                                                                                       <  2:465602/245‑1 (MQ=255)
                                                                                                                                                           tGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCa                                                                                                             >  1:573423/1‑226 (MQ=255)
                                                                                                                                                           tGTCGGGATGGCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCa                                                                                                             <  2:573423/226‑1 (MQ=255)
                                                                                                                                                              cGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGacca                                                                                       >  1:465602/1‑245 (MQ=255)
                                                                                                                                                                   tGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTcc                                                                            >  1:597874/1‑250 (MQ=255)
                                                                                                                                                                   tGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTc                                                                             >  1:156759/1‑250 (MQ=255)
                                                                                                                                                                    gtagatgtctgtttaggttaTGGTTTTACAACCCTGCACGTTTATACGTGGATGATGGAAATCATTCGTCCCCCGTTCCAATTCGCGCTCGCGTATCAGCCCGTGACCGATGGCGATCCTCGCTTCTGGCGTCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCGCCATGCTGTCa                                                                            <  2:22132/231‑1 (MQ=255)
                                                                                                                                                                      cgatccctggagcGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTACAATTCGCGCTCGCGCATATGCCCGTGACAGATGGCGATACGCGCTTCTGGCACCAGTTCAGCCAGCCGATCGGCGGCTTTCTGCATGTTTTCCAAATCGTTGTAGAGATAATGACCCGGTCCTCCGAGAAAAATTTAACGCGGGATCGACTCCAGGACCACAAAGCTGGCATa                                                                          >  2:85023/13‑250 (MQ=255)
                                                                                                                                                                      cGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTGCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATa                                                                          <  1:497714/250‑1 (MQ=255)
                                                                                                                                                                        atcacGGTTTCGATACTGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCACTTCGCGCTCGCGCATCTGCCCGTGACCGATGACGATCCGCGCTTCTGGCACCAGCTCTGCCAGCCGTTCGGCGGCTTTCTGCATGGTTTCAACATCCTTGGAGAGATAATAAACCTGTCCTCCGCCCAAAATTTCACGCAGGATCGCCACCCGGACCACCATTCTGGCATCCTc                                                                       >  2:156852/5‑251 (MQ=255)
                                                                                                                                                                                   gATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAgg                                                            <  1:134155/251‑1 (MQ=255)
                                                                                                                                                                                    aTAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGt                                                           >  1:68665/1‑251 (MQ=255)
                                                                                                                                                                                       aTGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCcac                                                                                                                                                        >  1:88390/1‑155 (MQ=255)
                                                                                                                                                                                       aTGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCGCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCcac                                                                                                                                                        <  2:88390/155‑1 (MQ=255)
                                                                                                                                                                                         ggTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTaa                                                      <  1:219970/251‑1 (MQ=255)
                                                                                                                                                                                          gTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAAc                                                     <  1:332122/251‑1 (MQ=255)
                                                                                                                                                                                                cAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCACGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGcc                                                 >  2:141411/1‑249 (MQ=255)
                                                                                                                                                                                                  aaCCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTAATACTCACGGCCAAAGGTTTTAACTGCCAgc                                              >  2:219297/1‑249 (MQ=255)
                                                                                                                                                                                                   aCCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAgacg                                            <  1:595217/251‑1 (MQ=255)
                                                                                                                                                                                                     cAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAgacgac                                          <  1:294264/251‑1 (MQ=255)
                                                                                                                                                                                                       gCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGAcgg                                        >  1:22551/1‑251 (MQ=255)
                                                                                                                                                                                                          cGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGAcgggcg                                     >  1:582360/1‑251 (MQ=255)
                                                                                                                                                                                                               aaccGTTGATGATGGAAATCATTCATCTCCCGGTCCAATTCGCGCTCGCGCATATGACCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCAAGACGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGTACCACCATGCTGTCATACGCCCGGCAAAAGGTTTTAACGACCAGACGCcggcagggagg                                >  2:2284/4‑243 (MQ=255)
                                                                                                                                                                                                                          atgGAAATCATTCATCACCCGTGCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCGCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATaa                     <  2:342285/251‑1 (MQ=255)
                                                                                                                                                                                                                           tggcaaTCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTACGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAAt                    >  1:613137/5‑251 (MQ=255)
                                                                                                                                                                                                                            ggAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGCCAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATaagc                   >  1:442976/1‑249 (MQ=255)
                                                                                                                                                                                                                            ggAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAAt                    <  1:528074/250‑1 (MQ=255)
                                                                                                                                                                                                                                  cattcatCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGGTAATCggccgg             >  2:124355/1‑246 (MQ=255)
                                                                                                                                                                                                                                      catcaCCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAGGTCAc         >  1:301493/1‑251 (MQ=255)
                                                                                                                                                                                                                                       atcaCCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCCTGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAGGGCACg        >  2:513554/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aCCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAACATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAGGTCACGCAt     >  1:146718/1‑251 (MQ=255)
                                                                                                                                                                                                                                               gtCCAATTCGCGCGCGCGCATCTGCCCGTGTCCGATGGCGATCCGCGTTTCTGGCACCAGTTATGCCAGCCGTTCGGCGGCTTTCTGAATGTTTGCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGAGCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAGGTCACGCATTcc  <  2:91767/248‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
TCAAGACGTTTTTGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAGTTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAGGTCACGCATTCC  >  NC_000913/1171173‑1171660

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: