Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,550,930 A→G V515A (GTG→GCG)  malP ← maltodextrin phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,550,9300AG98.9% 1164.9 / 0.5 280V515A (GTG→GCG) malPmaltodextrin phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (135/142);  total (137/143)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

Only 100 of 283 total aligned reads displayed.

GCAGCAGATTCAGGTGCTGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCACAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAATTTTTCCAGATTGATCAGCTG  >  NC_000913/3550832‑3551027
                                                                                                  |                                                                                                 
gcagcaGATTCAGGTGCTGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTcgca                                                                                                 <  2:5209850/101‑1 (MQ=255)
gcagcaGATTCAGGTGCTGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTcgca                                                                                                 >  1:5210213/1‑101 (MQ=255)
  agcagATTCAGGTGCTGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCaaa                                                                                               <  1:5209818/101‑1 (MQ=255)
  agcagATTCAGGTGCTGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCaaa                                                                                               >  1:5210216/1‑101 (MQ=255)
      gATTCAGGTGCTGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTcc                                                                                           >  1:5210222/1‑101 (MQ=255)
      gATTCAGGTGCTGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTcc                                                                                           >  2:5210221/1‑101 (MQ=255)
       aTTCAGGTGCTGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCg                                                                                          >  1:5210225/1‑101 (MQ=255)
        ttCAGGTGCTGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGc                                                                                         >  1:5210228/1‑101 (MQ=255)
         tCAGGTGCTGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGc                                                                                         <  1:5209721/100‑1 (MQ=255)
          cAGGTGCTGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCa                                                                                       >  2:5210232/1‑101 (MQ=255)
            ggTGCTGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGa                                                                                     >  1:5210235/1‑101 (MQ=255)
              tGCTGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACg                                                                                   >  2:5210237/1‑101 (MQ=255)
                cTGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGa                                                                                 >  1:5210240/1‑101 (MQ=255)
                cTGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGa                                                                                 <  1:5209834/101‑1 (MQ=255)
                cTGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGa                                                                                 <  1:5209638/101‑1 (MQ=255)
                 tGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGAc                                                                                <  2:5209725/101‑1 (MQ=255)
                 tGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGAc                                                                                >  2:5210243/1‑101 (MQ=255)
                   gCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACtt                                                                              >  2:5210246/1‑101 (MQ=255)
                    cGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACttt                                                                             <  2:5209681/101‑1 (MQ=255)
                     gTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTa                                                                            >  1:5210249/1‑101 (MQ=255)
                        tGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTc                                                                         >  1:5210253/1‑101 (MQ=255)
                        tGTACTCATGCAAACGTTTGATCTGAATACCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTc                                                                         <  1:5209986/101‑1 (MQ=255)
                         gTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCg                                                                        <  1:5209838/101‑1 (MQ=255)
                          tACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGc                                                                       >  1:5210256/1‑101 (MQ=255)
                           aCTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGcc                                                                      >  2:5210257/1‑101 (MQ=255)
                            cTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCt                                                                     <  2:5210020/101‑1 (MQ=255)
                             tCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTg                                                                    >  1:5210260/1‑101 (MQ=255)
                              cATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGc                                                                   <  1:5209943/101‑1 (MQ=255)
                              cATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGc                                                                   <  1:5209973/101‑1 (MQ=255)
                               aTGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCt                                                                     >  1:5210264/1‑98 (MQ=255)
                                tGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGa                                                                                 >  2:5210266/1‑85 (MQ=255)
                                tGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCtt                                                                 <  1:5209927/101‑1 (MQ=255)
                                 gCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTg                                                                <  1:5209853/101‑1 (MQ=255)
                                 gCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTg                                                                <  2:5209906/101‑1 (MQ=255)
                                   aaaCGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGAt                                                              >  1:5210271/1‑101 (MQ=255)
                                    aaCGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGAtc                                                             >  1:5210274/1‑101 (MQ=255)
                                     aCGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGAtct                                                            >  2:5210276/1‑101 (MQ=255)
                                     aCGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGAtct                                                            <  1:5209762/101‑1 (MQ=255)
                                      cGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGAtctc                                                           <  1:5209889/101‑1 (MQ=255)
                                      cGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGAtctc                                                           <  1:5209735/101‑1 (MQ=255)
                                      cGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGAtctc                                                           <  1:5209750/101‑1 (MQ=255)
                                       gTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTcg                                                          >  2:5210280/1‑101 (MQ=255)
                                       gTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTcg                                                          <  2:5209883/101‑1 (MQ=255)
                                        tttGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTcgc                                                         >  2:5210283/1‑101 (MQ=255)
                                        tttGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTcgc                                                         <  1:5209868/101‑1 (MQ=255)
                                        tttGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTcgc                                                         >  1:5210282/1‑101 (MQ=255)
                                           gATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGat                                                      >  1:5210286/1‑101 (MQ=255)
                                            aTCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGata                                                     <  2:5209575/101‑1 (MQ=255)
                                             tCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGatat                                                    >  1:5210288/1‑101 (MQ=255)
                                              cTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATAtt                                                   <  2:5209969/101‑1 (MQ=255)
                                              cTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATAtt                                                   <  2:5209934/101‑1 (MQ=255)
                                              cTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATAtt                                                   <  1:5209825/101‑1 (MQ=255)
                                                  atatCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTg                                               <  1:5209938/101‑1 (MQ=255)
                                                   tatCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGa                                              >  1:5210292/1‑101 (MQ=255)
                                                    atCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGAc                                             <  1:5209879/101‑1 (MQ=255)
                                                    atCGAAAAACGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTc                                                                                            <  2:5209457/54‑1 (MQ=39)
                                                     tCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACg                                            >  1:5210297/1‑101 (MQ=255)
                                                     tCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACg                                            >  2:5210296/1‑101 (MQ=255)
                                                        aaaaTCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGaa                                         <  1:5209579/101‑1 (MQ=255)
                                                        aaaaTCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGaa                                         >  1:5210302/1‑101 (MQ=255)
                                                         aaaTCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAAt                                        <  1:5210042/101‑1 (MQ=255)
                                                         aaaTCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAAt                                        <  2:5209771/101‑1 (MQ=255)
                                                          aaTCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAAtt                                       >  2:5210309/1‑101 (MQ=255)
                                                          aaTCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAAtt                                       >  2:5210305/1‑101 (MQ=255)
                                                          aaTCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAAtt                                       >  1:5210306/1‑101 (MQ=255)
                                                           aTCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAAttt                                      <  1:5209803/101‑1 (MQ=255)
                                                            tCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTc                                     <  2:5209621/101‑1 (MQ=255)
                                                            tCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTc                                     <  2:5209898/101‑1 (MQ=255)
                                                              gCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGc                                   <  2:5210095/101‑1 (MQ=255)
                                                              gCCTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGc                                   >  1:5210312/1‑101 (MQ=255)
                                                               ccTGTGGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGca                                  >  1:5210315/1‑101 (MQ=255)
                                                                  gtgGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGcatca                               <  2:5209786/101‑1 (MQ=255)
                                                                   tgGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGcatcat                              <  1:5209699/101‑1 (MQ=255)
                                                                   tgGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGcatcat                              >  1:5210319/1‑101 (MQ=255)
                                                                   tgGATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGcatcat                              <  1:5209459/101‑1 (MQ=255)
                                                                    ggATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCatcatc                             >  2:5210321/1‑101 (MQ=255)
                                                                    ggATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCatcatc                             <  1:5209961/101‑1 (MQ=255)
                                                                     gATTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCAtcatca                            <  2:5210002/101‑1 (MQ=255)
                                                                      aTTGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAg                           <  1:5210006/101‑1 (MQ=255)
                                                                        tGATCTCAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCg                         >  1:5210324/1‑101 (MQ=255)
                                                                           tctcAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAAt                      >  2:5210332/1‑101 (MQ=255)
                                                                           tctcAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAAt                      >  2:5210329/1‑101 (MQ=255)
                                                                           tctcAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAAt                      >  1:5210328/1‑101 (MQ=255)
                                                                            ctcAATACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAAtt                     <  1:5209913/101‑1 (MQ=255)
                                                                               aaTACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAAttttt                  >  1:5210335/1‑101 (MQ=255)
                                                                                aTACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAATTTTTc                 <  1:5210126/101‑1 (MQ=255)
                                                                                aTACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATGCGCCTGCTTGATCTCGCGATATTGCTGAc                                             <  2:5209872/73‑1 (MQ=255)
                                                                                 tACCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAATTTTTcc                >  2:5210337/1‑101 (MQ=255)
                                                                                  aCCGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAATTTTTCCa               >  2:5210340/1‑101 (MQ=255)
                                                                                   ccGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCtt                                                                 <  2:5209794/50‑1 (MQ=255)
                                                                                   ccGGTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAAt                      <  1:5210049/93‑1 (MQ=255)
                                                                                    cgttACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAATTTTTCCAGa             >  1:5210343/4‑101 (MQ=255)
                                                                                     ggTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAATTTTTCCAGAt            <  2:5210178/101‑1 (MQ=255)
                                                                                      gTACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAATTTTTCCAGAtt           <  1:5209952/101‑1 (MQ=255)
                                                                                       tACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAATTTTTCCAGATTg          <  1:5209989/101‑1 (MQ=255)
                                                                                       tACGAACTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAATTTTTCCAGATTg          >  2:5210345/1‑101 (MQ=255)
                                                                                           aaCTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAATTTTTCCAGATTGATCa      <  2:5209962/101‑1 (MQ=255)
                                                                                            aCTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAATTTTTCCAGATTGATCAg     <  1:5209787/101‑1 (MQ=255)
                                                                                             cTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAATTTTTCCAGATTGATCAGc    <  1:5209857/101‑1 (MQ=255)
                                                                                             cTTTCGCAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAATTTTTCCAGATTGATCAGc    >  2:5210348/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
GCAGCAGATTCAGGTGCTGGCGTTTGTACTCATGCAAACGTTTGATCTGAATATCGAAAATCGCCTGTGGATTGATCTCAATACCGGTACGAACTTTCACAAACTCCGCCAGACGGACTTTATTCGCCTGCTTGATCTCGCGATATTGCTGACGGAATTTCGCATCATCAGCGAATTTTTCCAGATTGATCAGCTG  >  NC_000913/3550832‑3551027

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