Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,906,275 A→T 7.1% L98Q (CTA→CAA)  bglG ← transcriptional antiterminator BglG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,906,2750AT 7.1% 259.4 / 5.6 82L98Q (CTA→CAA) bglGtranscriptional antiterminator BglG
Reads supporting (aligned to +/- strand):  ref base A (40/36);  new base T (2/4);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACGCCCAACCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACTGCTCAAGGCA  >  NC_000913/3906126‑3906418
                                                                                                                                                     |                                                                                                                                               
gCACGCCCAACCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTa                                                                                                                                                 >  1:531220/1‑150 (MQ=255)
   cGCCCAACCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGTATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAATCGTTTAATCGCTAACTGACCATTGCCAGTTTGCg                                                                                                                                              <  1:1577894/150‑1 (MQ=255)
    gCCCAACCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCga                                                                                                                                             <  2:862115/150‑1 (MQ=255)
     cccAACCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCgag                                                                                                                                            >  2:1332312/1‑150 (MQ=255)
      ccAACCGTTTAACAATAATGGTTAACGCTTCTTCCCATAGCTTCAACTCATTCGCGTAAAGCCGCTCGATATCCCACACCAACGGCTTGGGCAGCAACACCTTTAGCTGAAAGCGTTTAATCCCAAAATGGCAATGGTCAGTTAGCgaga                                                                                                                                           >  2:536091/1‑150 (MQ=255)
         aCCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCTAACTGACCATTGCCAGTTTGCgattaaa                                                                                                                                        <  1:346984/150‑6 (MQ=255)
         aCCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGatat                                                                                                                                        >  1:744721/1‑150 (MQ=255)
          ccGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGAtata                                                                                                                                       >  1:1358854/1‑150 (MQ=255)
           cGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCTAACTGACCATTGCCAGTTTGCgattaaacg                                                                                                                                      <  1:139915/150‑8 (MQ=255)
           cGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATaa                                                                                                                                      >  1:324061/1‑150 (MQ=255)
            gTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATaaa                                                                                                                                     >  2:852281/1‑150 (MQ=255)
                aTCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGCTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACt                                                                                                                                 >  1:1382349/1‑150 (MQ=255)
                   aataatGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTc                                                                                                                              >  2:1339815/1‑150 (MQ=255)
                       atGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGt                                                                                                                          <  2:1231924/150‑1 (MQ=255)
                         ggTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTaa                                                                                                                        >  1:1294563/1‑150 (MQ=255)
                            tAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAAttt                                                                                                                     >  2:871303/1‑150 (MQ=255)
                              aCGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTc                                                                                                                   <  2:425415/150‑1 (MQ=255)
                                gCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTccc                                                                                                                 <  2:92928/150‑1 (MQ=255)
                                 cttcttCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCATCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCa                                                                                                                <  2:1412602/150‑1 (MQ=255)
                                   tcttcCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAg                                                                                                              <  1:1373683/150‑1 (MQ=255)
                                     ttcCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAgcg                                                                                                            <  2:820596/150‑1 (MQ=255)
                                     ttcCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAgcg                                                                                                            <  1:247487/150‑1 (MQ=255)
                                       ccccTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCt                                                                                                          >  2:1207578/1‑150 (MQ=255)
                                       ccccTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCt                                                                                                          >  2:1289694/1‑150 (MQ=255)
                                         ccTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTcc                                                                                                        >  2:111610/1‑150 (MQ=255)
                                         ccTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTcc                                                                                                        >  2:365189/1‑150 (MQ=255)
                                           tAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTg                                                                                                      <  1:764149/150‑1 (MQ=255)
                                                   aCTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAag                                                                                              <  1:794663/150‑1 (MQ=255)
                                                     tctTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAagag                                                                                            <  1:1575273/150‑1 (MQ=255)
                                                     tctTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAagag                                                                                            <  2:1294563/150‑1 (MQ=255)
                                                      ctTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAgaga                                                                                           <  1:852281/150‑1 (MQ=255)
                                                      ctTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAgaga                                                                                           >  1:1099866/1‑150 (MQ=255)
                                                      ctTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAgaga                                                                                           >  2:1569452/1‑150 (MQ=255)
                                                          cGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGataat                                                                                       >  1:387012/1‑150 (MQ=255)
                                                          cGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATATTCCCAAGCGCTCCTGCGCTAAAGAGgtaat                                                                                       <  2:1437035/150‑1 (MQ=255)
                                                           gggTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGataata                                                                                      >  2:1022165/1‑150 (MQ=255)
                                                           gggTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGataata                                                                                      >  1:1076742/1‑150 (MQ=255)
                                                              tAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCTCATAGTTGCAATGTTCTGTTAGCGAGATATAAATACTGTCCTGTAATTTTGCCAATCTCTCCTGATCTAAATTGataattagt                                                                                   <  2:511278/150‑5 (MQ=255)
                                                                  gCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATcac                                                                               <  2:324061/150‑1 (MQ=255)
                                                                      cTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGt                                                                           <  1:521769/150‑1 (MQ=255)
                                                                       tGGATATCCCACAGCAACGTGTTGGTCAGCAACACGTTTTGCTGAAAGCGTTTAATAGCAAACTGGCAATGTTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACATGtt                                                                          <  1:247382/150‑1 (MQ=255)
                                                                        ggATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTg                                                                         >  2:811673/1‑150 (MQ=255)
                                                                             tCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATc                                                                    >  2:1520832/1‑150 (MQ=255)
                                                                              cccACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGAt                                                                                  >  2:1145087/1‑135 (MQ=255)
                                                                              cccACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGAt                                                                                  <  1:1145087/135‑1 (MQ=255)
                                                                              cccACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCa                                                                   >  1:975678/1‑150 (MQ=255)
                                                                              cccACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCa                                                                   >  2:660129/1‑150 (MQ=255)
                                                                                cacaGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCAcc                                                                 <  2:744721/150‑1 (MQ=255)
                                                                                cacaGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCAcc                                                                 >  1:1281691/1‑150 (MQ=255)
                                                                                   aGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCa                                                              <  1:365189/150‑1 (MQ=255)
                                                                                    gCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCaa                                                             >  2:1032103/1‑150 (MQ=255)
                                                                                     cAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAag                                                            >  2:942339/1‑150 (MQ=255)
                                                                                          ggTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGaa                                                       >  1:169478/1‑150 (MQ=255)
                                                                                          ggTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGaa                                                       >  2:468197/1‑150 (MQ=255)
                                                                                            ttGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAata                                                     <  2:1494155/150‑1 (MQ=255)
                                                                                              gggCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATg                                                   <  1:1332312/150‑1 (MQ=255)
                                                                                              gggCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATg                                                   <  1:1514171/150‑1 (MQ=255)
                                                                                               ggCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGa                                                  >  1:1374546/1‑150 (MQ=255)
                                                                                                gcagcaACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGAc                                                 <  2:1045199/150‑1 (MQ=255)
                                                                                                          gTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTc                                       >  1:20561/1‑150 (MQ=255)
                                                                                                            tttGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGc                                     <  2:422578/150‑1 (MQ=255)
                                                                                                                 tGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAAt                                >  1:137588/1‑150 (MQ=255)
                                                                                                                  gAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATc                               >  1:1481753/1‑150 (MQ=255)
                                                                                                                  gAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTAGCTTAAta                               >  2:1580919/1‑149 (MQ=255)
                                                                                                                     aGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGcc                            <  2:1358854/150‑1 (MQ=255)
                                                                                                                     aGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGcc                            <  2:531220/150‑1 (MQ=255)
                                                                                                                      gCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGccc                           <  2:144952/150‑1 (MQ=255)
                                                                                                                       cGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCg                          <  1:441210/150‑1 (MQ=255)
                                                                                                                        gTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCgt                         >  2:401877/1‑150 (MQ=255)
                                                                                                                           atctCGCTAACTGACCATTGCCAGTTTGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCgttca                      >  1:2964/4‑150 (MQ=255)
                                                                                                                            tctCGCTAACTGACCATTGCCAGTTTGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCACCCgttcag                     >  1:339163/3‑150 (MQ=255)
                                                                                                                               cGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCgttcagttc                  >  2:316803/1‑150 (MQ=255)
                                                                                                                               cGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCgttcagttc                  <  1:111610/150‑1 (MQ=255)
                                                                                                                               cGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCgttcagttc                  >  2:8697/1‑150 (MQ=255)
                                                                                                                                 cAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCAt                <  1:1339815/150‑1 (MQ=255)
                                                                                                                                  aaaCTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATg               <  2:1281691/150‑1 (MQ=255)
                                                                                                                                  aaaCTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATg               <  1:660129/150‑1 (MQ=255)
                                                                                                                                     cTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACt            >  2:352560/1‑150 (MQ=255)
                                                                                                                                         cAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACTGCTc        >  1:305482/1‑150 (MQ=255)
                                                                                                                                             tgccAGTTTGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACTGCTCAAgg    <  2:2964/147‑1 (MQ=255)
                                                                                                                                             ggTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACTGCTCAAgg    >  1:1190916/1‑150 (MQ=255)
                                                                                                                                              gTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACTGCTCAAGGc   <  1:1207578/150‑1 (MQ=255)
                                                                                                                                               tCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACTGCTCAAGGCa  <  2:1076742/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GCACGCCCAACCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACTGCTCAAGGCA  >  NC_000913/3906126‑3906418

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