Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,643,116 C→A 5.4% Q42H (CAG→CAT) ‡ ydfT ← Qin prophage; putative antitermination protein Q

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,643,1160CA 5.4% 245.9 / 2.9 74Q42H (CAG→CAT) ‡ydfTQin prophage; putative antitermination protein Q
Reads supporting (aligned to +/- strand):  ref base C (32/37);  new base A (2/2);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

ATCTGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTA  >  NC_000913/1642970‑1643257
                                                                                                                                                  |                                                                                                                                             
atCTGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGc                                                                                                                                            >  1:1180060/1‑150 (MQ=255)
atCTGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGc                                                                                                                                            <  2:211617/150‑1 (MQ=255)
   tGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCAcc                                                                                                                                         >  2:243205/1‑150 (MQ=255)
       aCCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCAGTCGAAAGCACCactt                                                                                                                                     >  2:803657/1‑147 (MQ=255)
        ccGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCAGTCGAAAGCACCacttg                                                                                                                                    <  1:126605/150‑5 (MQ=255)
        ccGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCAGTCGAAAGCACCacttg                                                                                                                                    <  1:803657/150‑5 (MQ=255)
          ggATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACCTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCt                                                                                                                                  >  1:76017/1‑150 (MQ=255)
              ttcttcACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTgcgc                                                                                                                              <  1:1423987/150‑1 (MQ=255)
                cttcACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTg                                                                                                                            >  2:522466/1‑150 (MQ=255)
                 ttcACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGc                                                                                                                           <  2:86607/150‑1 (MQ=255)
                  tcACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGcc                                                                                                                          <  1:151555/150‑1 (MQ=255)
                  tcACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGcc                                                                                                                          <  2:918099/150‑1 (MQ=255)
                   cACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCa                                                                                                                         <  2:366659/150‑1 (MQ=255)
                   cACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCa                                                                                                                         <  2:885654/150‑1 (MQ=255)
                   cACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCa                                                                                                                         >  2:653338/1‑150 (MQ=255)
                   cACTCAATTTGTTCAATGTCGGGTTATCCAGAGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCa                                                                                                                         >  1:805804/1‑150 (MQ=255)
                   cACTCAATTTGTTCAATGACGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCa                                                                                                                         >  2:690942/1‑150 (MQ=255)
                      tCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCCCCATCCCCTGCGCTGCCATTa                                                                                                                      >  1:714335/1‑150 (MQ=255)
                      tCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTa                                                                                                                      <  1:243205/150‑1 (MQ=255)
                       cAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTAc                                                                                                                     >  2:261085/1‑150 (MQ=255)
                            tGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGTGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCg                                                                                                                <  2:367210/150‑1 (MQ=255)
                             gTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGt                                                                                                               >  1:1116143/1‑150 (MQ=255)
                              ttCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTa                                                                                                              >  2:1469366/1‑150 (MQ=255)
                                cAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAAt                                                                                                            <  1:188756/150‑1 (MQ=255)
                                   tGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGACCACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGt                                                                                                         >  1:1425619/1‑150 (MQ=255)
                                     tCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCATCTTTTGCCAGAAATAAGGCCAACCCAAGTGTTGCTTTCGACTGCACCATCCCCTTCGCTGCCATTACATACGTAATTGTTa                                                                                                       <  2:158838/150‑1 (MQ=255)
                                     tCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACCTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTa                                                                                                       <  2:76017/150‑1 (MQ=255)
                                      cGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTaa                                                                                                      <  1:1463776/150‑1 (MQ=255)
                                        ggTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTGCATCCGTAATTGTTAaac                                                                                                    >  2:1257080/1‑150 (MQ=255)
                                        ggTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAaac                                                                                                    <  1:768749/150‑1 (MQ=255)
                                          ttATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACCCCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAaacaa                                                                                                  >  1:1225745/1‑150 (MQ=255)
                                            aTCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAAcc                                                                                                >  2:677183/1‑150 (MQ=255)
                                            aTCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTGAAACAAcc                                                                                                >  2:1270921/1‑150 (MQ=255)
                                             tCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCg                                                                                               >  1:261338/1‑150 (MQ=255)
                                              ccAGTGCCATCGAGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAAACCAAGTGGTGCATTCGACTGAACAAACCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCgg                                                                                              >  2:331457/1‑150 (MQ=255)
                                                  tGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGcc                                                                                          <  2:1400179/150‑1 (MQ=255)
                                                  tGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGcc                                                                                          <  1:653338/150‑1 (MQ=255)
                                                    ccATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTg                                                                                        <  1:1400406/150‑1 (MQ=255)
                                                         gcgTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCt                                                                                   >  2:1407975/1‑150 (MQ=255)
                                                           gTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTgg                                                                                 >  2:1138547/1‑150 (MQ=255)
                                                              aTTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGc                                                                              >  2:653325/1‑150 (MQ=255)
                                                                   gCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCAt                                                                         >  2:795931/1‑150 (MQ=255)
                                                                    cAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATc                                                                        <  1:559737/150‑1 (MQ=255)
                                                                     aGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCg                                                                       <  1:261085/150‑1 (MQ=255)
                                                                     aGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGGGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCg                                                                       <  1:208355/150‑1 (MQ=255)
                                                                            cAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAAtt                                                                >  2:1070831/1‑150 (MQ=255)
                                                                              aTCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAAtttt                                                              >  1:863735/1‑150 (MQ=255)
                                                                                cGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTgg                                                            <  1:493737/150‑1 (MQ=255)
                                                                                 gCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCGCTGCTCTGCCATTACATCCGTAATTGTTAAACTACCGGTGCCTGTCGGTGGAGCGTCATCGCTCAATTTTGga                                                           <  2:69490/150‑1 (MQ=255)
                                                                                  cAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGATTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGATCAATTTTGgaa                                                          >  2:963587/1‑149 (MQ=255)
                                                                                    aaCTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGAt                                                        <  1:1138547/150‑1 (MQ=255)
                                                                                          gggTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTggg                                                  <  1:1410695/150‑1 (MQ=255)
                                                                                            gTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGa                                                <  1:622043/150‑1 (MQ=255)
                                                                                                tGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTaa                                            >  1:1371041/1‑150 (MQ=255)
                                                                                                  aaCACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAata                                          <  1:1261796/150‑1 (MQ=255)
                                                                                                  aaCACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGACATTACATCCGTAATTGTTATACAAAATGTTACTGTCGCTGGAGCGTCATCGATCAATTTTGTAGATTTTGGTTATTAata                                          <  2:651362/150‑1 (MQ=255)
                                                                                                      ccAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATAtttt                                      >  1:457789/1‑150 (MQ=255)
                                                                                                        aaCTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGt                                                                                       <  1:1355224/99‑1 (MQ=255)
                                                                                                        aaCTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGt                                                                                       >  2:1355224/1‑99 (MQ=255)
                                                                                                        aaCTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTgg                                    <  2:1116143/150‑1 (MQ=255)
                                                                                                        aaCTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCAGCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTgg                                    >  1:1314957/1‑150 (MQ=255)
                                                                                                             ttGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAgg                               <  2:878332/150‑1 (MQ=255)
                                                                                                               gCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCt                             >  2:436292/1‑150 (MQ=255)
                                                                                                                 cAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTc                           >  2:1132710/1‑150 (MQ=255)
                                                                                                                             cAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCt               <  1:522466/150‑1 (MQ=255)
                                                                                                                               aCCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCg             <  2:1213701/150‑1 (MQ=255)
                                                                                                                               aCCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCg             <  2:805804/150‑1 (MQ=255)
                                                                                                                               aCCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCg             <  2:1425619/150‑1 (MQ=255)
                                                                                                                                     atggtgCAGTCGAAAGCACCATCCCCTGCGCTGCCATTACATCCATAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTcc       >  2:650617/2‑150 (MQ=255)
                                                                                                                                      gggtgCTTTCGAGTGGACCATCCCCTGCGCTGCCATTACATCCGTATTTGTTAAACAACGGGTGCCTGTCGCTGGTGCGTGATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAATGCTTCAAGGTTCATGCTCGTTCTCCa      <  1:331457/149‑1 (MQ=255)
                                                                                                                                       ggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCAc     <  1:653325/150‑1 (MQ=255)
                                                                                                                                       ggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCAc     <  1:1469366/150‑1 (MQ=255)
                                                                                                                                        gtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACt    <  1:690942/150‑1 (MQ=255)
                                                                                                                                          gCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTa  >  2:603144/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
ATCTGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTA  >  NC_000913/1642970‑1643257

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