Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,316,740 A→G 7.8% S263P (TCT→CCT)  nusA ← transcription termination/antitermination L factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,316,7400AG 7.8% 243.8 / 9.0 92S263P (TCT→CCT) nusAtranscription termination/antitermination L factor
Reads supporting (aligned to +/- strand):  ref base A (34/46);  new base G (4/3);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-01

CCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGC  >  NC_000913/3316592‑3316888
                                                                                                                                                    |                                                                                                                                                    
ccGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGac                                                                                                                                                     <  2:843904/150‑1 (MQ=255)
 cGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGacc                                                                                                                                                    >  2:1722462/1‑149 (MQ=255)
   gCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATATGCTCGCCACCCAGTTCAGTAGCCAcc                                                                                                                                                  >  2:1111776/1‑150 (MQ=255)
   gCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACAcc                                                                                                                                                  <  1:1267696/150‑1 (MQ=255)
     ttCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATAACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACCCTCGCCCCCCAGTTCAGTAGACACCGc                                                                                                                                                >  1:1567765/1‑150 (MQ=255)
       cAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGGCACCGCCt                                                                                                                                              >  1:470462/1‑150 (MQ=255)
        aaCGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTg                                                                                                                                             >  2:1667949/1‑150 (MQ=255)
          cGGCGATATCCATGGAGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGaa                                                                                                                                           >  1:1476081/1‑150 (MQ=255)
              gATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACAcg                                                                                                                                       <  2:1046472/150‑1 (MQ=255)
                tatCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACAcgcg                                                                                                                                     <  1:137279/150‑1 (MQ=255)
                       ggTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTCCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCACGCCACCCCGTTCAGTAGGCACCGCCTGCACCCGCGCGCCACg                                                                                                                              >  1:1609688/1‑150 (MQ=255)
                                 tCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTGCAGTAGACACCGCCTGAACACGCGCGCCCCGCATACCTACg                                                                                                                    >  1:620542/1‑150 (MQ=255)
                                 tCTTCATCAACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACTTACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCTGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACg                                                                                                                    <  1:746656/150‑1 (MQ=255)
                                     cATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAg                                                                                                                >  2:1642436/1‑150 (MQ=255)
                                     cATCCACAACGATAGAAGCAACGTCTGCCGGTGCCATTGAGTTAATCACGAACTGCGACGGGTTATCTTCCCACAGGACGATATCGATACGCTCGCCACCCAGTTAAGTAGAAACCGCCTGAACACGCGAGCCACGCATACATACGAAAg                                                                                                                >  2:951035/1‑150 (MQ=255)
                                      gtCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGc                                                                                                               >  2:1449050/2‑150 (MQ=255)
                                      aTCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTTCGTTAATCACTAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGc                                                                                                               <  2:1623830/150‑1 (MQ=255)
                                        ccaccaCGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCAc                                                                                                             >  2:1216765/1‑150 (MQ=255)
                                         caccacGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGGCACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCAcc                                                                                                            <  1:1112722/150‑1 (MQ=255)
                                          accacGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCt                                                                                                           <  1:1019887/150‑1 (MQ=255)
                                            cacGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTAc                                                                                                         <  2:1296860/150‑1 (MQ=255)
                                             acGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTAcc                                                                                                        <  2:1578494/150‑1 (MQ=255)
                                                aTAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGa                                                                                                     >  1:705244/1‑150 (MQ=255)
                                                 tAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGAt                                                                                                    >  1:591912/1‑150 (MQ=255)
                                                  aGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATc                                                                                                   <  2:1759633/150‑1 (MQ=255)
                                                   gAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCg                                                                                                  <  1:1560728/150‑1 (MQ=255)
                                                   gAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCg                                                                                                  >  1:769781/1‑150 (MQ=255)
                                                   gAAGCAACGTCTGCCGGTGCCATTGCGATAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCg                                                                                                  >  2:1548609/1‑150 (MQ=255)
                                                   gAAGCAACGTCTGCCGGTGCCATCGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCg                                                                                                  >  2:1570118/1‑150 (MQ=255)
                                                       cAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCCGGTCAGTAGCCACCGCCTGAACACGCGCGCCCCGCATACCTACGCCAGCCCCTTCCGGAACGAtcc                                                                                              >  2:296795/1‑148 (MQ=255)
                                                       cAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGGCACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATAc                                                                                              <  1:136305/150‑1 (MQ=255)
                                                        aaCGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACg                                                                                             <  1:1581266/150‑1 (MQ=255)
                                                          cGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGGCACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGtt                                                                                           >  2:1396757/1‑150 (MQ=255)
                                                          cGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACTCAAGCACCTACCGGATCGATACGtt                                                                                           <  2:1389107/150‑1 (MQ=255)
                                                             cTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTAt                                                                                        <  1:1352304/150‑1 (MQ=255)
                                                             cTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTAt                                                                                        <  1:1359724/150‑1 (MQ=255)
                                                              tGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCATGCACCTACCGGATCGATACGTTTATc                                                                                       <  1:1732880/150‑1 (MQ=255)
                                                                 cGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCTGATCGATACGTTTATCGtt                                                                                    <  1:1224343/150‑1 (MQ=255)
                                                                    tGCCATTGCGTTAATCACGAACTGCGCCGGGGTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGt                                                                                 <  1:670651/150‑1 (MQ=255)
                                                                       cATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGtttt                                                                              >  2:913977/1‑150 (MQ=255)
                                                                        aTTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTc                                                                             <  1:1667949/150‑1 (MQ=255)
                                                                         ttGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCa                                                                            <  1:654627/150‑1 (MQ=255)
                                                                            cGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGGCACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACcg                                                                         <  1:128340/150‑1 (MQ=255)
                                                                               tAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACcgcgt                                                                      >  2:1066628/1‑149 (MQ=255)
                                                                                aaTCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAt                                                                     >  2:753304/1‑150 (MQ=255)
                                                                                 aTCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAtt                                                                    >  1:1531640/1‑150 (MQ=255)
                                                                                 aTCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAtt                                                                    >  1:1634877/1‑150 (MQ=255)
                                                                                 aTCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCAGATCGATACGTTTATCGTTGGTTTTCACCGCGAtt                                                                    >  2:1582836/1‑150 (MQ=255)
                                                                                  tCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAttt                                                                   >  2:393290/1‑150 (MQ=255)
                                                                                    aCGAACTGCGCCGGGTTATCATCGCACAGCACGATATCGATACGCTCGCCACCCAGTTCCGTAGCCACCGCCTGAACACGCGCGCCACGCATACCTACGCCAGCACCTACCGGAGCGATACGTGTATCGTTGGTTTTCACCGCGATTTTc                                                                 <  1:1816522/150‑1 (MQ=255)
                                                                                    aCGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGCCACCGCCTGAACACGCGCGCCACGCATACCTACGCCAGCACCTACCGGATCGATACGTGTATCGTTGGTTTTCCCCGCGATTTTc                                                                 >  2:641449/1‑150 (MQ=255)
                                                                                     cGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCg                                                                <  1:1096882/150‑1 (MQ=255)
                                                                                     cGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCg                                                                <  2:1529633/150‑1 (MQ=255)
                                                                                      gAACTGCGCCGGGTTATCATCGAACATTACGATATCGATACGCTCGCCACCCAGTTCAGTAGACTCCGCCTGAACACGCGCGTCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTTTGGTTGGGTTTCACCGCGTTTTTCGc                                                               <  2:188150/150‑1 (MQ=255)
                                                                                      gAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAATCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGc                                                               <  1:1169381/150‑1 (MQ=255)
                                                                                      gAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGc                                                               <  1:1769312/150‑1 (MQ=255)
                                                                                      gAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGc                                                               <  1:1776438/150‑1 (MQ=255)
                                                                                        aCTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGccc                                                             >  1:171826/1‑148 (MQ=255)
                                                                                         cTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACg                                                            <  1:998788/150‑1 (MQ=255)
                                                                                           gcgcCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACgag                                                          >  2:487573/1‑150 (MQ=255)
                                                                                            cgcCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACgaga                                                         >  2:38588/1‑150 (MQ=255)
                                                                                                  gTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTTTCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCgg                                                   <  2:1507587/150‑1 (MQ=255)
                                                                                                   ttATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTGTTCACCGCGATGTTCGCACGAGAACCCGGa                                                  >  2:1294730/1‑150 (MQ=255)
                                                                                                         tcCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGa                                            <  2:1451019/150‑1 (MQ=255)
                                                                                                           ccACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGc                                          <  2:1421328/150‑1 (MQ=255)
                                                                                                           ccACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCAACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGc                                          <  1:1310334/150‑1 (MQ=255)
                                                                                                            cacaGGAAGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGcc                                         <  2:58479/150‑1 (MQ=255)
                                                                                                               aGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCgct                                      <  2:460541/150‑1 (MQ=255)
                                                                                                                  aCGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCgctgct                                   <  2:850160/150‑1 (MQ=255)
                                                                                                                       atCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAAt                              <  1:1216765/150‑1 (MQ=255)
                                                                                                                       atCGATACGCGCGCCACCCAGTTCCGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCCAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAAt                              <  1:304011/150‑1 (MQ=255)
                                                                                                                        tCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAAtt                             >  1:1669755/1‑150 (MQ=255)
                                                                                                                         cGATACGCTCGCCACCCAGTTCAGAAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAAttt                            <  1:145794/150‑1 (MQ=255)
                                                                                                                          gATACGCTCGCCACCCAGTTCAGTAGCCACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGAAACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGCGAACCCGGATCGCGAGCCGCTGGTTTAATATc                           >  1:1615344/1‑150 (MQ=255)
                                                                                                                          gATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTc                           <  1:881146/150‑1 (MQ=255)
                                                                                                                          gATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTc                           <  1:786624/150‑1 (MQ=255)
                                                                                                                            tACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCaa                         <  1:1370215/150‑1 (MQ=255)
                                                                                                                              cGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATc                       <  2:1500690/150‑1 (MQ=255)
                                                                                                                              cGCTCGCCACCCAGATCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATAGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCGGGATCGCGAGCCGCTGCTTTAATTTCAATc                       <  2:1567765/150‑1 (MQ=255)
                                                                                                                               gCTCGACACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGCTACGTTTTTCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATcc                      >  1:886836/1‑149 (MQ=255)
                                                                                                                                 tCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCAct                    >  1:1022150/1‑150 (MQ=255)
                                                                                                                                   gCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCActtc                  >  2:1506392/1‑150 (MQ=255)
                                                                                                                                    ccACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCActtct                 >  1:1210686/1‑150 (MQ=255)
                                                                                                                                    ccACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACAAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCActtct                 >  1:1662907/1‑150 (MQ=255)
                                                                                                                                      aCCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTCCGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCCCCGCGATTTTCGCACGAGAACCCGGATCGCGCGCCGCTGCTTTAATTTCAATCACttcttc               >  2:1425893/1‑150 (MQ=255)
                                                                                                                                       cccAGTTCAGTAGGCACCGCCTGAACACGCGCGCCACGCATACCTACGCCAGCACCTACCGGATCGAAACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCGTTAATTTCAATCACTTCTTCg              >  2:236960/1‑150 (MQ=255)
                                                                                                                                       cccAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCg              >  2:1673318/1‑150 (MQ=255)
                                                                                                                                         cAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGcc            <  1:284342/150‑1 (MQ=255)
                                                                                                                                          aGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCg           <  2:1751583/150‑1 (MQ=255)
                                                                                                                                                   gACACCGCCTGAACACGCGCGCCACTCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGc  <  1:327148/150‑1 (MQ=255)
                                                                                                                                                   gACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGCTACGTTTATCGTTGGTTTTCACCGCGATTTTCGCCCGAGCACCCGGATCGCGCGCCGCTGCTGTAATTTCAATCACTTCTTCGCCGATGTCTGGc  >  1:1727830/1‑150 (MQ=255)
                                                                                                                                                   gACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGc  <  1:1722462/150‑1 (MQ=255)
                                                                                                                                                   gACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGc  <  2:1531640/150‑1 (MQ=255)
                                                                                                                                                   gACACCGCCTGAACACGCGCGCCACGCATAACTACGCAAGCACCTACCGGATCGCTACGTTTATCGTTGGTTTTCACCGCGATTTACGCACGAGACCCCGGATCGCGAGCCGCTGCTTTAAGTTCAATCAATTCTTCGCCGATTTCTGGc  >  2:1171446/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGC  >  NC_000913/3316592‑3316888

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: