Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 136,532 G→C 7.4% intergenic (‑68/+38) speE ← / ← yacC spermidine synthase/putative lipoprotein YacC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913136,5320GC 7.4% 267.7 / 9.9 81intergenic (‑68/+38)speE/yacCspermidine synthase/putative lipoprotein YacC
Reads supporting (aligned to +/- strand):  ref base G (42/33);  new base C (3/3);  total (45/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 = 4.90e-01

GATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAGGG  >  NC_000913/136386‑136679
                                                                                                                                                  |                                                                                                                                                   
gATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATg                                                                                                                                                  <  2:941545/150‑1 (MQ=255)
gATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATg                                                                                                                                                  >  1:1628228/1‑150 (MQ=255)
  tACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCa                                                                                                                                                >  1:289216/1‑150 (MQ=255)
  tACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCAAACATTGACCGGGGGTGcc                                                                                                                                                >  2:1128800/1‑149 (MQ=255)
    cAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACg                                                                                                                                              >  1:611521/1‑150 (MQ=255)
    cAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACg                                                                                                                                              >  1:267694/1‑150 (MQ=255)
        aCGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTACCACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCa                                                                                                                                          >  2:1263082/1‑150 (MQ=255)
         cGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCaa                                                                                                                                         >  1:1181665/1‑150 (MQ=255)
           ttATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGTGCATCCGCGGTCAatg                                                                                                                                       <  2:284809/150‑2 (MQ=255)
             aTCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATAtt                                                                                                                                     <  1:794850/150‑1 (MQ=255)
             aTCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATAtt                                                                                                                                     >  2:706633/1‑150 (MQ=255)
               cTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGTGCATCCGCGGTCAATGttag                                                                                                                                   <  1:82242/150‑3 (MQ=255)
                 aCCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTa                                                                                                                                 <  1:1443704/150‑1 (MQ=255)
                 aCCGCAAAGTACTGCCCAAACTGGTCATGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGTGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACAGCGGATGCACGGTCAATATTTCTa                                                                                                                                 >  2:855961/1‑150 (MQ=255)
                       aaGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACaaaaa                                                                                                                           <  1:1525643/150‑1 (MQ=255)
                         gTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGTACGGTCAATATTTCTACAAAAAgg                                                                                                                         <  1:1424817/150‑1 (MQ=255)
                         gTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAgg                                                                                                                         <  2:1168615/150‑1 (MQ=255)
                         gTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGACATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAgg                                                                                                                         >  2:156553/1‑150 (MQ=255)
                          tACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGt                                                                                                                        >  1:664967/1‑150 (MQ=255)
                          tACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGt                                                                                                                        <  1:657971/150‑1 (MQ=255)
                             tGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATc                                                                                                                     >  1:605754/1‑150 (MQ=255)
                              gCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCa                                                                                                                    >  2:1227721/1‑150 (MQ=255)
                                   aaCTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTGGGACATGGGTTTATACCTGCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCTCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGat                                                                                                               <  2:603471/150‑1 (MQ=255)
                                   aaCTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGat                                                                                                               >  1:835990/1‑150 (MQ=255)
                                   aaCTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGat                                                                                                               >  1:830306/1‑150 (MQ=255)
                                         tCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTg                                                                                                         <  1:1369838/150‑1 (MQ=255)
                                          cGTGTAGGGTTTCATGCCACTGTTTTTTTTCGGGCCTTGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGTGCGCAACATAATAGCTAACATTGACCGCGGATGCACTGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGa                                                                                                        <  1:405432/150‑1 (MQ=255)
                                           gtgtAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGAc                                                                                                       <  2:1390639/150‑1 (MQ=255)
                                                 cGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGc                                                                                                 <  2:1013727/150‑1 (MQ=255)
                                                  gTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCa                                                                                                <  2:516777/150‑1 (MQ=255)
                                                      cATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCa                                                                                            >  1:538238/1‑150 (MQ=255)
                                                          ccACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCt                                                                                        <  1:1321929/150‑1 (MQ=255)
                                                             cTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTaa                                                                                     <  1:1227721/150‑1 (MQ=255)
                                                               gTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAgg                                                                                   <  2:289216/150‑1 (MQ=255)
                                                                 tttttttCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCCGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGCaa                                                                                 >  2:289387/1‑150 (MQ=255)
                                                                          gCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATcgc                                                                             >  2:1059988/1‑145 (MQ=255)
                                                                          gCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATcgc                                                                             <  1:1059988/145‑1 (MQ=255)
                                                                            cATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACAGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAg                                                                      >  2:1037928/1‑150 (MQ=255)
                                                                             aTGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAgg                                                                     >  2:679570/1‑150 (MQ=255)
                                                                              tGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAggg                                                                    >  1:38976/1‑150 (MQ=255)
                                                                                     aTACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACTTTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCa                                                             <  1:37559/150‑1 (MQ=255)
                                                                                      tACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCAt                                                            >  2:973464/1‑150 (MQ=255)
                                                                                      tACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCAt                                                            >  2:672932/1‑150 (MQ=255)
                                                                                        cctcctTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGGCGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCAttt                                                          >  1:182840/1‑150 (MQ=255)
                                                                                           cctTTGTTAACACCCGTAAAAAAAGGGCGCAACATAGTAGCTAACATTGGCCGCGGATGCACGGTCAATATTTCGACTAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTGATGCCATCGCCGGCCAGGGCTTTGCCttgttt                                                       >  2:1444684/1‑150 (MQ=255)
                                                                                           cctTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCAtttttt                                                       >  2:1075281/1‑150 (MQ=255)
                                                                                           cctTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCAtttttt                                                       >  2:1406632/1‑150 (MQ=255)
                                                                                            ctTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTa                                                      >  2:935455/1‑150 (MQ=255)
                                                                                             tttGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAg                                                     <  2:664967/150‑1 (MQ=255)
                                                                                                 ttAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGAc                                                 >  1:1585443/1‑150 (MQ=255)
                                                                                                   aaCACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACgg                                               <  1:987882/150‑1 (MQ=255)
                                                                                                       cccGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaat                                           >  1:1184275/1‑150 (MQ=255)
                                                                                                       cccGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaat                                           >  2:358206/1‑150 (MQ=255)
                                                                                                       cccGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaat                                           >  2:1441634/1‑150 (MQ=255)
                                                                                                         cGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgc                                         >  1:1476743/1‑150 (MQ=255)
                                                                                                          gTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgcc                                        <  1:672932/150‑1 (MQ=255)
                                                                                                          gTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgcc                                        <  1:1472194/150‑1 (MQ=255)
                                                                                                             aaaaaaGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgccaat                                     >  2:832268/1‑150 (MQ=255)
                                                                                                               aaaaGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgccaatgc                                   >  1:723113/1‑150 (MQ=255)
                                                                                                                    ggCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTg                              >  1:1618205/1‑150 (MQ=255)
                                                                                                                    ggCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTg                              >  1:19594/1‑150 (MQ=255)
                                                                                                                    ggCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTg                              <  2:19594/150‑1 (MQ=255)
                                                                                                                      cgcAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTgag                            <  2:1098267/150‑1 (MQ=255)
                                                                                                                        cAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGAt                          >  2:1169902/1‑150 (MQ=255)
                                                                                                                              aataGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGt                    <  2:538238/150‑1 (MQ=255)
                                                                                                                                 tagaaataTTGACCGTGCATCCGCGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGc                 >  1:7310/8‑150 (MQ=255)
                                                                                                                                 tagaaataTTGACCGTGCATCCGCGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGc                 <  2:1319608/143‑1 (MQ=255)
                                                                                                                                 aGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGc                 <  1:1263082/150‑1 (MQ=255)
                                                                                                                                  gCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGAGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCt                >  2:103704/1‑150 (MQ=255)
                                                                                                                                  agaaataTTGACCGTGCATCCGCGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCt                >  2:943535/7‑150 (MQ=255)
                                                                                                                                   cTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTg               <  2:550162/150‑1 (MQ=255)
                                                                                                                                     aaCATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTc             <  1:1128800/150‑1 (MQ=255)
                                                                                                                                       cATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCtt           >  2:1055824/1‑150 (MQ=255)
                                                                                                                                        aTTGACCGTGCATCCGCGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTACCGACGCGAATGCCAATGCCGCTGAGATCGCGGTAGCTGACtta          >  2:565099/1‑149 (MQ=255)
                                                                                                                                        aTTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTc          <  1:156553/150‑1 (MQ=255)
                                                                                                                                          tGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCAc        <  1:1005951/150‑1 (MQ=255)
                                                                                                                                           gACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCgg                     <  1:296409/136‑1 (MQ=255)
                                                                                                                                           gACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCgg                     >  2:296409/1‑136 (MQ=255)
                                                                                                                                               gcgGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGTAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAgg   >  2:1247025/1‑150 (MQ=255)
                                                                                                                                               gcgGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCTGGCGAGGGCTTTTGATTTTTTAGCGACGGGAATGCGAATGCGTCTGAGATCGCGGTAGCTGTCTTCACCGAgg   <  2:491639/150‑1 (MQ=255)
                                                                                                                                                cgGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAggg  <  1:1317083/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAGGG  >  NC_000913/136386‑136679

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