Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,278,673 T→A 10.2% C65* (TGT→TGA) ‡ ghxP → guanine/hypoxanthine transporter GhxP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,278,6730TA10.2% 261.7 / 17.9 98C65* (TGT→TGA) ‡ghxPguanine/hypoxanthine transporter GhxP
Reads supporting (aligned to +/- strand):  ref base T (43/44);  new base A (3/7);  total (47/51)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

TAAAATTTCACAACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTT  >  NC_000913/4278526‑4278822
                                                                                                                                                   |                                                                                                                                                     
tAAAATTTCACAACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCt                                                                                                                                                     <  1:88559/150‑1 (MQ=255)
 aaaaTTTCACAACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTg                                                                                                                                                    <  2:141422/150‑1 (MQ=255)
 aaaaTTTCACAACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTg                                                                                                                                                    <  1:1009039/150‑1 (MQ=255)
 aaaaTTTCACAACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTg                                                                                                                                                    <  1:371069/150‑1 (MQ=255)
      ttCACAACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGc                                                                                                                                               <  1:1462830/150‑1 (MQ=255)
       tCACAACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCAGACAGGTTGca                                                                                                                                              <  2:559137/150‑2 (MQ=255)
        cacaACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCg                                                                                                                                             <  1:42246/150‑1 (MQ=255)
         acaACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCgg                                                                                                                                            <  1:778641/150‑1 (MQ=255)
         acaACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCAGACAGGTTGcaac                                                                                                                                            <  2:72464/150‑4 (MQ=255)
         acaACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCAGACAGGTTGcaac                                                                                                                                            <  1:82744/150‑4 (MQ=255)
         acaACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCAGACAGGTTGcaac                                                                                                                                            <  1:638746/150‑4 (MQ=255)
         acaACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCAGACAGGTTGcaac                                                                                                                                            <  1:1092274/150‑4 (MQ=255)
         acaACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCAGACAGGTTGcaac                                                                                                                                            <  1:1179720/150‑4 (MQ=255)
          caACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACCACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAACGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGGCTGGTTGCCggc                                                                                                                                           >  2:1303914/1‑149 (MQ=255)
          caACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCAGACAGGTTGcaacg                                                                                                                                           <  2:1175949/150‑5 (MQ=255)
                ggAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGt                                                                                                                                     <  2:173558/150‑1 (MQ=255)
                  aaGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTc                                                                                                                                   >  1:345910/1‑150 (MQ=255)
                      aCTGTCCGTCAGGAAGTGGGTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGGCATCGGCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTGCTATc                                                                                                                               <  2:166259/150‑1 (MQ=255)
                       cTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCg                                                                                                                              <  2:751062/150‑1 (MQ=255)
                       cTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCg                                                                                                                              >  1:134234/1‑150 (MQ=255)
                        tGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGt                                                                                                                             >  2:1132929/1‑150 (MQ=255)
                          tCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGa                                                                                                                           <  2:803153/150‑1 (MQ=255)
                           ccGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCg                                                                                                                                       <  1:1104613/137‑1 (MQ=255)
                           ccGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCg                                                                                                                                       >  2:1104613/1‑137 (MQ=255)
                                 ggAAGTTGTTGACGAGTTAAAAACTTTTATAGCGATTGTATAATCGGTCAACCACTTTCAAGATATATTTAGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCt                                                                                                                    >  2:1202560/1‑150 (MQ=255)
                                 ggAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCCTCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCt                                                                                                                    >  2:1422728/1‑150 (MQ=255)
                                    aGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCtgtg                                                                                                                 <  1:1217995/150‑1 (MQ=255)
                                     gTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGCAATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTgg                                                                                                                >  2:1396719/1‑150 (MQ=255)
                                      tGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTgg                                                                                                                <  2:771758/149‑1 (MQ=255)
                                      tGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTgg                                                                                                                >  1:771758/1‑149 (MQ=255)
                                       gggtGCCGGGTTAACAACGTTTCTGGCGATGTTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGc                                                                                                              <  2:1302013/147‑1 (MQ=255)
                                       ggTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGc                                                                                                              >  2:540978/1‑150 (MQ=255)
                                        gTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCt                                                                                                             >  2:1346556/1‑150 (MQ=255)
                                          tGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTaa                                                                                                           >  1:695344/1‑150 (MQ=255)
                                           gCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAAt                                                                                                          >  1:499537/1‑150 (MQ=255)
                                                 ttAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCGAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCg                                                                                                    <  2:1120787/150‑1 (MQ=255)
                                                   aacaacGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGtt                                                                                                  >  1:789131/1‑150 (MQ=255)
                                                   aacaacGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGtt                                                                                                  >  2:943949/1‑150 (MQ=255)
                                                    acaacGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTg                                                                                                 <  1:943949/150‑1 (MQ=255)
                                                       acGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCg                                                                                              <  2:464377/150‑1 (MQ=255)
                                                            tCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTgg                                                                                         <  2:750779/150‑1 (MQ=255)
                                                               ggCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGAAGCGACTGTGGGCTAATCTAACGTTTGCAattgattg                                                                                      >  2:792277/1‑150 (MQ=255)
                                                                   aTGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGcc                                                                                  <  1:1538605/150‑1 (MQ=255)
                                                                    tGGTCTACTCGGTCATCGGCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCa                                                                                 <  2:737147/150‑1 (MQ=255)
                                                                      gTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCAtc                                                                               >  1:1018223/1‑150 (MQ=255)
                                                                      gTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCAtc                                                                               <  2:134234/150‑1 (MQ=255)
                                                                       gcTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCAtct                                                                              <  1:1468037/149‑1 (MQ=255)
                                                                        cTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCAtctc                                                                             <  1:911257/150‑1 (MQ=255)
                                                                        cTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCAtctc                                                                             >  2:667020/1‑150 (MQ=255)
                                                                          tcTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCTCCATCTccc                                                                           <  2:42954/149‑1 (MQ=255)
                                                                          aCTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTccc                                                                           >  1:246465/1‑150 (MQ=255)
                                                                           cTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCt                                                                          <  1:167653/150‑1 (MQ=255)
                                                                           cTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCt                                                                          >  1:1189307/1‑150 (MQ=255)
                                                                           cTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCt                                                                          >  1:1194814/1‑150 (MQ=255)
                                                                           cTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCt                                                                          <  1:540978/150‑1 (MQ=255)
                                                                            tCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTg                                                                         >  1:190473/1‑150 (MQ=255)
                                                                                 cATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGc                                                                    >  2:953689/1‑150 (MQ=255)
                                                                                 cATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGc                                                                    >  2:1541123/1‑150 (MQ=255)
                                                                                  aTCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCg                                                                   >  1:801931/1‑150 (MQ=255)
                                                                                  aTCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCg                                                                   >  1:1183525/1‑150 (MQ=255)
                                                                                    cgtcgtTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGtt                                                                 <  1:84701/150‑1 (MQ=255)
                                                                                    cggcgtTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGGAGGTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGtt                                                                 <  2:68061/147‑1 (MQ=255)
                                                                                           ccAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCa                                                          <  2:802094/150‑1 (MQ=255)
                                                                                           ccAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCa                                                          >  2:929310/1‑150 (MQ=255)
                                                                                              ggTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCAtct                                                                              <  2:459938/127‑1 (MQ=255)
                                                                                              ggTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCAtct                                                                              >  1:459938/1‑127 (MQ=255)
                                                                                                tATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAg                                                     <  2:105223/150‑1 (MQ=255)
                                                                                                  tGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGcc                                                   >  2:885460/1‑150 (MQ=255)
                                                                                                   gTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCt                                                  >  1:1170798/1‑150 (MQ=255)
                                                                                                    ttGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTg                                                 <  2:1451400/150‑1 (MQ=255)
                                                                                                          aaaGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTg                                           <  1:667020/150‑1 (MQ=255)
                                                                                                           aaGCAGGCTTCCCGCCTGCGGGAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTgg                                          <  1:322473/150‑1 (MQ=255)
                                                                                                           aaGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTgg                                          >  2:17709/1‑150 (MQ=255)
                                                                                                           aaGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTgg                                          <  1:1202560/150‑1 (MQ=255)
                                                                                                            aGCAGGCTTCCCGCCTGCTGCAGTTTGCGTTGCATCCTGTCTGGTTGGCGGACTCGGTTCTATCGTGATGGGTCTGTGGTCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGGCTTGTGCTggg                                         <  2:612659/150‑1 (MQ=255)
                                                                                                            aGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTggg                                         <  2:1183525/150‑1 (MQ=255)
                                                                                                             gCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTgggg                                        >  2:1592708/1‑150 (MQ=255)
                                                                                                             gCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTgggg                                        <  2:345910/150‑1 (MQ=255)
                                                                                                                ggCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCaa                                     <  2:703129/150‑1 (MQ=255)
                                                                                                                 gCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAAc                                    <  2:499537/150‑1 (MQ=255)
                                                                                                                       cGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTa                              >  2:1241482/1‑150 (MQ=255)
                                                                                                                             cGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTAc                        >  2:990011/1‑150 (MQ=255)
                                                                                                                                cAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTg                     <  1:1132929/150‑1 (MQ=255)
                                                                                                                                    tttCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTcgc                 >  1:122148/1‑150 (MQ=255)
                                                                                                                                      tCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCgcgc               >  2:1123491/1‑150 (MQ=255)
                                                                                                                                      tCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCgcgc               >  2:1128945/1‑150 (MQ=255)
                                                                                                                                      gtccggCAACCAGACAGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCgcgc               >  1:1107543/6‑150 (MQ=255)
                                                                                                                                      gtccggCAACCAGACAGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCgcgc               >  1:1101951/6‑150 (MQ=255)
                                                                                                                                       tccggCAACCAGACAGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCt              >  1:744469/5‑150 (MQ=255)
                                                                                                                                        gTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTg             >  1:1038425/1‑150 (MQ=255)
                                                                                                                                        gTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTg             >  1:1076276/1‑150 (MQ=255)
                                                                                                                                          tGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTggg           <  1:1088344/150‑1 (MQ=255)
                                                                                                                                          tGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTggg           <  1:1422728/150‑1 (MQ=255)
                                                                                                                                          tGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTggg           >  2:314482/1‑150 (MQ=255)
                                                                                                                                           gCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGATAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGt          >  2:496107/1‑150 (MQ=255)
                                                                                                                                            cAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTg         >  1:1587800/1‑150 (MQ=255)
                                                                                                                                            cAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTg         >  1:1588044/1‑150 (MQ=255)
                                                                                                                                               ccTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCGGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCg      >  1:809091/1‑150 (MQ=255)
                                                                                                                                                   tCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGtt  <  2:1588044/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TAAAATTTCACAACGTGGAAGCACTGTCCGTCAGGAAGTGGTTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTT  >  NC_000913/4278526‑4278822

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