Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,906,269 T→C 6.1% D100G (GAC→GGC)  bglG ← transcriptional antiterminator BglG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,906,2690TC 6.1% 244.9 / 5.0 82D100G (GAC→GGC) bglGtranscriptional antiterminator BglG
Reads supporting (aligned to +/- strand):  ref base T (42/35);  new base C (2/3);  total (44/38)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

TGCACGCCCAACCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACTGCTCAAGGCA  >  NC_000913/3906125‑3906418
                                                                                                                                                |                                                                                                                                                     
tGCACGCCCAACCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGtt                                                                                                                                                  >  1:545719/1‑150 (MQ=255)
 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:1289694/1‑150 (MQ=255)
                                        ccccTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCt                                                                                                          >  2:1207578/1‑150 (MQ=255)
                                          ccTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTcc                                                                                                        >  2:365189/1‑150 (MQ=255)
                                          ccTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTcc                                                                                                        >  2:111610/1‑150 (MQ=255)
                                            tAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTg                                                                                                      <  1:764149/150‑1 (MQ=255)
                                                    aCTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAag                                                                                              <  1:794663/150‑1 (MQ=255)
                                                      tctTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAagag                                                                                            <  2:1294563/150‑1 (MQ=255)
                                                      tctTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAagag                                                                                            <  1:1575273/150‑1 (MQ=255)
                                                       ctTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAgaga                                                                                           <  1:852281/150‑1 (MQ=255)
                                                       ctTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAgaga                                                                                           >  2:1569452/1‑150 (MQ=255)
                                                       ctTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAgaga                                                                                           >  1:1099866/1‑150 (MQ=255)
                                                           cGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGataat                                                                                       >  1:387012/1‑150 (MQ=255)
                                                           cGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATATTCCCAAGCGCTCCTGCGCTAAAGAGgtaat                                                                                       <  2:1437035/150‑1 (MQ=255)
                                                            gggTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGataata                                                                                      >  1:1076742/1‑150 (MQ=255)
                                                            gggTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGataata                                                                                      >  2:1022165/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                                                                                  <  1:1145087/135‑1 (MQ=255)
                                                                               cccACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGAt                                                                                  >  2:1145087/1‑135 (MQ=255)
                                                                               cccACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCa                                                                   >  1:975678/1‑150 (MQ=255)
                                                                               cccACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCa                                                                   >  2:660129/1‑150 (MQ=255)
                                                                                 cacaGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCAcc                                                                 >  1:1281691/1‑150 (MQ=255)
                                                                                 cacaGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCAcc                                                                 <  2:744721/150‑1 (MQ=255)
                                                                                    aGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCa                                                              <  1:365189/150‑1 (MQ=255)
                                                                                     gCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCaa                                                             >  2:1032103/1‑150 (MQ=255)
                                                                                      cAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAag                                                            >  2:942339/1‑150 (MQ=255)
                                                                                           ggTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGaa                                                       >  2:468197/1‑150 (MQ=255)
                                                                                           ggTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGaa                                                       >  1:169478/1‑150 (MQ=255)
                                                                                             ttGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAata                                                     <  2:1494155/150‑1 (MQ=255)
                                                                                               gggCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATg                                                   <  1:1514171/150‑1 (MQ=255)
                                                                                               gggCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATg                                                   <  1:1332312/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                  <  1:111610/150‑1 (MQ=255)
                                                                                                                                cGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCgttcagttc                  >  2:316803/1‑150 (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)
                                                                                                                                              ggTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACTGCTCAAgg    >  1:1190916/1‑150 (MQ=255)
                                                                                                                                               gTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACTGCTCAAGGc   <  1:1207578/150‑1 (MQ=255)
                                                                                                                                                tCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACTGCTCAAGGCa  <  2:1076742/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TGCACGCCCAACCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACTGCTCAAGGCA  >  NC_000913/3906125‑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: