Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 821,875 C→T 5.5% Q208Q (CAG→CAA)  ybhN ← UPF0104 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913821,8750CT 5.5% 227.9 / 2.1 72Q208Q (CAG→CAA) ybhNUPF0104 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (47/21);  new base T (2/2);  total (49/23)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

CCGATCCCCGCCGGAATATGCACGATGACGCCAGCAATACTACTAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAATCATCAGTAAGCCAATGCCGAGAATACGCAGCGTAGTTTGATCGA  >  NC_000913/821729‑822023
                                                                                                                                                  |                                                                                                                                                    
ccGATCCCCGCCGGAATATGCACGATGACGCCAGCAATACTACTAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACAATCCAGTTAACGCTGGAAATCAGCATCAGa                                                                                                                                                   >  1:365522/1‑150 (MQ=255)
  gATCCCCGCCGGAATATGCACGATGACGCCAGCAATACTACTAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCGCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGc                                                                                                                                                 >  2:174652/1‑150 (MQ=255)
  gATCCCCGCCGGAATATGCACGATGACGCCAGCAATACTACTAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGc                                                                                                                                                 >  1:683245/1‑150 (MQ=255)
     ccccGCCGGAATATGCACGATGACGCCAGCAATACTACTAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAg                                                                                                                                              <  2:251591/150‑1 (MQ=255)
        cgccgGAATATGCACGATGACGCCAGCAATACTACTAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGc                                                                                                                                           >  1:650481/1‑150 (MQ=255)
         gccgGAATATGCACGATGACGCCAGCAATACTACTAACCAGTACCACGCCGAGTACAAAGAAATAGTTCACGCTTTGACCAAGAAACAGCCAGATAAACGCCCCCATTACCATCAAGATAACGCTCGAAATCAGCTTCTGAGCGacgccg                                                                                                                                          >  2:435653/1‑145 (MQ=255)
              aaTATGCACGATGACGCCAGCAATACTACTAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAAttt                                                                                                                                     >  2:284587/1‑150 (MQ=255)
                 atGCACGATGACGCCAGCAATACTACTAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATAGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCa                                                                                                                                  >  2:609469/1‑150 (MQ=255)
                 atGCACGATGACGCCAGCAATACTACTAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATAGCCCCCATTACCATCCAGTTAACGCTGGAAAAGAGCATCTGCGCGAGGGCGAATTTCCa                                                                                                                                  >  2:821975/1‑150 (MQ=255)
                   gCACGATGACGCCAGCATTACTACTAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGCCCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATg                                                                                                                                >  1:957886/1‑150 (MQ=255)
                   gCACGATGACGCCAGCAATACTACTAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATg                                                                                                                                <  1:1244493/150‑1 (MQ=255)
                    cACGATGACGCCAGCAATACTACTAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGa                                                                                                                               >  1:347507/1‑150 (MQ=255)
                      cGATGACGCCAGCAATACTACTAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTCAGCGAGGGCGAATTTCCATTAAg                                                                                                                             >  1:615705/1‑150 (MQ=255)
                        aTGACGCCAGCAATACTACTAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGATGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGgc                                                                                                                           >  1:940582/1‑150 (MQ=255)
                                  cAATACTACTAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATTTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCt                                                                                                                 >  1:1101165/1‑150 (MQ=255)
                                  cAATACTACTAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATTTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCt                                                                                                                 >  1:1330349/1‑150 (MQ=255)
                                    aaactactAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCttt                                                                                                               <  1:844100/148‑1 (MQ=255)
                                      actactAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTg                                                                                                             >  2:32281/1‑150 (MQ=255)
                                       ctactaACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTATCAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGCGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGc                                                                                                            >  2:735210/1‑150 (MQ=255)
                                        tactaACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGcc                                                                                                           <  2:477653/150‑1 (MQ=255)
                                        tactaACCAGTAACACGCCCAGTACAAAGAAATAGATCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGcc                                                                                                           <  2:388933/150‑1 (MQ=255)
                                         actaCCCAGTAACACGCCCAGTAGAAAGAAATTGTTGCCGCTTTGCCCAAGTAACAGCCAGATAATCGCCCCCATTTCCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGA‑GGCGAATTTCCATGAAGGCAGCACCCGCTTTTGcccg                                                                                                         >  2:992269/1‑149 (MQ=255)
                                         actaACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGccc                                                                                                          >  1:11423/1‑150 (MQ=255)
                                           taACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCtt                                                                                                        >  1:650616/1‑150 (MQ=255)
                                             aCCAGTAACACGCCCCGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCACATAATCGCCCCCATTACCATCCAGTTAACGCAGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCGTTTGCCCTTTg                                                                                                      <  1:80180/150‑1 (MQ=255)
                                                   aaCACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTc                                                                                                >  2:314292/1‑150 (MQ=255)
                                                     cacGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCat                                                                                              >  1:445880/1‑150 (MQ=255)
                                                     cacGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCATTTGATGGTCat                                                                                              >  2:86934/1‑150 (MQ=255)
                                                         cccAGTACAAAGAAATAGTTCACGCTTTGACCATGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCCCCAGCTTTTGCCCGTTGATGGTCATATgg                                                                                          >  1:44827/1‑150 (MQ=255)
                                                         cccAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGAATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATgg                                                                                          >  2:1194388/1‑150 (MQ=255)
                                                           cAGTACAAAGAAATAGTTCCCGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTGCATGAGGGCAGCACCAGCTTTTGCCCTGTGATGGTCATATGGCg                                                                                        >  2:142093/1‑150 (MQ=255)
                                                               aCAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGc                                                                                    <  1:862859/150‑1 (MQ=255)
                                                                       aTAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGaa                                                                            <  2:748971/150‑1 (MQ=255)
                                                                            tCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAccgc                                                                       >  1:797591/1‑150 (MQ=255)
                                                                             aaCGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCa                                                                      <  1:123923/149‑1 (MQ=255)
                                                                                gCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGAAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCaaaa                                                                   >  2:854656/1‑150 (MQ=255)
                                                                                 cTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCCCCAGCTTTTGCCCTTTGATGGTCTTATGGCGGTGCTTCGCGAAAGCGCAAAAc                                                                  >  2:379824/1‑150 (MQ=255)
                                                                                 cTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAAc                                                                  <  2:650616/150‑1 (MQ=255)
                                                                                 cTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAAc                                                                  >  1:430913/1‑150 (MQ=255)
                                                                                  tttGACCAAGTAACAGCCAGATACTCGCCCCCATTACCATCCCGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCATTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAAcc                                                                 <  1:834126/150‑1 (MQ=255)
                                                                                     gACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACcaca                                                              <  2:579454/150‑1 (MQ=255)
                                                                                      aCCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAATTCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACaa                                                             >  2:511224/1‑150 (MQ=255)
                                                                                      aCCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACaa                                                             >  1:1296256/1‑150 (MQ=255)
                                                                                        cAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAAt                                                           >  1:1273740/1‑150 (MQ=255)
                                                                                         aaGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGGCATATGGCGGTGCTGCGCGAAAGCGCAAGACCACCAATa                                                          >  2:1015319/1‑150 (MQ=255)
                                                                                         aaGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACAAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATa                                                          >  2:1021097/1‑150 (MQ=255)
                                                                                              aCAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTCAGCGAGGGCGAATATCCATGAAGGCATCACCAGCTTTTGCCCTGTGAAGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACcg                                                     >  1:131782/1‑150 (MQ=255)
                                                                                               cAGCCAGATACTCGCCCCGATTACCATCCAGTTAAGGCTGGAAATCAGCATCTGAGCGAGTGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACcgc                                                    <  2:992209/150‑1 (MQ=255)
                                                                                                 gCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATTTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGa                                                  <  2:1330349/150‑1 (MQ=255)
                                                                                                 gCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATTTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGa                                                  <  2:1101165/150‑1 (MQ=255)
                                                                                                 gCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAAGTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCCTCGCGAAAGCGCAAAACCACAAATAAACCGCGa                                                  >  1:646553/1‑150 (MQ=255)
                                                                                                  ccAGATAATCGCCCCCATTACCATCCAGTTAACGCAGGAAATCAGTATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGAt                                                 <  1:58574/150‑1 (MQ=255)
                                                                                                   cAGAAAATCGCCCCCTATACCATCCCGTTAACGCAGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCCCCAGCTTTTGCCCTTTGATGGTCAAATGGCGGGGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATa                                                <  2:271014/150‑1 (MQ=255)
                                                                                                         aTCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAatcatc                                          <  1:839775/150‑1 (MQ=255)
                                                                                                               cccATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAATCATCAGTAAg                                    >  1:431319/1‑150 (MQ=255)
                                                                                                                 cATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAATCATCAGTAAGcc                                  <  2:80510/150‑1 (MQ=255)
                                                                                                                 cATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAATCATCAGTAAGcc                                  <  2:997250/150‑1 (MQ=255)
                                                                                                                    tACCATCCAGTTAACGCTGGTAATCAGCATCTGAGCGAGGGCGAATTTCCATGAATGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAATCATCAGTAAGCCAAt                               <  1:1218902/150‑1 (MQ=255)
                                                                                                                      ccATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAATCATCAGTAAGCCAATGc                             >  2:757192/1‑150 (MQ=255)
                                                                                                                      ccATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAATCATCAGTAAGCCAATGc                             >  1:309348/1‑150 (MQ=255)
                                                                                                                        aTCCAGTTAACGATGGAAATCAGCATCTGAGCGAGGGAGAATTTCAATGAAGGCAGAACCAGCTTTTGACCTTTGATGGTCATAAGGCGGTGCTTCGCGAAAGCGCAAACCCACAAAAAAACCGCGATAATCAACAGTAAGCCAATGCCg                           >  2:356796/1‑150 (MQ=255)
                                                                                                                            aGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCCGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAATCATCAGTAAGCCAATGCCGAGaa                       >  2:228184/1‑150 (MQ=255)
                                                                                                                                   gCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAATCATCAGTAAGCCAATGCCGAGAATACGCAg                >  1:1160951/1‑150 (MQ=255)
                                                                                                                                    cTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAATCATCAGTAAGCCAATGCCGAGAATACGCAGc               >  1:411928/1‑150 (MQ=255)
                                                                                                                                     tGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAATCATCAGTAAGCCAATGCCGAGAATACGCAGCg              >  2:1015897/1‑150 (MQ=255)
                                                                                                                                      ggAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAATCATCAGTAAGCCAATGCCGAGAATACGCAGCGt             >  1:154893/1‑150 (MQ=255)
                                                                                                                                          aTCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAATCATCAGTAAGCCAATGCCGAGAATACGCAGCGTAGtt         >  1:1113985/1‑150 (MQ=255)
                                                                                                                                            cAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTTCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAATCATCAGTAAGCCAATGCCGAGAATACGCAGCGTAGTTTg       <  2:683245/150‑1 (MQ=255)
                                                                                                                                            cAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAATCATCAGTAAGCCAATGCCGAGAATACGCAGCGTAGTTTg       >  2:656894/1‑150 (MQ=255)
                                                                                                                                            cAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAATCATCAGTAAGCCAATGCCGAGAATACGCAGCGTAGTTTg       >  2:673016/1‑150 (MQ=255)
                                                                                                                                               cATCTGAGCGAGGGCGAATTACCAAGAAGCCAGCACCAGCTGTTGCCCTTTTAAGGTCATATGGCGGTGCTTCGCGAACGCGCAATACCACAAATAACCCGCGATCATCATCAGTAACCCAATGACGAGAATACGCAGCGTAGTCTGATc    <  2:651273/150‑1 (MQ=255)
                                                                                                                                                 tCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAATCATCAGTAAGCCAATGCCGAGAATACGCAGCGTAGTTTGATCGa  >  2:63369/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCGATCCCCGCCGGAATATGCACGATGACGCCAGCAATACTACTAACCAGTAACACGCCCAGTACAAAGAAATAGTTCACGCTTTGACCAAGTAACAGCCAGATAATCGCCCCCATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAACCACAAATAAACCGCGATAATCATCAGTAAGCCAATGCCGAGAATACGCAGCGTAGTTTGATCGA  >  NC_000913/821729‑822023

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