Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 136,532 G→C 5.8% intergenic (‑68/+38) speE ← / ← yacC spermidine synthase (putrescine aminopropyltransferase)/PulS_OutS family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913136,5320GC 5.8% 292.4 / 4.3 95intergenic (‑68/+38)speE/yacCspermidine synthase (putrescine aminopropyltransferase)/PulS_OutS family protein
Reads supporting (aligned to +/- strand):  ref base G (56/33);  new base C (3/3);  total (59/36)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAGGGCT  >  NC_000913/136384‑136681
                                                                                                                                                    |                                                                                                                                                     
aTGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTGGGCAATGGGTTGTAACCACGTTGTTTAACACCCGAAAAAAAAGGGGCCAAAAAAATAGATAAATTTAACCCCggt                                                                                                                                                      >  1:1275597/1‑149 (MQ=255)
      cAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACACGTAAAAAAAGTGCGCAACATAATAGCTAACATTGACCGCGGATGCACg                                                                                                                                                >  2:1116777/1‑150 (MQ=255)
         aaCGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTc                                                                                                                                             <  2:931651/150‑1 (MQ=255)
         aaCGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTc                                                                                                                                             <  1:362154/150‑1 (MQ=255)
          aCGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCa                                                                                                                                            <  2:237612/150‑1 (MQ=255)
             ttATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGTGCATCCGCGGTCAatg                                                                                                                                         <  2:1216296/150‑2 (MQ=255)
              tATCTACCGCACAGTACTGCCGAAACTGGTCGTGTCGCGTTTCATGCCAATTTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGCCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAatat                                                                                                                                        <  2:260627/150‑1 (MQ=255)
              tATCTACCGCAAAGTAATGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGTGCATCCGCGGTCAatgt                                                                                                                                        <  1:501453/150‑3 (MQ=255)
               aTCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATAtt                                                                                                                                       <  2:811303/150‑1 (MQ=255)
               aTCTACCGCAAAGTACTGACCAAACTGGTCGTGTAGCGTTTCATTCCACAGTTTTTTTTCGGCCAAGGGTTGATACCTACTTTGTTAAAACCCGGAAAAAAAGGGCGCAACATAATAGCAAACATTTACAGCTGATGCACGGTCAataat                                                                                                                                       >  2:545198/1‑148 (MQ=255)
                 cTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTTTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAAAATTGACCGCGGATGCACGGTAAATAttta                                                                                                                                     >  1:1622400/1‑149 (MQ=255)
                 atACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGTTCAATATTTc                                                                                                                                     <  1:545198/149‑1 (MQ=255)
                  tACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCt                                                                                                                                    <  2:1059239/150‑1 (MQ=255)
                   gccGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTa                                                                                                                                   >  1:379789/2‑150 (MQ=255)
                   aCCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAATATTGACCGTGCATCCGCGGTCAATGttagct                                                                                                                                   <  2:1206818/150‑5 (MQ=255)
                   aCCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTa                                                                                                                                   >  2:649803/1‑150 (MQ=255)
                   aCCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTAGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTa                                                                                                                                   >  2:1288054/1‑150 (MQ=255)
                    ccGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTAc                                                                                                                                  >  1:862108/1‑150 (MQ=255)
                      gCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCCTGCCACTGTTTTTTTTCGGCCATGGGTTGATAACTCCTTTGTTAACACCCGTAAAAAAAGGGCGAAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACaa                                                                                                                                >  2:1282238/1‑150 (MQ=255)
                      gCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATTGCTAACATTGACCGCGGATGCACGGGCAATTTTTCTACaa                                                                                                                                >  2:1671610/1‑150 (MQ=255)
                        aaaGTACTGCCCAAACTGGTCGTGTAGCTTTTCATGCCACTGTTTTTTTTCGTCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCTGATGCACGGTCAATATTTCTACaaaa                                                                                                                              <  1:1269193/150‑1 (MQ=255)
                           ttacttaccAAACTGGTCGTGTATCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAgg                                                                                                                           <  2:1622400/143‑1 (MQ=255)
                             aCTGACCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTa                                                                                                                         <  1:1329651/150‑1 (MQ=255)
                              cTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTAt                                                                                                                        <  2:1407362/150‑1 (MQ=255)
                                gCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACAATTACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTAtaa                                                                                                                      >  2:593194/1‑148 (MQ=255)
                                 cccAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCATAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAg                                                                                                                     >  1:418142/1‑150 (MQ=255)
                                  ccAAACTGGTCGTGTAGCGTTTCATGCCACTGGTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAgg                                                                                                                    <  2:302881/150‑1 (MQ=255)
                                     aaCTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAATAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTTCAAAAAGGTATCAGGGat                                                                                                                 <  1:1516165/150‑1 (MQ=255)
                                       cTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGatta                                                                                                               <  2:1011431/150‑1 (MQ=255)
                                         ggTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGattatt                                                                                                             >  2:1362052/1‑150 (MQ=255)
                                         ggTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGattatt                                                                                                             >  2:516605/1‑150 (MQ=255)
                                         ggTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACAACCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGattatt                                                                                                             >  2:396867/1‑150 (MQ=255)
                                           tCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGTGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCTCGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTg                                                                                                           <  1:1528971/150‑1 (MQ=255)
                                           gcGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTGGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGATCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTg                                                                                                           <  2:446221/149‑1 (MQ=255)
                                             gtgtAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGAACGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGAc                                                                                                         >  1:1492273/1‑150 (MQ=255)
                                                   cGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGc                                                                                                   >  1:710260/1‑150 (MQ=255)
                                                         aTGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCCCGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAg                                                                                             >  2:422706/1‑150 (MQ=255)
                                                          tGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAgg                                                                                            <  2:465534/150‑1 (MQ=255)
                                                            ccACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCt                                                                                          <  1:1446134/150‑1 (MQ=255)
                                                            ccACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCt                                                                                          <  2:1653638/150‑1 (MQ=255)
                                                              aCTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTa                                                                                        >  2:780079/1‑150 (MQ=255)
                                                               cTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTaa                                                                                       >  1:855604/1‑150 (MQ=255)
                                                                  ttttttttCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGa                                                                                    >  2:1299160/1‑150 (MQ=255)
                                                                  ttttttttCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGGCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCCGGATTAAGGa                                                                                    >  2:852882/1‑150 (MQ=255)
                                                                     tttttCGGCCATGGGTTGATACCTCCTTTGTTACCACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCCCGGTCAATATTTCTACAAAAAGGTATCACGGATTATTTGACGTAGGCAATCAGGCTTAAGGAATc                                                                                 >  1:1469193/1‑150 (MQ=255)
                                                                     tttttCGGACATGGGTTGATACCTCCTTTGTTAACACCAGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATc                                                                                 <  1:99941/150‑1 (MQ=255)
                                                                      ttttCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATcg                                                                                >  2:629973/1‑150 (MQ=255)
                                                                       tttCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATcgc                                                                               <  1:649803/150‑1 (MQ=255)
                                                                       tttCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATcgc                                                                               >  1:743379/1‑150 (MQ=255)
                                                                         tCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCgg                                                                             >  1:1608488/1‑150 (MQ=255)
                                                                         tCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAACGTCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCgg                                                                             <  1:482139/150‑1 (MQ=255)
                                                                              cATGGGTTGATTCCTCCTTTGTTAACACCCGTAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGCCTTAAGGAATCGCGGGCCAgg                                                                       >  1:709363/1‑150 (MQ=255)
                                                                              cATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGCCCGCGGATGCACGGTCAATATTTCTACATAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGATTCGCGGGCCAg                                                                        >  1:822386/1‑150 (MQ=255)
                                                                                tGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAggg                                                                      <  1:1014552/150‑1 (MQ=255)
                                                                                 gggTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGTGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGGCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGc                                                                     <  1:1085754/150‑1 (MQ=255)
                                                                                  ggTTGATACCTCCTTTGTTAACCCCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCAGGGCtt                                                                   >  1:665140/1‑150 (MQ=255)
                                                                                   gTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCtt                                                                   >  2:1553574/1‑150 (MQ=255)
                                                                                     tGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGACAATATTTCTACAAAAGGTATAAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGATTTGc                                                                >  2:1649702/1‑150 (MQ=255)
                                                                                     tGATACCTCCTTTGATAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTg                                                                 <  2:1192472/150‑1 (MQ=255)
                                                                                     tGATACCCCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACAAAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATATGACGTAGGCAAGCAGGCTTAAGGACTCGCGGGCCAGGGCTATg                                                                 >  2:776374/1‑150 (MQ=255)
                                                                                      gATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGc                                                                <  2:862108/150‑1 (MQ=255)
                                                                                         aCCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCAtt                                                             >  2:978813/1‑150 (MQ=255)
                                                                                         aCCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCAtt                                                             >  2:191232/1‑150 (MQ=255)
                                                                                          cctcctTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCAttt                                                            >  2:1208967/1‑150 (MQ=255)
                                                                                             cctTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCAtttttt                                                         <  2:433488/150‑1 (MQ=255)
                                                                                             cctTTGTTAACACACGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCAtttttt                                                         <  1:676331/150‑1 (MQ=255)
                                                                                                 tGTTAACACCCGTAAAAAAAGGGCGCAACATGATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCg                                                     >  1:1587916/1‑150 (MQ=255)
                                                                                                   ttAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCgcc                                                   >  1:1549614/1‑148 (MQ=255)
                                                                                                    tAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACg                                                  >  2:200548/1‑150 (MQ=255)
                                                                                                     aaCACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACgg                                                 >  2:42179/1‑150 (MQ=255)
                                                                                                      acacCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACggg                                                >  1:1340223/1‑150 (MQ=255)
                                                                                                         cccGTAAAAAAAGGGCGCATCATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCACGGCTTTGCATTTTTTAGCGACGGGaat                                             >  2:1193713/1‑150 (MQ=255)
                                                                                                           cGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgc                                           <  2:1305020/150‑1 (MQ=255)
                                                                                                               aaaaaaGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgccaat                                       >  1:521780/1‑150 (MQ=255)
                                                                                                               aaaaaaGGGCGCAACATAATAGCTAACATTGACAGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgccaat                                       >  1:1595798/1‑150 (MQ=255)
                                                                                                                 aaaaGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGTCAAGCAGGCTTAAGGAATCGCTGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgccaatgc                                     >  1:256337/1‑150 (MQ=255)
                                                                                                                  aaaGGGCGTAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgccaatgcc                                    >  1:946182/1‑150 (MQ=255)
                                                                                                                  aaaGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgccaatgcc                                    <  1:629973/150‑1 (MQ=255)
                                                                                                                    aGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTGTGCATTTTTTAGCGACGGGAATGCCAATGCCGc                                  >  1:1233555/1‑150 (MQ=255)
                                                                                                                      ggCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGTCGCTg                                >  2:657439/1‑150 (MQ=255)
                                                                                                                              ataataGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATcgcg                        <  2:1275597/150‑1 (MQ=255)
                                                                                                                               aaaaaGCTAACATTCACCGCGGATGCACGGTCAATATTACTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCgg                       <  2:1119807/146‑1 (MQ=255)
                                                                                                                                aataGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGt                      >  2:1690153/1‑150 (MQ=255)
                                                                                                                                   aGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGTGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAga                   >  1:936677/1‑149 (MQ=255)
                                                                                                                                   aGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGc                   >  2:782481/1‑150 (MQ=255)
                                                                                                                                   aGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGc                   >  2:886629/1‑150 (MQ=255)
                                                                                                                                   aGCTAACATTGACCGCGGAAGCACCGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGACTCGCGGGCCAGGGCCTGGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGAGCGCGGTAGc                   <  2:1250829/150‑1 (MQ=255)
                                                                                                                                    gCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCTGGCTTATGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCAGCTGAGAACGCGGTAgat                  >  1:126591/1‑148 (MQ=255)
                                                                                                                                    agaaataTTGACCGTGCATCCGCGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCt                  >  1:1186015/7‑150 (MQ=255)
                                                                                                                                     cTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTg                 >  2:1158356/1‑150 (MQ=255)
                                                                                                                                     cTAACATTGACCGCGGATGCACGGTCAAGATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCATGCTGGCTTAAGGAATCGCGGGCCAGGGATTTGGAGTTTTTTGCGACGGGAATGCCATTGCGGGTGAGATCGGGGTAGCTg                 >  2:1543167/1‑150 (MQ=255)
                                                                                                                                      tAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGt                >  1:737776/1‑150 (MQ=255)
                                                                                                                                         cATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCtt             >  1:368161/1‑150 (MQ=255)
                                                                                                                                          aTTGACCGTGCATCCGCGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACTTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTAAATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGAGGGAGCTGTCtta            >  2:71677/1‑149 (MQ=255)
                                                                                                                                           ttGACCGTGCATCCGCGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCa           >  2:532229/1‑150 (MQ=255)
                                                                                                                                             gACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCAcc         >  1:663909/1‑150 (MQ=255)
                                                                                                                                                    gATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAGGGCt  >  2:478403/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAGGGCT  >  NC_000913/136384‑136681

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