Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,472,426 (A)8→7 coding (108/135 nt) pyrL ← pyrBI operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,472,4190A.99.1% 1641.1 / 5.9 319coding (115/135 nt)pyrLpyrBI operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base . (213/103);  total (213/106)
Fisher's exact test for biased strand distribution p-value = 3.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 363 total aligned reads displayed.

GCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATG  >  NC_000913/4472319‑4472515
                                                                                                    |                                                                                                
gCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCa                                                                                                  <  2:7381050/101‑1 (MQ=255)
gCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCa                                                                                                  >  2:7381504/1‑101 (MQ=255)
gCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCa                                                                                                  <  2:7381150/101‑1 (MQ=255)
  cAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCaaa                                                                                                >  2:7381508/1‑101 (MQ=255)
  cAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCaaa                                                                                                <  1:7381061/101‑1 (MQ=255)
   aGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCaaaa                                                                                               >  1:7381511/1‑101 (MQ=255)
    gCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCaaaaa                                                                                              <  1:7381076/101‑1 (MQ=255)
     cACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCaaaaaa                                                                                             >  2:7381516/1‑101 (MQ=255)
     cACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCaaaaaa                                                                                             >  1:7381518/1‑101 (MQ=255)
      aCCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCaaaaaaa                                                                                            >  2:7381520/1‑101 (MQ=255)
       ccAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCaaaaaaag                                                                                           <  1:7381153/101‑2 (MQ=255)
          gATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCaaaaaaagccc                                                                                        >  1:7381525/1‑97 (MQ=255)
            ttAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCt                                                                                     <  2:7380961/101‑1 (MQ=255)
             tAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTc                                                                                    <  1:7381090/101‑1 (MQ=255)
             tAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTc                                                                                    <  2:7380998/101‑1 (MQ=255)
              aaGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCg                                                                                   >  1:7381531/1‑101 (MQ=255)
               aGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGa                                                                                  >  2:7381534/1‑101 (MQ=255)
                  tCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTg                                                                               <  1:7381135/101‑1 (MQ=255)
                  tCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTg                                                                               <  2:7381035/101‑1 (MQ=255)
                   cATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGa                                                                              >  2:7381538/1‑101 (MQ=255)
                    aTCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAg                                                                             >  2:7381543/1‑101 (MQ=255)
                    aTCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAg                                                                             <  2:7381108/101‑1 (MQ=255)
                      cgcgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAggg                                                                           <  2:7381221/101‑1 (MQ=255)
                       gcgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAgggg                                                                          >  2:7381547/1‑101 (MQ=255)
                        cgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGc                                                                         >  1:7381549/1‑101 (MQ=255)
                        cgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGc                                                                         <  2:7381254/101‑1 (MQ=255)
                         gACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCa                                                                                                  >  1:7381557/1‑76 (MQ=255)
                         gACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCt                                                                        >  2:7381554/1‑101 (MQ=255)
                           cTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCaaaaaaa                                                                                            >  1:7381563/1‑80 (MQ=255)
                            tAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTggg                                                                     >  2:7381564/1‑101 (MQ=255)
                              aGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGaa                                                                   >  2:7381568/1‑101 (MQ=255)
                                gTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATg                                                                 <  1:7381186/101‑1 (MQ=255)
                                gTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATg                                                                 <  1:7381023/101‑1 (MQ=255)
                                 tCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATgg                                                                >  1:7381575/1‑101 (MQ=255)
                                   gTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGt                                                              >  2:7381576/1‑101 (MQ=255)
                                     ttATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATggg                                                               >  2:7381584/1‑98 (MQ=255)
                                     ttATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGa                                                            >  1:7381582/1‑101 (MQ=255)
                                       aTGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATc                                                          >  2:7381590/1‑101 (MQ=255)
                                       aTGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATc                                                          >  1:7381591/1‑101 (MQ=255)
                                        tGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCa                                                         >  2:7381595/1‑101 (MQ=255)
                                          gAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAAc                                                       <  2:7381075/101‑1 (MQ=255)
                                           aaaTGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACg                                                      >  2:7381602/1‑101 (MQ=255)
                                           aaaTGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACg                                                      >  1:7381604/1‑101 (MQ=255)
                                            aaTGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACgg                                                     <  2:7381425/101‑1 (MQ=255)
                                             aTGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACggg                                                    >  1:7381609/1‑101 (MQ=255)
                                              tGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGAGGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTTGGAATGGGTGa                                                            <  2:7380863/92‑1 (MQ=255)
                                               gATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCgct                                                                                 >  1:7381614/1‑68 (MQ=255)
                                               gATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaa                                                  <  2:7381295/101‑1 (MQ=255)
                                               gATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaa                                                  >  2:7381612/1‑101 (MQ=255)
                                                 tatGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaaga                                                >  2:7381618/1‑101 (MQ=255)
                                                 tatGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaaga                                                <  2:7381194/101‑1 (MQ=255)
                                                  atGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaagaa                                               >  2:7381622/1‑101 (MQ=255)
                                                    gTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTAAGGGGCTGGGAATGGGTGATCAACGGGAAGaaaa                                             <  1:7380788/101‑1 (MQ=255)
                                                     tttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaaaa                                            <  1:7381327/101‑1 (MQ=255)
                                                     tttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaaaa                                            >  1:7381626/1‑101 (MQ=255)
                                                      ttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAAc                                           >  2:7381633/1‑101 (MQ=255)
                                                      ttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAAc                                           >  1:7381631/1‑101 (MQ=255)
                                                        cTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACgg                                         >  2:7381640/1‑101 (MQ=255)
                                                        cTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACgg                                         <  2:7381185/101‑1 (MQ=255)
                                                        cTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACgg                                         >  1:7381638/1‑101 (MQ=255)
                                                         tGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGc                                        <  1:7381211/101‑1 (MQ=255)
                                                          gATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCa                                       >  2:7381646/1‑101 (MQ=255)
                                                          gATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCACCTCGATTGAgg                                                                            >  1:7381648/1‑64 (MQ=255)
                                                            tataGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAgg                                     >  2:7381651/1‑101 (MQ=255)
                                                             ataGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGc                                    >  1:7381655/1‑101 (MQ=255)
                                                              taGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGcc                                   >  2:7381660/1‑101 (MQ=255)
                                                               aGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGcccc                                                                                      >  1:7381667/1‑49 (MQ=39)
                                                               aGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCa                                  <  2:7380923/101‑1 (MQ=255)
                                                               aGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCa                                  <  2:7381364/101‑1 (MQ=255)
                                                               aGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCa                                  >  2:7381666/1‑101 (MQ=255)
                                                                 cGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGc                                <  1:7380978/101‑1 (MQ=255)
                                                                  ggATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCg                               >  2:7381671/1‑101 (MQ=255)
                                                                   gATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGt                              >  1:7381675/1‑101 (MQ=255)
                                                                     ttAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGGAt                                                                  >  1:7381679/1‑63 (MQ=38)
                                                                      tAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtt                           >  2:7381681/1‑101 (MQ=255)
                                                                        gCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGcccc                                                                                      >  1:7381686/1‑40 (MQ=38)
                                                                         ccATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCttttt                        >  2:7381688/1‑101 (MQ=255)
                                                                          cATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGtttttttt                       <  1:7380910/101‑1 (MQ=255)
                                                                          cATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtttttt                       >  1:7381691/1‑101 (MQ=255)
                                                                           aTCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGGCAGCg                               >  2:7381697/1‑92 (MQ=255)
                                                                            tCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCa                     <  1:7381296/101‑1 (MQ=255)
                                                                             cTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAg                    >  1:7381699/1‑101 (MQ=255)
                                                                              ttttATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCa                     <  1:7381160/99‑1 (MQ=255)
                                                                               tttATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAc                  >  1:7381707/1‑101 (MQ=255)
                                                                               tttATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAc                  >  2:7381706/1‑101 (MQ=255)
                                                                                ttATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAg                                 >  2:7381713/1‑85 (MQ=255)
                                                                                ttATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAcg                 <  2:7381087/101‑1 (MQ=255)
                                                                                ttATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAcg                 >  1:7381711/1‑101 (MQ=255)
                                                                                 tATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAcgc                >  1:7381717/1‑101 (MQ=255)
                                                                                   tctcCTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCgg              >  2:7381721/1‑101 (MQ=255)
                                                                                      ccTGACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTaa           >  1:7381726/1‑101 (MQ=255)
                                                                                         gACGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGAc        >  1:7381730/1‑101 (MQ=255)
                                                                                          aCGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACa       <  1:7381363/101‑1 (MQ=255)
                                                                                          aCGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACa       >  2:7381735/1‑101 (MQ=255)
                                                                                           cGCCTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaa      >  1:7381736/1‑101 (MQ=255)
                                                                                              cTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCggga            >  2:7381740/1‑90 (MQ=255)
                                                                                              cTGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAAt   <  2:7381202/101‑1 (MQ=255)
                                                                                               tGGGCAAAAAAAGCCCCTCGATTGAGTGGCTGGGAATGGGTGATCa                                                         <  1:7381241/46‑1 (MQ=37)
                                                                                               tGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCg                               >  2:7381746/1‑72 (MQ=255)
                                                                                               tGGGCAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATg  >  1:7381743/1‑101 (MQ=255)
                                                                                                    |                                                                                                
GCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATG  >  NC_000913/4472319‑4472515

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