Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,193,286 G→C 8.8% D97E (GAC→GAG) ‡ thiF ← ThiS adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,193,2860GC 8.8% 252.4 / 12.9 80D97E (GAC→GAG) ‡thiFThiS adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (44/28);  new base C (4/3);  total (48/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

TGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGC  >  NC_000913/4193147‑4193431
                                                                                                                                           |                                                                                                                                                 
tgtTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCTCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCa                                                                                                                                         >  1:1243399/1‑150 (MQ=255)
  ttGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAAc                                                                                                                                       <  1:1595321/150‑1 (MQ=255)
   tGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACt                                                                                                                                      >  2:1528870/1‑150 (MQ=255)
    gAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACtg                                                                                                                                     <  1:646962/150‑1 (MQ=255)
    gAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATCCCGACATTCAACtg                                                                                                                                     <  1:863354/150‑1 (MQ=255)
       tGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTgtgt                                                                                                                                  >  1:211330/1‑150 (MQ=255)
        gCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGACAGTTGAATCCCGACATTCAACtgacag                                                                                                                                 <  2:951753/150‑5 (MQ=255)
         ccACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCa                                                                                                                                >  2:1380743/1‑150 (MQ=255)
         ccACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCa                                                                                                                                >  1:224810/1‑150 (MQ=255)
          cACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAg                                                                                                                               >  1:1501357/1‑150 (MQ=255)
           aCGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATCCCGACATTCAACtgacagcat                                                                                                                              <  2:550372/150‑8 (MQ=255)
                    gcATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTgg                                                                                                                     <  1:749150/150‑1 (MQ=255)
                            ctcCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTg                                                                                                             <  2:189175/150‑1 (MQ=255)
                            ctcCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTg                                                                                                             <  2:211330/150‑1 (MQ=255)
                                 ggCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGAtt                                                                                                        >  2:1161817/1‑150 (MQ=255)
                                 ggCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACACGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGCCCTGCGAtt                                                                                                        >  2:1394409/1‑150 (MQ=255)
                                   agAGTCGCCATATTGTAGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTATTAACGCCACACCCGTTACCCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAAGAGTGTCAGTCGCTGTTGGCTGACCTGCGAttta                                                                                                      >  2:700216/2‑149 (MQ=255)
                                    gAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCg                                                                                                     <  2:752821/150‑1 (MQ=255)
                                      gTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCggg                                                                                                   <  1:1383864/150‑1 (MQ=255)
                                       tCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCCCCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGCCCTGCGATTTCGGGc                                                                                                  >  2:1453784/1‑150 (MQ=255)
                                       tCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGc                                                                                                  <  2:1501357/150‑1 (MQ=255)
                                         gCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcga                                                                                                >  2:1401476/1‑150 (MQ=255)
                                           cATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcgatc                                                                                              >  2:249636/1‑150 (MQ=255)
                                             tatTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcgatcga                                                                                            >  2:377214/1‑150 (MQ=255)
                                             tatTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcgatcga                                                                                            >  1:500009/1‑150 (MQ=255)
                                               ttGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGata                                                                                          <  2:1263731/150‑1 (MQ=255)
                                                    ggTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTc                                                                                     >  1:1426955/1‑150 (MQ=255)
                                                       aCAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGt                                                                                  <  2:690493/150‑1 (MQ=255)
                                                         aGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAAGCGTTGTTGTAATGCTGTGAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTgg                                                                                <  1:1247752/150‑1 (MQ=255)
                                                         aGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATTTTGTTATTATACTGTGTCATTCGATTTTGTCTGAATTGGTATTTCGGGCGATAGATTTCTTTAGTgg                                                                                <  2:539184/150‑1 (MQ=255)
                                                                  ccacATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGa                                                                       >  1:297249/1‑150 (MQ=255)
                                                                  ccacATCGGCCCGTGCAACCGCATCTTTTAACCCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGa                                                                       >  1:1320577/1‑150 (MQ=255)
                                                                          gCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGt                                                               <  1:834091/150‑1 (MQ=255)
                                                                           cccGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGtt                                                              <  1:882563/150‑1 (MQ=255)
                                                                           cccGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGtt                                                              >  1:137745/1‑150 (MQ=255)
                                                                              gTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCa                                                           <  2:342520/150‑1 (MQ=255)
                                                                               tGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAg                                                          <  1:997451/150‑1 (MQ=255)
                                                                                      gCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCtt                                                   >  1:254176/1‑150 (MQ=255)
                                                                                      gCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCtt                                                   >  2:629550/1‑150 (MQ=255)
                                                                                       cATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTa                                                  >  2:976238/1‑150 (MQ=255)
                                                                                         tCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTaaa                                                >  1:909043/1‑150 (MQ=255)
                                                                                         tCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTaaa                                                <  1:1161817/150‑1 (MQ=255)
                                                                                          cTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAAt                                               >  2:1314749/1‑150 (MQ=255)
                                                                                           ttttAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATg                                              <  2:22411/150‑1 (MQ=255)
                                                                                            tttAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGc                                             <  1:1380743/150‑1 (MQ=255)
                                                                                            tttAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGc                                             >  1:487891/1‑150 (MQ=255)
                                                                                             ttAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGca                                            >  1:1040801/1‑150 (MQ=255)
                                                                                               aaCGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATg                                              <  1:1025548/146‑1 (MQ=255)
                                                                                               aaCGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATg                                              >  2:1025548/1‑146 (MQ=255)
                                                                                               aaCGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGcaca                                          >  2:1622053/1‑150 (MQ=255)
                                                                                               aaCGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGcaca                                          <  1:249636/150‑1 (MQ=255)
                                                                                                aCGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGcacaa                                         >  1:1494003/1‑149 (MQ=255)
                                                                                                    cTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACAtcgtc                                     >  1:1530847/1‑150 (MQ=255)
                                                                                                    cTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACAtcgtc                                     >  1:28119/1‑150 (MQ=255)
                                                                                                     tCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACAtcgtcg                                    <  2:1460675/150‑1 (MQ=255)
                                                                                                        cccGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGcac                                           >  1:808514/1‑140 (MQ=255)
                                                                                                        cccGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGcac                                           <  2:808514/140‑1 (MQ=255)
                                                                                                        cccGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCa                                 >  1:1185892/1‑150 (MQ=255)
                                                                                                        cccGGTAACCGGTGTTGTAATGCTGTCAGTTGAATGTCGGGAGTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCa                                 <  1:733197/150‑1 (MQ=255)
                                                                                                          cGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATc                               >  2:811178/1‑150 (MQ=255)
                                                                                                             tAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGc                            >  1:1621013/1‑150 (MQ=255)
                                                                                                             tAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGCATTTGTCGTTGCAGATTGGTTAAATGCACATCGTCGTCATCtcc                            >  1:447336/1‑148 (MQ=255)
                                                                                                              aaCCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGcc                           <  1:1528870/150‑1 (MQ=255)
                                                                                                               aCCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCa                          <  1:629550/150‑1 (MQ=255)
                                                                                                                 cGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGt                        >  1:1628916/1‑150 (MQ=255)
                                                                                                                       tgtAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGc                  >  2:1160711/1‑150 (MQ=255)
                                                                                                                       tgtAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGc                  >  2:1161621/1‑150 (MQ=255)
                                                                                                                        gtAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCg                 <  1:147872/150‑1 (MQ=255)
                                                                                                                        gtAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCg                 >  2:805566/1‑150 (MQ=255)
                                                                                                                         tAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTg                             <  1:857610/137‑1 (MQ=255)
                                                                                                                         tAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTg                             >  2:857610/1‑137 (MQ=255)
                                                                                                                               acacAGTTGAATCCCGACATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGAc          >  2:364029/4‑150 (MQ=255)
                                                                                                                                gTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACg         >  1:581457/1‑150 (MQ=255)
                                                                                                                                cacAGTTGAATCCCGACATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACg         >  2:343214/3‑150 (MQ=255)
                                                                                                                                  cAGTTGAATCCCGACATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGcc       >  2:511753/1‑150 (MQ=255)
                                                                                                                                  cAGTTGAATCCCGACATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCAGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGcc       >  1:1224425/1‑150 (MQ=255)
                                                                                                                                      tGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAgcg   >  1:977774/1‑150 (MQ=255)
                                                                                                                                       gAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAgcgc  >  2:903161/1‑150 (MQ=255)
                                                                                                                                          tGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACg         <  1:653953/140‑1 (MQ=255)
                                                                                                                                          tGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACg         >  2:653953/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
TGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGC  >  NC_000913/4193147‑4193431

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: ‑