Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,303,611 T→A 5.7% I55N (ATT→AAT)  yhbW → putative luciferase‑like monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,303,6110TA 5.7% 334.1 / 5.2 115I55N (ATT→AAT) yhbWputative luciferase‑like monooxygenase
Reads supporting (aligned to +/- strand):  ref base T (40/68);  new base A (2/5);  total (42/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

ATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGC  >  NC_000913/3303463‑3303760
                                                                                                                                                    |                                                                                                                                                     
   gCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCATAACACCACAATATGACTGGCATTGcc                                                                                                                                                   <  2:1376798/150‑1 (MQ=255)
   gCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGcc                                                                                                                                                   <  1:1667389/150‑1 (MQ=255)
   gCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGcc                                                                                                                                                   <  2:1247587/150‑1 (MQ=255)
    cGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGCACACCACAATATGACTGGCATTGCCa                                                                                                                                                  >  1:1127111/1‑150 (MQ=255)
     gTTTTCGCTACTCGATCTGGCCCCCATTCCCGATGGTTCTTCAGCGCGAGAAGCATTCTACCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCGATCATCGCTACTGGCTGGCAGAACACACCAATATGACTGGCATTGCCAg                                                                                                                                                 >  1:1035970/1‑150 (MQ=255)
        ttCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCTCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAACCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGtgc                                                                                                                                              <  2:333003/150‑1 (MQ=255)
        ttCGCTACTCGATCGGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGACCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCAATGCCAGtca                                                                                                                                              <  2:1005686/150‑3 (MQ=255)
         tCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGtgctg                                                                                                                                            <  2:581171/150‑1 (MQ=255)
          cGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAACGCGCGGCTATCATCGCTACTGGCTGGCAGAACCGCTCAATATGACTGGCTTTGCCAGtgctg                                                                                                                                            >  1:204346/1‑150 (MQ=255)
           gCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGCACACCACAATATGACTGGCATTGCCAGtgctgc                                                                                                                                           >  2:1266545/1‑150 (MQ=255)
           gCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGACCACCACAATATGACTGGCATTGCCAGtgctgc                                                                                                                                           >  1:1157452/1‑150 (MQ=255)
           gCTACTCGATAGTGCACCCATTCCCGAAGGTTCTTCAGCGCGAGAAGAATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTCGCAGAACACCACAATCTGCCTGGCAATGCCAGtcatat                                                                                                                                           <  1:121676/150‑6 (MQ=255)
            cTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCAATGCCAGtcatatt                                                                                                                                          <  1:1692821/150‑7 (MQ=255)
             tACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCAATGCCAGtcatattg                                                                                                                                         <  1:722074/150‑8 (MQ=255)
              aCTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCAc                                                                                                                                        >  1:545985/1‑150 (MQ=255)
               cTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACg                                                                                                                                       <  2:1344170/150‑1 (MQ=255)
               cTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACg                                                                                                                                       <  2:1307443/150‑1 (MQ=255)
               cTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACg                                                                                                                                       >  2:901383/1‑150 (MQ=255)
                tCGATCTGGCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTc                                                                                                                                     >  1:1432128/1‑150 (MQ=255)
                 cGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTc                                                                                                                                     <  1:5862/150‑1 (MQ=255)
                 cGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTc                                                                                                                                     >  1:1030033/1‑150 (MQ=255)
                      tGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTAt                                                                                                                                <  2:43716/150‑1 (MQ=255)
                         cccccATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGa                                                                                                                             >  2:38154/1‑150 (MQ=255)
                           cccATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATc                                                                                                                           >  2:1371671/1‑150 (MQ=255)
                           cccATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATc                                                                                                                           >  1:1047771/1‑150 (MQ=255)
                           cccATTCCCGAAGGTTCTTCAGAGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATc                                                                                                                           <  2:1048033/150‑1 (MQ=255)
                            ccATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACTCCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCg                                                                                                                          <  2:910505/150‑1 (MQ=255)
                            ccATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACCCCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCg                                                                                                                          >  1:1504992/1‑150 (MQ=255)
                             cATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGGCGGTATTGATCgg                                                                                                                         <  1:933302/150‑1 (MQ=255)
                              aTTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGc                                                                                                                        <  1:1662970/150‑1 (MQ=255)
                                tCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTa                                                                                                                      <  2:1444620/150‑1 (MQ=255)
                                 cccGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTTCTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTAt                                                                                                                     <  1:1041043/150‑1 (MQ=255)
                                   cGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCCCAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTAGTGATCGGCTATCt                                                                                                                   >  1:155997/1‑150 (MQ=255)
                                   cGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCt                                                                                                                   >  1:1430431/1‑150 (MQ=255)
                                       ggTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCGTTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTggcg                                                                                                               >  1:1021873/1‑150 (MQ=255)
                                        gTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTggcgg                                                                                                              <  2:692683/150‑1 (MQ=255)
                                        gTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCGGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTggcgg                                                                                                              <  2:4279/150‑1 (MQ=255)
                                           cttcAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGa                                                                                                           <  2:1494806/150‑1 (MQ=255)
                                             tcAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAAt                                                                                                         <  1:1241553/150‑1 (MQ=255)
                                             tcAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAAt                                                                                                         <  1:283146/150‑1 (MQ=255)
                                              cAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATa                                                                                                        <  1:244778/150‑1 (MQ=255)
                                                gcgcGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATacc                                                                                                      <  1:367318/150‑1 (MQ=255)
                                                gcgcGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATacc                                                                                                      <  2:1186731/150‑1 (MQ=255)
                                                   cgAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATaccacc                                                                                                   >  1:1231916/1‑150 (MQ=255)
                                                     cgaAGCATGCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTTCAGGCTGCCAGACCACCACAATATGACTCGCATAGCCAGTGCTTCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACcaccac                                                                                                 <  1:213196/149‑1 (MQ=255)
                                                       aaGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGc                                                                                               <  1:1455955/150‑1 (MQ=255)
                                                         gCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTg                                                                                             <  1:1204474/150‑1 (MQ=255)
                                                           aTTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCa                                                                                           <  1:1008667/150‑1 (MQ=255)
                                                             tctcCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATc                                                                                         >  1:205678/1‑150 (MQ=255)
                                                               tcCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTg                                                                                       <  2:1131563/150‑1 (MQ=255)
                                                                cccACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTgg                                                                                      <  1:582453/150‑1 (MQ=255)
                                                                cccACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTgg                                                                                      <  2:772134/150‑1 (MQ=255)
                                                                  cACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTgggg                                                                                    <  2:151471/150‑1 (MQ=255)
                                                                  cACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTgggg                                                                                    <  2:1608126/150‑1 (MQ=255)
                                                                   aCTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGt                                                                                   >  1:454011/1‑150 (MQ=255)
                                                                   aCTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGt                                                                                   >  2:693208/1‑150 (MQ=255)
                                                                      ctctcGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTg                                                                                >  1:1367568/1‑150 (MQ=255)
                                                                        ctcGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTggc                                                                              >  2:697690/1‑150 (MQ=255)
                                                                             tctcGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGGCTGGCGGCGt                                                                         <  2:178946/150‑1 (MQ=255)
                                                                              ctcGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTg                                                                        <  2:1047771/150‑1 (MQ=255)
                                                                                  cccGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGt                                                                    >  2:1351915/1‑150 (MQ=255)
                                                                                    cGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTg                                                                  <  1:594356/150‑1 (MQ=255)
                                                                                     gTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGc                                                                 <  2:889542/150‑1 (MQ=255)
                                                                                           cTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAAcc                                                           <  1:901383/150‑1 (MQ=255)
                                                                                           cTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGGCTGGCGGCGTGATGTTGCCTAAcc                                                           <  2:1087844/150‑1 (MQ=255)
                                                                                           cTGAAAAGCGCGGCTATCATCGCTACTGGCAGGCAGAACACCACAATATGACTGGCATTGCCAGAGCTGCCACGTCGGAATTGATCGGCTATCTGGCGGCGAATACAAACACGCTGCATCAGGGGTCGGGCGGCGTGATGTAGCCTAAcc                                                           >  1:1615941/1‑150 (MQ=255)
                                                                                                 aGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCAc                                                     <  1:1589982/150‑1 (MQ=255)
                                                                                                    gcgGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGt                                                  <  2:813744/150‑1 (MQ=255)
                                                                                                        cTATCCTCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTCTCTGGCGGCGAATACCACCACGCTGCATCTCGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGt                                              >  1:260416/1‑150 (MQ=255)
                                                                                                          atcatcGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCa                                            <  1:38154/150‑1 (MQ=255)
                                                                                                              tcGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGc                                        <  1:663042/150‑1 (MQ=255)
                                                                                                               cGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCa                                       <  1:946722/150‑1 (MQ=255)
                                                                                                                gCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAg                                      <  1:1327197/150‑1 (MQ=255)
                                                                                                                gCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAg                                      <  1:1335602/150‑1 (MQ=255)
                                                                                                                 cTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGa                                     >  2:1399963/1‑150 (MQ=255)
                                                                                                                   aCTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAAc                                   <  2:660753/150‑1 (MQ=255)
                                                                                                                   aCTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAAc                                   <  1:1408725/150‑1 (MQ=255)
                                                                                                                     tggctggcAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAg                                 <  2:1348944/150‑1 (MQ=255)
                                                                                                                        ctggcAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGGTGGTCATTGCAGAACAGTTc                              <  2:1573379/150‑1 (MQ=255)
                                                                                                                          ggcAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCgg                            >  2:1427852/1‑150 (MQ=255)
                                                                                                                          ggcAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCAGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCgg                            >  2:1332578/1‑150 (MQ=255)
                                                                                                                           gcAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGc                           <  2:1465482/150‑1 (MQ=255)
                                                                                                                           gcAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGc                           <  2:247319/150‑1 (MQ=255)
                                                                                                                                acacCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCt                      <  2:1667754/150‑1 (MQ=255)
                                                                                                                                   ccacAATATGACTGGCATTGCCAGTGATGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTaa                   >  2:397149/1‑150 (MQ=255)
                                                                                                                                    cacaATATGACTTGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAAt                  <  2:155997/150‑1 (MQ=255)
                                                                                                                                    cacaATATGACTGGCATTGCCAGTGCTGCCCCGTCGGTATAGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTTCCTAACCACTCACCGTTTGTCATTGCAGAACAGTTCGGCACGCTTAAt                  >  2:256350/1‑150 (MQ=255)
                                                                                                                                    cacaATATGACTGGCATTGCCAGTGCGGCCACGTCGGTATTGATCGGATATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAAt                  <  1:155278/150‑1 (MQ=255)
                                                                                                                                     acaATATGACTCGCAATGCCAGTGCTGCCTAGTCGGTATTGTTCGGCTTTCTGGCGGCGAATACCACCACGCTGCATCTGTTGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGATCAGTTCGGCACGCTTAATa                 <  2:448715/150‑1 (MQ=255)
                                                                                                                                      caATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATac                >  2:1324424/1‑150 (MQ=255)
                                                                                                                                       aaTATGACTGGCATTGCCAGTGCTGCCACGGCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGTGTCTGGCGGCGTTATGTTGCCTAACCACTCACCGTTTGTCATTGCAGAACAGTTCGGCACGCTTAATaca               <  1:1344798/150‑1 (MQ=255)
                                                                                                                                        gcagcaCTGGCAATGCCAGTGCTGCCACGTCGGTATTGCACGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGAGAAGTTGCCTAACCACTCCCCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATacac              >  2:252048/6‑150 (MQ=255)
                                                                                                                                         tatGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACAct             <  2:347152/150‑1 (MQ=255)
                                                                                                                                         cagcaCTGGCAATGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACTCTTAATACAct             >  2:830505/5‑150 (MQ=255)
                                                                                                                                           tGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACActct           >  2:533493/1‑150 (MQ=255)
                                                                                                                                             aCTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTAt         <  1:196896/150‑1 (MQ=255)
                                                                                                                                              cTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATc        >  1:405028/1‑150 (MQ=255)
                                                                                                                                               tGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCGGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATcc       <  1:64536/150‑1 (MQ=255)
                                                                                                                                                  cATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCCCGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCggg    >  2:1335159/1‑150 (MQ=255)
                                                                                                                                                    ttGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATTCGGGGc  <  2:1699285/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGC  >  NC_000913/3303463‑3303760

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