Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,554,731 C→T 50.9% M313I (ATG→ATA maeA ← malate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,554,7310CT50.9% ‑5.7 / 133.0 89M313I (ATG→ATAmaeAmalate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme
Reads supporting (aligned to +/- strand):  ref base C (20/23);  new base T (26/19);  total (47/42)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

GGCATCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAGTGTG  >  NC_000913/1554633‑1554830
                                                                                                  |                                                                                                   
ggCATCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCa                                                                                                    >  1:2337729/1‑100 (MQ=255)
  cATCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGCTTAtt                                                                                                  >  2:1340399/1‑100 (MQ=255)
     attGTCAGTCAGCAAGCCAAAGCGATCGACGATCAAGACTTTCGTCCGCGCCGCTTCCTCGCTTAATCATCCGAGCTGGGTCTGGGAGGTCATAATTTCg                                                                                               >  2:463486/2‑100 (MQ=255)
      ttGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCgg                                                                                              <  1:420933/100‑1 (MQ=255)
       tGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGc                                                                                             <  2:517578/100‑1 (MQ=255)
       tGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGc                                                                                             <  1:1203439/100‑1 (MQ=255)
        gtcagtcaGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCa                                                                                            >  2:1424064/1‑100 (MQ=255)
         tcagtcagCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCaa                                                                                           >  2:681471/1‑100 (MQ=255)
           agtcagCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATg                                                                                         >  2:765221/1‑100 (MQ=255)
            gtcagCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGc                                                                                        >  2:1096550/1‑100 (MQ=255)
            gtcagCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGc                                                                                        <  1:389021/100‑1 (MQ=255)
            gtcagCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGc                                                                                        >  1:1472012/1‑100 (MQ=255)
             tcagCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGATATGATCATTTCGGCAATGcc                                                                                       <  1:2247129/100‑1 (MQ=255)
             tcagCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGcc                                                                                       <  1:1984717/100‑1 (MQ=255)
               agCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGc                                                                                     >  1:95535/1‑100 (MQ=255)
               agCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGc                                                                                     >  2:356362/1‑100 (MQ=255)
                   aGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATcc                                                                                 <  1:305952/100‑1 (MQ=255)
                       aaaGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCt                                                                             <  2:424057/100‑1 (MQ=255)
                       aaaGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCt                                                                             <  1:689544/100‑1 (MQ=255)
                       aaaGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCt                                                                             >  2:987679/1‑100 (MQ=255)
                          gCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGaa                                                                          >  1:1642810/1‑100 (MQ=255)
                          gCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGaa                                                                          <  1:733113/100‑1 (MQ=255)
                           cGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAAc                                                                         <  1:233797/100‑1 (MQ=255)
                            gATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCTGCATCCCGCTGAAcc                                                                        <  2:2337707/100‑1 (MQ=255)
                            gATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAAcc                                                                        <  1:289565/100‑1 (MQ=255)
                              tCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTg                                                                      >  1:514737/1‑100 (MQ=255)
                              tCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTg                                                                      <  2:267434/100‑1 (MQ=255)
                               cGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTgc                                                                     >  2:1698502/1‑100 (MQ=255)
                                gACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTgcg                                                                    <  1:1722578/100‑1 (MQ=255)
                                gACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTgcg                                                                    <  2:2208188/100‑1 (MQ=255)
                                gACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTgcg                                                                    >  2:852281/1‑100 (MQ=255)
                                 aCCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTgcgc                                                                   >  1:1404896/1‑100 (MQ=255)
                                  ccATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGcc                                                                  <  1:1272702/100‑1 (MQ=255)
                                    aTAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCaa                                                                >  2:1760580/1‑100 (MQ=255)
                                      aaaGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAgg                                                              >  2:362410/1‑100 (MQ=255)
                                      aaaGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAgg                                                              <  2:110126/100‑1 (MQ=255)
                                         gccTTTCTGCCGCGCCTCTTCCTCGCTTATTCCTTCACGCTGGGTCTGGGAGATGATCATTTCGGCAATGACGCATCCCGCTGACCCTGCGCCAAGGa                                                             <  2:1146174/96‑1 (MQ=255)
                                         gACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAg                                                           <  2:564448/100‑1 (MQ=255)
                                           cTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGAc                                                         >  2:1120163/1‑100 (MQ=255)
                                             ttCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGa                                                       <  2:1093420/100‑1 (MQ=255)
                                               cTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATACCGCTGAACCTGCGCCAAGGAAGACgct                                                      >  2:717027/1‑97 (MQ=255)
                                                tGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGAttt                                                    >  1:742466/1‑100 (MQ=255)
                                                  ccgcgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGAttttt                                                  <  1:356364/100‑1 (MQ=255)
                                                   ngcgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGAtttttt                                                 >  2:21671/2‑100 (MQ=255)
                                                   cgcgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGAtttttt                                                 >  2:1220865/1‑100 (MQ=255)
                                                    gcgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCTATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGAttttttt                                                >  2:2344330/1‑100 (MQ=255)
                                                     ttccgactCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTtc                                               <  1:2129941/93‑1 (MQ=255)
                                                     cgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTtc                                               <  2:1048009/100‑1 (MQ=255)
                                                     cgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTtc                                               <  2:2136627/100‑1 (MQ=255)
                                                     cgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTtc                                               <  1:765224/100‑1 (MQ=255)
                                                     cgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTtc                                               <  1:286426/100‑1 (MQ=255)
                                                      gccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTtct                                              <  1:1424072/100‑1 (MQ=255)
                                                       ccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTtctc                                             <  2:138908/100‑1 (MQ=255)
                                                       ccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTtctc                                             <  1:1220870/100‑1 (MQ=255)
                                                        cgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCg                                            >  2:1951937/1‑100 (MQ=255)
                                                        cgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCg                                            >  2:1952917/1‑100 (MQ=255)
                                                         gcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGc                                           >  2:2058075/1‑100 (MQ=255)
                                                          cTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCt                                          >  1:826456/1‑100 (MQ=255)
                                                            tCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTa                                        >  2:1226615/1‑100 (MQ=255)
                                                             ccTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTtcagatcgga                                       >  2:2129923/1‑92 (MQ=255)
                                                              cTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAAc                                      >  1:1478078/1‑100 (MQ=255)
                                                              cTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAAc                                      >  1:54951/1‑100 (MQ=255)
                                                               tCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACt                                     >  1:1704140/1‑100 (MQ=255)
                                                               tCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACt                                     <  1:501055/100‑1 (MQ=255)
                                                                 gCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGa                                   >  1:1184613/1‑100 (MQ=255)
                                                                 gCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGa                                   <  2:1235249/100‑1 (MQ=255)
                                                                     aaTCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGaccac                               >  2:503848/1‑100 (MQ=255)
                                                                     aaTCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGaccac                               >  2:712195/1‑100 (MQ=255)
                                                                      aTCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGaccacc                              >  1:2273094/1‑100 (MQ=255)
                                                                         cTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGc                           >  1:1035541/1‑100 (MQ=255)
                                                                         cTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGc                           >  2:776520/1‑100 (MQ=255)
                                                                         cTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGc                           <  1:1534794/100‑1 (MQ=255)
                                                                           tCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCcg                         <  2:2001837/100‑1 (MQ=255)
                                                                            cgcgCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCcgc                        >  1:1550580/1‑100 (MQ=255)
                                                                            cgcgCTGGGTCTGGGAGATGATTATTTCGGCAATGACGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCcgc                        <  2:104177/100‑1 (MQ=255)
                                                                             gcgcTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCcgcg                       <  1:1016366/100‑1 (MQ=255)
                                                                                 tGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGc                   <  1:636014/100‑1 (MQ=255)
                                                                                  gggTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCt                  <  1:1465008/100‑1 (MQ=255)
                                                                                  gggTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCt                  <  1:2044773/100‑1 (MQ=255)
                                                                                    gTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTtg                >  2:2343926/1‑100 (MQ=255)
                                                                                     tCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTtgc               >  1:63321/1‑100 (MQ=255)
                                                                                      cTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTtgct              >  2:497011/1‑100 (MQ=255)
                                                                                       tGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTtgctg             >  2:1712659/1‑100 (MQ=255)
                                                                                          gagaTGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGa          <  2:1642799/100‑1 (MQ=255)
                                                                                           agaTGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGAt         <  1:630283/100‑1 (MQ=255)
                                                                                              tgatTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAg      >  1:282988/1‑100 (MQ=255)
                                                                                               gatTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCTCTTATCTGACCACCTGCCGCGCGGCTTGCTTCGATCAgt     >  2:1373403/1‑100 (MQ=255)
                                                                                               gatCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAgt     <  1:1053236/100‑1 (MQ=255)
                                                                                                atCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAgtg    >  1:751268/1‑100 (MQ=255)
                                                                                                  cATTTCGTCAATGCTGCATCCCGCTGAACCTGCGCCAAGGAAGCCGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAGtgtg  <  2:1022990/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GGCATCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAGTGTG  >  NC_000913/1554633‑1554830

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