Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,906,262 G→A 6.1% C102C (TGC→TGT) ‡ bglG ← transcriptional antiterminator BglG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,906,2620GA 6.1% 250.6 / 2.3 81C102C (TGC→TGT) ‡bglGtranscriptional antiterminator BglG
Reads supporting (aligned to +/- strand):  ref base G (40/35);  new base A (2/3);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTCGGTAACTGCACGCCCAACCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACT  >  NC_000913/3906114‑3906408
                                                                                                                                                    |                                                                                                                                                  
cTTTCGGTAACTGCACGCCCAACCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGc                                                                                                                                                   <  1:64129/150‑1 (MQ=255)
       tAACTGCAAGCCCAACCGTTTATCAATAATGGATAACGATTATTCCCATAGCTGGAACTCTTTAGGGTAAAGCCCCTCGTAATCCCACAGAAACTGGTTCGGCAGCAACACGTTTTGCTTAAAGCTTTTAATCGAAAACTGCCAATgcta                                                                                                                                            >  2:448050/1‑147 (MQ=255)
           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: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                                                                                 >  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                                                                            >  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                                                         >  2:660129/1‑150 (MQ=255)
                                                                                          cccACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCa                                                         >  1:975678/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                                             >  2:468197/1‑150 (MQ=255)
                                                                                                      ggTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGaa                                             >  1:169478/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:531220/150‑1 (MQ=255)
                                                                                                                                 aGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGcc                  <  2:1358854/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:8697/1‑150 (MQ=255)
                                                                                                                                           cGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCgttcagttc        >  2:316803/1‑150 (MQ=255)
                                                                                                                                           cGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCgttcagttc        <  1:111610/150‑1 (MQ=255)
                                                                                                                                             cAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCAt      <  1:1339815/150‑1 (MQ=255)
                                                                                                                                              aaaCTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATg     <  1:660129/150‑1 (MQ=255)
                                                                                                                                              aaaCTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATg     <  2:1281691/150‑1 (MQ=255)
                                                                                                                                                 cTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACt  >  2:352560/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTTTCGGTAACTGCACGCCCAACCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACT  >  NC_000913/3906114‑3906408

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: