Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,201,469 A→G E88G (GAG→GGG)  ygiM → putative signal transduction protein (SH3 domain)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,201,4690AG99.7% 1311.9 / ‑6.9 318E88G (GAG→GGG) ygiMputative signal transduction protein (SH3 domain)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (175/142);  minor base C (1/0);  total (176/142)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 318 total aligned reads displayed.

GTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACACCTGGAATCAGCGCACC  >  NC_000913/3201372‑3201566
                                                                                                 |                                                                                                 
  gACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCaa                                                                                              >  1:5405071/1‑101 (MQ=255)
   aCCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAg                                                                                             <  2:5404885/101‑1 (MQ=255)
    ccTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGc                                                                                            >  2:5405074/1‑101 (MQ=255)
    ccTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGc                                                                                            >  1:5405075/1‑101 (MQ=255)
     cTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGcc                                                                                           <  2:5404562/101‑1 (MQ=255)
      ttattaCAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGccc                                                                                          <  1:5404584/101‑2 (MQ=255)
      ttattaCAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCt                                                                                          <  2:5404583/101‑1 (MQ=255)
       tattaCAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTg                                                                                         >  1:5405081/1‑101 (MQ=255)
         ttaCAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTgcg                                                                                       <  1:5404424/101‑1 (MQ=255)
         ttaCAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTgcg                                                                                       <  2:5404651/101‑1 (MQ=255)
         ttaCAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTgcg                                                                                       <  2:5404458/101‑1 (MQ=255)
         ttaCAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTgcg                                                                                       <  1:5404855/101‑1 (MQ=255)
          taCAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTgcgc                                                                                      <  2:5404388/101‑1 (MQ=255)
          taCAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTgcgc                                                                                      >  2:5405084/1‑101 (MQ=255)
            cAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTc                                                                                    >  1:5405085/1‑101 (MQ=255)
               aCTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCg                                                                                 >  2:5405090/1‑101 (MQ=255)
               aCTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCg                                                                                 >  1:5405091/1‑101 (MQ=255)
                  gACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCgtgt                                                                              >  1:5405097/1‑101 (MQ=255)
                  gACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTACCACTGGGCCAAGCCTGCGCTCCCgtgt                                                                              <  2:5404665/101‑1 (MQ=255)
                  gACGCCAACACCAATTATGCCCAGGTGAAAGAAAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCgtgt                                                                              >  2:5405099/1‑101 (MQ=255)
                    cGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGc                                                                            >  2:5405103/1‑101 (MQ=255)
                      ccAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCa                                                                          <  1:5404561/101‑1 (MQ=255)
                       cAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAg                                                                         <  2:5404377/101‑1 (MQ=255)
                        aaCACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGa                                                                        <  2:5404774/101‑1 (MQ=255)
                          cacCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATc                                                                      >  2:5405111/1‑101 (MQ=255)
                          cacCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATc                                                                      >  1:5405110/1‑101 (MQ=255)
                           acCAATTATGCCCAGGTGGAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTCGGCCAAGCCTGCGCTCCCGTGTGCCAGATCt                                                                     <  1:5404667/101‑1 (MQ=255)
                           acCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCt                                                                     <  2:5404671/101‑1 (MQ=255)
                            ccAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTg                                                                    >  2:5405116/1‑101 (MQ=255)
                              aaTTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGa                                                                  >  1:5405118/1‑101 (MQ=255)
                               aTTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGaa                                                                 <  1:5404632/101‑1 (MQ=255)
                                ttATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGaaa                                                                >  1:5405122/1‑101 (MQ=255)
                                ttATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGaaa                                                                <  2:5404639/101‑1 (MQ=255)
                                 tATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGaaaa                                                               >  2:5405124/1‑101 (MQ=255)
                                  aTGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAAt                                                              >  1:5405126/1‑101 (MQ=255)
                                   tGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATc                                                             <  1:5404504/101‑1 (MQ=255)
                                     cccAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAg                                                           <  1:5404688/101‑1 (MQ=255)
                                     cccAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAg                                                           <  1:5404189/101‑1 (MQ=255)
                                      ccAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAgg                                                          >  2:5405131/1‑101 (MQ=255)
                                        aGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTc                                                        >  1:5405134/1‑101 (MQ=255)
                                         ggTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCa                                                       <  1:5404701/101‑1 (MQ=255)
                                          gTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCaa                                                      >  2:5405138/1‑101 (MQ=255)
                                           tGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCaaa                                                     >  1:5405139/1‑101 (MQ=255)
                                            gAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCaaaa                                                    >  1:5405147/1‑101 (MQ=255)
                                            gAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCaaaa                                                    >  2:5405145/1‑101 (MQ=255)
                                             aaaGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAAc                                                   >  2:5405149/1‑101 (MQ=255)
                                              aaGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAAcc                                                  >  1:5405154/1‑101 (MQ=255)
                                               aGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAAccc                                                 <  1:5404648/101‑1 (MQ=255)
                                               aGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAAccc                                                 >  2:5405155/1‑101 (MQ=255)
                                                gACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCt                                                <  2:5404783/101‑1 (MQ=255)
                                                 aCAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTg                                               <  2:5404415/101‑1 (MQ=255)
                                                   aGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGAc                                             >  2:5405159/1‑101 (MQ=255)
                                                    gCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGAcc                                            <  2:5404754/101‑1 (MQ=255)
                                                     ctctGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCg                                           <  1:5404771/101‑1 (MQ=255)
                                                      tctGGCCGTACCGCCTTGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGa                                          >  1:5405164/1‑101 (MQ=255)
                                                       ctGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGAt                                         >  2:5405167/1‑101 (MQ=255)
                                                         ggCCGTACCGTCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATaa                                       <  2:5404524/101‑1 (MQ=255)
                                                         ggCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATaa                                       >  2:5405170/1‑101 (MQ=255)
                                                         ggCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATaa                                       <  1:5404720/101‑1 (MQ=255)
                                                         ggCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATaa                                       <  2:5404624/101‑1 (MQ=255)
                                                          gCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATaaa                                      <  2:5404602/101‑1 (MQ=255)
                                                           ccGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAAc                                     <  2:5404686/101‑1 (MQ=255)
                                                           ccGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAAc                                     >  2:5405173/1‑101 (MQ=255)
                                                           ccGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAAc                                     <  1:5405042/101‑1 (MQ=255)
                                                            cGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACt                                    >  2:5405176/1‑101 (MQ=255)
                                                             gTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTc                                   >  1:5405177/1‑101 (MQ=255)
                                                              tACCGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCa                                  >  2:5405180/1‑101 (MQ=255)
                                                                ccGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCAcc                                <  2:5404617/101‑1 (MQ=255)
                                                                ccGCCTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCAc                                 >  2:5405184/1‑100 (MQ=255)
                                                                    cTGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAata                            <  2:5404606/101‑1 (MQ=255)
                                                                     tGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAatat                           >  2:5405190/1‑101 (MQ=255)
                                                                     tGGATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAatat                           >  1:5405189/1‑101 (MQ=255)
                                                                       gATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCg                         <  2:5405022/101‑1 (MQ=255)
                                                                       gATCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCg                         >  2:5405193/1‑101 (MQ=255)
                                                                        aTCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGa                        >  2:5405197/1‑101 (MQ=255)
                                                                         tCCCGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGAt                       >  1:5405202/1‑101 (MQ=255)
                                                                          cccGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATa                      >  2:5405204/1‑101 (MQ=255)
                                                                          cccGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATa                      <  1:5404996/101‑1 (MQ=255)
                                                                           ccGTTGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATaa                     >  1:5405206/1‑101 (MQ=255)
                                                                              ttGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAacac                  >  1:5405212/1‑101 (MQ=255)
                                                                              ttGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAacac                  <  1:5404589/101‑1 (MQ=255)
                                                                              ttGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAacac                  >  2:5405211/1‑101 (MQ=255)
                                                                               tGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACAcc                 >  1:5405217/1‑101 (MQ=255)
                                                                               tGAAACAACTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACAcc                 <  1:5404748/101‑1 (MQ=255)
                                                                                   acaacTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACACCTGGa             <  2:5404797/101‑1 (MQ=255)
                                                                                   acaacTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACACCTGGa             <  2:5404719/101‑1 (MQ=255)
                                                                                    caacTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACACCTGGaa            >  2:5405221/1‑101 (MQ=255)
                                                                                     aacTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACACCTGGAAt           >  1:5405222/1‑101 (MQ=255)
                                                                                      acTTAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACACCTGGAATc          >  2:5405226/1‑101 (MQ=255)
                                                                                        ttAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACACCTGGAATCAg        <  2:5404843/101‑1 (MQ=255)
                                                                                        ttAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGCCAAAActc                                                 >  2:5405229/1‑58 (MQ=255)
                                                                                         tAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACACCTGGAATCAgc       >  1:5405231/1‑101 (MQ=255)
                                                                                         tAGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACACCTGGAATCAgc       <  1:5404540/101‑1 (MQ=255)
                                                                                          aGCACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGACCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACACCTGGAATCAgcg      >  2:5405234/1‑101 (MQ=255)
                                                                                            cACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCGCCAATATCGATAaca                   <  1:5404732/86‑1 (MQ=255)
                                                                                            cACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACACCTGGAATCAGCGCa    <  2:5404903/101‑1 (MQ=255)
                                                                                            cACTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACACCTGGAATCAGCGCa    <  1:5404495/101‑1 (MQ=255)
                                                                                             aCTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACACCTGGAATCAGCGCAc   <  1:5405103/101‑1 (MQ=255)
                                                                                              cTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACACCTGGAATCAGCGCAcc  <  1:5404877/101‑1 (MQ=255)
                                                                                              cTGGGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACACCCGGAATCAGCGCAcc  >  2:5405239/1‑101 (MQ=255)
                                                                                                 |                                                                                                 
GTGACCTTATTACAAACTGACGCCAACACCAATTATGCCCAGGTGAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACACCTGGAATCAGCGCACC  >  NC_000913/3201372‑3201566

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: ‑