Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,906,264 A→C 6.2% C102G (TGC→GGC) ‡ bglG ← transcriptional antiterminator BglG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,906,2640AC 6.2% 243.9 / 7.7 81C102G (TGC→GGC) ‡bglGtranscriptional antiterminator BglG
Reads supporting (aligned to +/- strand):  ref base A (42/34);  new base C (2/3);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

TAACTGCACGCCCAACCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACTGCTC  >  NC_000913/3906121‑3906412
                                                                                                                                               |                                                                                                                                                    
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: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:1099866/1‑150 (MQ=255)
                                                           ctTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAgaga                                                                                     >  2:1569452/1‑150 (MQ=255)
                                                           ctTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAgaga                                                                                     <  1:852281/150‑1 (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                                                                            >  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                                                 >  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: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: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)
                                                                                                                                              cAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACTGCTc  >  1:305482/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TAACTGCACGCCCAACCGTTTATCAATAATGGTTAACGCTTCTTCCCCTAGCTGGAACTCTTTCGGGTAAAGCCGCTGGATATCCCACAGCAACGGGTTGGGCAGCAACACGTTTTGCTGAAAGCGTTTAATCGCAAACTGGCAATGGTCAGTTAGCGAGATATAAATACTGTCCTGTAATTTTCCCAAGCGCTCCTGCGCTAAAGAGATAATACGATCACAGGTTGCCATCACCTCAAGAGGAATATGACTTAAGAGTTCGCTTAATCGCCCGTTCAGTTCATGACTGCTC  >  NC_000913/3906121‑3906412

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: