Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,282,906 C→A 7.2% intergenic (+105/+34) rplY → / ← yejK 50S ribosomal subunit protein L25/spermidine nucleoid‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,282,9060CA 7.2% 367.3 / 13.8 129intergenic (+105/+34)rplY/yejK50S ribosomal subunit protein L25/spermidine nucleoid‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (59/60);  new base A (4/6);  total (63/66)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

Only 100 of 132 total aligned reads displayed.

CCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCG  >  NC_000913/2282757‑2283055
                                                                                                                                                     |                                                                                                                                                     
ccGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTc                                                                                                                                                       >  2:209844/1‑150 (MQ=255)
ccGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTc                                                                                                                                                       <  2:31742/150‑1 (MQ=255)
 cGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCa                                                                                                                                                      <  1:1246344/150‑1 (MQ=255)
  gTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAc                                                                                                                                                     >  2:361446/1‑150 (MQ=255)
   tACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGACTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcg                                                                                                                                                    >  2:691711/1‑150 (MQ=255)
   tACAAACCGAAGCTGCAGCACATCGACTGCGATCGCGCTAAAATGCTGAATAATTTGTAAAAAACCCCGCTCCGGCGGGTTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcg                                                                                                                                                    <  1:1005686/150‑1 (MQ=255)
     cacaCCCAAGCTGCAGCACCTCGTCTTCGTTCGCGCTTAATGGCTGAATAAGTAGTAAAAAACCCCGCGCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcgcc                                                                                                                                                  <  2:321313/147‑1 (MQ=255)
     cAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcgcc                                                                                                                                                  <  1:1268250/150‑1 (MQ=255)
      aaaCCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGAAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACTGTATTGCTATTTTTTGCAGGCAAGAAAAGGCGTTCAcgccg                                                                                                                                                 >  1:186188/1‑150 (MQ=255)
       aaCCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCCGGGTTTTTTGTATCTGCAGATTATCCCTGATAACTGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTAAcgccgc                                                                                                                                                >  2:1728877/1‑150 (MQ=255)
          cGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATc                                                                                                                                             >  2:558925/1‑150 (MQ=255)
            aaGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTg                                                                                                                                           <  1:1432331/150‑1 (MQ=255)
              gCTGCAGCACATCGACTTCGTTCGCGATTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTGGc                                                                                                                                         <  2:1193763/150‑1 (MQ=255)
                tGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTGGcct                                                                                                                                       <  1:945967/150‑3 (MQ=255)
                  cagcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTc                                                                                                                                     <  2:1195097/150‑1 (MQ=255)
                  cagcaCATCGACGTCGTTCGCTCTTAATTGCTGAATAAGTTGTAAATAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTGGcctgc                                                                                                                                     <  1:345835/150‑5 (MQ=255)
                   agcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCTGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCa                                                                                                                                    <  1:1082572/150‑1 (MQ=255)
                   agcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCa                                                                                                                                    <  2:1228572/150‑1 (MQ=255)
                   agcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCCATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTGGcctgca                                                                                                                                    <  1:431386/150‑6 (MQ=255)
                    gcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCATTCACGCCGCATCTGGCATTCAg                                                                                                                                   >  1:1562713/1‑150 (MQ=255)
                     cacaTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGa                                                                                                                                  <  1:566680/150‑1 (MQ=255)
                       caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATTCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAac                                                                                                                                >  1:48711/1‑150 (MQ=255)
                       caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAac                                                                                                                                >  2:345946/1‑150 (MQ=255)
                       caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAac                                                                                                                                >  1:1329811/1‑150 (MQ=255)
                       caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAac                                                                                                                                <  1:654034/150‑1 (MQ=255)
                       caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAac                                                                                                                                <  2:1772853/150‑1 (MQ=255)
                         tCGACTTCGGTGGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTTCGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAacac                                                                                                                              <  2:1760336/150‑1 (MQ=255)
                               tcgttcgCGCTTAATTGCTGATTAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGAATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAAcc                                                                                                                        <  2:233419/150‑1 (MQ=255)
                               tcgttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGTGTTTTTTGTATCTGCAGATTATGCCTGATTACGGAATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAAcc                                                                                                                        <  2:369178/150‑1 (MQ=255)
                               tcgttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAAcc                                                                                                                        >  2:431985/1‑150 (MQ=255)
                                   tcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGAtt                                                                                                                    <  2:1522421/150‑1 (MQ=255)
                                   tcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGAtt                                                                                                                    >  1:62720/1‑150 (MQ=255)
                                   tcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGAtt                                                                                                                    <  2:1052962/150‑1 (MQ=255)
                                   tcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCAGTCAGAACACAAAACCGAtt                                                                                                                    >  1:1229288/1‑150 (MQ=255)
                                         ttAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTgc                                                                                                              >  1:452530/1‑150 (MQ=255)
                                                cTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGa                                                                                                       <  2:1164262/150‑1 (MQ=255)
                                                 tGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCGTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGAc                                                                                                      >  1:1312737/1‑150 (MQ=255)
                                                  gAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGCCATTCAGAACACAAAACCGATTAATTGCCGCCAGACg                                                                                                     >  2:188955/1‑150 (MQ=255)
                                                    aTAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTc                                                                                                   <  2:1157703/150‑1 (MQ=255)
                                                        gTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGgc                                                                                               >  2:1662998/1‑150 (MQ=255)
                                                          tGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGgcgc                                                                                             >  1:744565/1‑150 (MQ=255)
                                                           gTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGAGTATGCCTGATTACGGTATTGCGAGTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCCTTCAGCACACAACACCGAGTAATTGCCGCCAGACGTCCGGCGCt                                                                                            >  1:313106/1‑150 (MQ=255)
                                                             aaaaaaCCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGc                                                                                          >  2:1960991/1‑150 (MQ=255)
                                                              aaaaaCCCCGCGCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTGCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCa                                                                                         <  1:176441/150‑1 (MQ=255)
                                                                  aCCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTg                                                                                     <  2:692298/150‑1 (MQ=255)
                                                                   ccccGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTgg                                                                                    >  1:215165/1‑150 (MQ=255)
                                                                    cccGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGt                                                                                   <  2:403246/150‑1 (MQ=255)
                                                                        cTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTcgcg                                                                               <  1:213631/150‑1 (MQ=255)
                                                                        cTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTcgcg                                                                               <  2:751537/150‑1 (MQ=255)
                                                                                gggTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTcg                                                                       <  1:715750/150‑1 (MQ=255)
                                                                                  gTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTcggc                                                                     <  2:1524648/150‑1 (MQ=255)
                                                                                  gTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTcggc                                                                     >  2:904909/1‑150 (MQ=255)
                                                                                    tttttGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTcggcgg                                                                   <  1:1855337/150‑1 (MQ=255)
                                                                                     ttttGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCggcggc                                                                  <  2:1713031/150‑1 (MQ=255)
                                                                                       ttGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGgcggcgg                                                                >  2:564225/1‑149 (MQ=255)
                                                                                         gTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGc                                                              >  2:566593/1‑150 (MQ=255)
                                                                                            tCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCTGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCtt                                                           >  2:1618155/1‑150 (MQ=255)
                                                                                             cTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCTGCCAGACGTCCGGCGCTGCAATTGGTCGCGAAAATTCGGCGGCGTGCCttt                                                          >  2:1503682/1‑150 (MQ=255)
                                                                                               gCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTaa                                                        <  1:209844/150‑1 (MQ=255)
                                                                                                cAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAAt                                                       <  1:398924/150‑1 (MQ=255)
                                                                                                 aGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATg                                                      <  1:1776759/150‑1 (MQ=255)
                                                                                                  gATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATgg                                                     >  2:1212459/1‑150 (MQ=255)
                                                                                                  gATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATgg                                                     <  1:1927166/150‑1 (MQ=255)
                                                                                                    ttATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTc                                                   <  1:55716/150‑1 (MQ=255)
                                                                                                      aTGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCaa                                                 <  1:1024680/150‑1 (MQ=255)
                                                                                                       tGCCTGATTACGGAATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAACCCGATAAATTGCCGCCAGACGGCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAg                                                <  2:261247/150‑1 (MQ=255)
                                                                                                        gCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGAGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAgg                                               >  2:659750/1‑150 (MQ=255)
                                                                                                           tGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTAt                                            >  2:1792066/1‑150 (MQ=255)
                                                                                                            gATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTCGCATTCAGAACACAACACCGATTAATTGCCGCCAGACGACCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGACAAGGTATc                                           >  2:1342687/1‑150 (MQ=255)
                                                                                                               tACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCTCTTAATTGCCGCCAGACGTCCGGCGATGCAATTGGTCGCGCAAAGTCGGCGGCGTGCCTTTAATGTGCAAGGTATCggg                                        >  1:321147/1‑149 (MQ=255)
                                                                                                                  ggTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCTGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGc                                     <  1:361446/150‑1 (MQ=255)
                                                                                                                    tATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTg                                   <  1:1548007/150‑1 (MQ=255)
                                                                                                                      ttGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGa                                 <  1:345946/150‑1 (MQ=255)
                                                                                                                        gCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATc                               >  1:1735931/1‑150 (MQ=255)
                                                                                                                         cTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGCCGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATcc                              >  1:1740306/1‑150 (MQ=255)
                                                                                                                         cTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATcc                              >  1:930290/1‑150 (MQ=255)
                                                                                                                            ttttttGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCag                           <  2:1103415/150‑1 (MQ=255)
                                                                                                                            ttttttGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCGGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCag                           <  2:1298861/150‑1 (MQ=255)
                                                                                                                            tttcttGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTAAGAACACACAACCGATTCATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATTGTCAAGGTATCGGTTGCTGGATCCCag                           <  1:108659/150‑1 (MQ=255)
                                                                                                                             tttttGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCaga                          >  2:899840/1‑150 (MQ=255)
                                                                                                                             tttttGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCaga                          <  2:48711/150‑1 (MQ=255)
                                                                                                                               tttGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCagaag                        <  1:63725/150‑1 (MQ=255)
                                                                                                                                ttGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCagaaga                       <  2:1329811/150‑1 (MQ=255)
                                                                                                                                ttGCAGGCCAGATAAGGCGTTCACGCAGCATCTGGAATTCAGAACACAAACCCGATTAATTGCAGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGTATCCAagaaga                       >  1:1591152/1‑150 (MQ=255)
                                                                                                                                 tGCAGGCCAGATAAGGCGTTAACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGCCGTCCGGCGATGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGAt                      <  2:186188/150‑1 (MQ=255)
                                                                                                                                  gCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATg                     >  2:1819995/1‑150 (MQ=255)
                                                                                                                                  gCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATg                     <  2:1666432/150‑1 (MQ=255)
                                                                                                                                   cAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGc                    >  2:499632/1‑150 (MQ=255)
                                                                                                                                   cAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGc                    >  2:74819/1‑150 (MQ=255)
                                                                                                                                    atgCCAGATGCGGCGTGAACGCCTTATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGAGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCTGAAGATGCg                   >  1:247635/3‑150 (MQ=255)
                                                                                                                                    aGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCg                   <  2:746029/150‑1 (MQ=255)
                                                                                                                                     ggCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGt                  >  1:1804455/1‑150 (MQ=255)
                                                                                                                                     ggCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGt                  >  2:847968/1‑150 (MQ=255)
                                                                                                                                     ggCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGt                  >  1:421803/1‑150 (MQ=255)
                                                                                                                                        cAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGGTCCCAGAAGATGCGTTCa               >  2:1073029/1‑150 (MQ=255)
                                                                                                                                          gATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAAGTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCAcc             <  1:413062/150‑1 (MQ=255)
                                                                                                                                             gcggCGTGAACGCCTTATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGACGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTaa          >  2:1882350/3‑150 (MQ=255)
                                                                                                                                                  gTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGGTGCTGGATCCCAGAAGAGGCGTTCACCTAACAgca     >  2:1331215/1‑150 (MQ=255)
                                                                                                                                                    tCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTAGGCGGCGTGCCTTTAATGGTCAATGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAgcatc   <  1:274532/150‑1 (MQ=255)
                                                                                                                                                     cACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAgcatcg  >  1:1672439/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCG  >  NC_000913/2282757‑2283055

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