Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 79,855 A→G V337A (GTG→GCG)  leuC ← 3‑isopropylmalate isomerase (dehydratase) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091379,8550AG99.4% 475.0 / NA 154V337A (GTG→GCG) leuC3‑isopropylmalate isomerase (dehydratase) subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (80/73);  minor base T (0/1);  total (80/74)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 159 total aligned reads displayed.

ACTTTTCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGATCGGCAAACGAAGCC  >  NC_000913/79765‑79944
                                                                                          |                                                                                         
aCTTTTCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCt                                                                                         <  2:4214313/93‑1 (MQ=255)
aCTTTTCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCt                                                                                         >  1:1145207/1‑93 (MQ=255)
aCTTTTCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCt                                                                                         >  1:1941158/1‑93 (MQ=255)
 cTTTTCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCtt                                                                                        >  1:3095842/1‑93 (MQ=255)
 cTTTTCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCtt                                                                                        >  1:433067/1‑93 (MQ=255)
  ttttCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTc                                                                                       >  2:2533188/1‑93 (MQ=255)
   tttCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCg                                                                                      <  1:4153569/93‑1 (MQ=255)
      cGCTCTTTGGCGATCTCCGCTGCCGATCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTc                                                                                   <  1:2330304/93‑1 (MQ=255)
       gctgtttGGCGACCTCCGCTGCCGCGCGTAAATCTCCAATGCGGGAGTTGGTACAGCAACCGATAAACGCTTTGTCGATAGCCGCTTCGGTCa                                                                                  <  2:44250/89‑1 (MQ=255)
       gCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCa                                                                                  >  1:2061710/1‑93 (MQ=255)
         ccTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGc                                                                                >  1:4157706/1‑93 (MQ=255)
           tttGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCgg                                                                              >  1:4280639/1‑93 (MQ=255)
            ttGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGa                                                                             >  2:1286559/1‑93 (MQ=255)
            ttGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGACGCTTCGGTCAGCGGa                                                                             >  1:1300009/1‑93 (MQ=255)
            tcggCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGa                                                                             <  2:2574128/91‑1 (MQ=255)
             tGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGaa                                                                            <  2:1145207/93‑1 (MQ=255)
             tGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGaa                                                                            <  1:3279382/93‑1 (MQ=255)
              ggCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAAt                                                                           <  1:4252800/93‑1 (MQ=255)
              ggCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAAt                                                                           >  2:2013116/1‑93 (MQ=255)
              ggCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAAt                                                                           <  2:940038/93‑1 (MQ=255)
                cGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATAc                                                                         <  1:3165675/93‑1 (MQ=255)
                cGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATAc                                                                         <  2:63192/93‑1 (MQ=255)
                    ctcCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCgg                                                                     <  1:653440/93‑1 (MQ=255)
                    ctcCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCgg                                                                     >  2:4025007/1‑93 (MQ=255)
                    ctcCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCgg                                                                     >  2:3704632/1‑93 (MQ=255)
                     tcCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGt                                                                    >  2:2115870/1‑93 (MQ=255)
                     tcCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGt                                                                    >  2:410284/1‑93 (MQ=255)
                     tcCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGt                                                                    >  2:2879019/1‑93 (MQ=255)
                        gCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTc                                                                 <  2:235601/93‑1 (MQ=255)
                        gCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTc                                                                 >  2:3553957/1‑93 (MQ=255)
                         cTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCa                                                                <  2:3024712/93‑1 (MQ=255)
                          tGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAg                                                               >  2:2601821/1‑93 (MQ=255)
                           gCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGc                                                              <  2:35896/93‑1 (MQ=255)
                           gCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGc                                                              >  2:3888509/1‑93 (MQ=255)
                             cgcgcgTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGccc                                                            <  1:331250/93‑1 (MQ=255)
                             cgcgcgTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGccc                                                            >  1:3901481/1‑93 (MQ=255)
                               cgcgTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCa                                                          <  2:2012748/93‑1 (MQ=255)
                               cgcgTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCa                                                          <  1:2765692/93‑1 (MQ=255)
                                gcgTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCat                                                         >  1:2818610/1‑93 (MQ=255)
                                 cgTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCata                                                        <  2:1751208/93‑1 (MQ=255)
                                 cgTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCata                                                        <  1:4536950/93‑1 (MQ=255)
                                 cgTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCata                                                        <  1:3515313/93‑1 (MQ=255)
                                  gTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCatat                                                       >  2:550580/1‑93 (MQ=255)
                                   tAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCAtata                                                      >  2:1459425/1‑93 (MQ=255)
                                   tAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCAtata                                                      >  1:26615/1‑93 (MQ=255)
                                     aaTCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAgg                                                    >  1:3730200/1‑93 (MQ=255)
                                      aTCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGc                                                   >  1:181935/1‑93 (MQ=255)
                                       tCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCTCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGcc                                                  <  1:3052558/93‑1 (MQ=255)
                                        tttCAATGCGCTAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCa                                                 <  2:4445124/92‑1 (MQ=255)
                                         ttCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAg                                                <  1:1594903/93‑1 (MQ=255)
                                         ttCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAg                                                >  1:2253363/1‑93 (MQ=255)
                                          tCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAgc                                               >  2:3964702/1‑93 (MQ=255)
                                             aTGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCt                                            <  1:3172098/93‑1 (MQ=255)
                                               gcgcGAGTTGGTACAGAAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATGCCCGGTTTCAGCCCCATATAGGCCAGCGCttt                                          <  1:2143508/93‑1 (MQ=255)
                                                cgcgAGGTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCtttt                                         <  2:451575/93‑1 (MQ=255)
                                                 gcgAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCttttt                                        >  1:528522/1‑93 (MQ=255)
                                                 gcgAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCttttt                                        <  2:736943/93‑1 (MQ=255)
                                                 gcgAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCttttt                                        >  2:3217330/1‑93 (MQ=255)
                                                  cgAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTc                                       <  1:650268/93‑1 (MQ=255)
                                                  cgAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTc                                       <  1:2501706/93‑1 (MQ=255)
                                                  cgAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTc                                       <  2:1418837/93‑1 (MQ=255)
                                                  cgAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTc                                       <  2:4543508/93‑1 (MQ=255)
                                                   gAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCt                                      >  2:207123/1‑93 (MQ=255)
                                                   gAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCt                                      <  2:3156819/93‑1 (MQ=255)
                                                   gAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCt                                      <  2:10873/93‑1 (MQ=255)
                                                   gAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCt                                      >  1:1258995/1‑93 (MQ=255)
                                                       tGGTACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCg                                  <  2:3480928/93‑1 (MQ=255)
                                                       tGGTACAGGAACAGATAAACACTTGGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCg                                  <  1:47814/93‑1 (MQ=255)
                                                          tACAGGAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGAcg                               >  1:3462249/1‑93 (MQ=255)
                                                              ggAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGAcgcgcg                           >  1:4171740/1‑93 (MQ=255)
                                                              ggAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGAcgcgcg                           >  2:3163111/1‑93 (MQ=255)
                                                               gAACCGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATAACCGGTTTCAGCCCCATATAGTCCAGCGCTTTTTATTCCGACGCGCGt                          >  1:37477/1‑93 (MQ=255)
                                                                   cGATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCaa                      >  2:3901688/1‑93 (MQ=255)
                                                                    gATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAAc                     >  1:3438386/1‑93 (MQ=255)
                                                                    gATAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAAc                     <  2:3422111/93‑1 (MQ=255)
                                                                     aTAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAAcc                    <  1:84858/93‑1 (MQ=255)
                                                                     aTAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAAcc                    >  2:2988202/1‑93 (MQ=255)
                                                                      tAAACACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCg                   <  1:1025348/93‑1 (MQ=255)
                                                                      tAAACACTTTGTCGATAGACGCTTCGGTCAGCGGAATACCCGGGTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCg                   >  2:185059/1‑93 (MQ=255)
                                                                      tAAACACTTCGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCg                   <  2:3763957/93‑1 (MQ=255)
                                                                        aaCACTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGa                 >  1:1591735/1‑93 (MQ=255)
                                                                         acacTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGAt                <  1:217835/93‑1 (MQ=255)
                                                                         acacTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGAt                <  2:3744029/93‑1 (MQ=255)
                                                                          cacTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGATc               <  2:4234240/93‑1 (MQ=255)
                                                                           acTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGATCg              >  1:2245187/1‑93 (MQ=255)
                                                                           acTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGATCg              >  2:2395325/1‑93 (MQ=255)
                                                                           acTTTGTCGATAGACGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCCTATAGGCCCGCGCTTTTTCTGCCGACGCGCGTTCAACCCGATCg              >  1:302821/1‑93 (MQ=255)
                                                                            cTTTGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGATCgg             >  2:751671/1‑93 (MQ=255)
                                                                              ttGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGATCGGCa           >  1:1033948/1‑93 (MQ=255)
                                                                               tGTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGATCGGCaa          >  1:3351156/1‑93 (MQ=255)
                                                                                gTCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGATCGGCaaa         <  1:1505397/93‑1 (MQ=255)
                                                                                 tCGATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGATCGGCAAAc        >  1:2593656/1‑93 (MQ=255)
                                                                                 tCGATAGCCGCGTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGATCGGCAAAc        <  2:4661541/93‑1 (MQ=255)
                                                                                  cGATAGCCGCTTCGGTCAGCGGAATCCCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGATCGGCAAACg       <  2:3385711/93‑1 (MQ=255)
                                                                                   gATAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGATCGGCAAACGa      <  1:2453477/93‑1 (MQ=255)
                                                                                    aTAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGATCGGCAAACGaa     <  1:1150354/93‑1 (MQ=255)
                                                                                     tAGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGATCGGCAAACGAAg    >  2:3651151/1‑93 (MQ=255)
                                                                                      tgCCGCTTCGGTCAGCGTAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCACCCGGATCGGCAAACGAAGc   <  2:567212/92‑1 (MQ=255)
                                                                                      aGCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGATCGGCAAACGAAGc   >  1:4330287/1‑93 (MQ=255)
                                                                                       gCCGCTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGATCGGCAAACGAAGcc  >  1:386650/1‑93 (MQ=255)
                                                                                          |                                                                                         
ACTTTTCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAGCGCTTTTTCTGCCGACGCGCGTTCAACCGGATCGGCAAACGAAGCC  >  NC_000913/79765‑79944

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