Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 136,537 A→G 6.2% intergenic (‑73/+33) 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,5370AG 6.2% 269.1 / 7.2 98intergenic (‑73/+33)speE/yacCspermidine synthase (putrescine aminopropyltransferase)/PulS_OutS family protein
Reads supporting (aligned to +/- strand):  ref base A (56/33);  new base G (3/3);  total (61/37)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

CAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAGGGCTTTC  >  NC_000913/136390‑136684
                                                                                                                                                   |                                                                                                                                                   
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                                                                                             <  2:1653638/150‑1 (MQ=255)
                                                      ccACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCt                                                                                             <  1:1446134/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:743379/1‑150 (MQ=255)
                                                                 tttCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATcgc                                                                                  <  1:649803/150‑1 (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:886629/1‑150 (MQ=255)
                                                                                                                             aGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGc                      >  2:782481/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)
                                                                                                                                                 gCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAGGGCTTTc  <  2:710260/150‑1 (MQ=255)
                                                                                                                                                 gCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAGGGCTTTc  >  1:1054114/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAGGGCTTTC  >  NC_000913/136390‑136684

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: