Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,906,256 A→T 6.0% F104L (TTT→TTA bglG ← transcriptional antiterminator BglG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,906,2560AT 6.0% 245.9 / 3.5 79F104L (TTT→TTAbglGtranscriptional antiterminator BglG
Reads supporting (aligned to +/- strand):  ref base A (40/34);  new base T (2/3);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCTTTCGGTAACTGCACGCCCAACCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATG  >  NC_000913/3906111‑3906405
                                                                                                                                                 |                                                                                                                                                     
cATCTTTCGGTAACTGCACGCCCAACCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACt                                                                                                                                                   <  1:887857/150‑1 (MQ=255)
   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                                                                                               <  1:247487/150‑1 (MQ=255)
                                                    ttcCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAgcg                                                                                               <  2:820596/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                                                                               <  2:1294563/150‑1 (MQ=255)
                                                                    tctTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAagag                                                                               <  1:1575273/150‑1 (MQ=255)
                                                                     ctTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAgaga                                                                              >  2:1569452/1‑150 (MQ=255)
                                                                     ctTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAgaga                                                                              <  1:852281/150‑1 (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                                                      >  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                                          >  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     >  2:8697/1‑150 (MQ=255)
                                                                                                                                              cGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCgttcagttc     <  1:111610/150‑1 (MQ=255)
                                                                                                                                                cAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCAt   <  1:1339815/150‑1 (MQ=255)
                                                                                                                                                 aaaCTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATg  <  2:1281691/150‑1 (MQ=255)
                                                                                                                                                 aaaCTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATg  <  1:660129/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CATCTTTCGGTAACTGCACGCCCAACCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATG  >  NC_000913/3906111‑3906405

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: