Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,007,907 +GG intergenic (+109/‑11) pldB → / → b4555 lysophospholipase L(2)/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,007,9071.G92.6% 1875.5 / 130.0 579intergenic (+109/‑11)pldB/b4555lysophospholipase L(2)/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (23/20);  new base G (256/279);  total (279/300)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01
*NC_0009134,007,9072.G92.6% 1854.8 / 130.0 579intergenic (+109/‑11)pldB/b4555lysophospholipase L(2)/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (23/20);  new base G (256/279);  total (279/300)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-01

Only 100 of 625 total aligned reads displayed.

GTTGTTGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGC‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTACCT  >  NC_000913/4007815‑4007996
                                                                                               ||                                                                                         
gttgttGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑Cgcgcgc                                                                                             >  1:8385909/1‑93 (MQ=255)
gttgttGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑Cgcgcgc                                                                                             >  2:1619612/1‑93 (MQ=255)
gttgttGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑Cgcgcgc                                                                                             >  2:7386157/1‑93 (MQ=255)
gttgttGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑Cgcgcgc                                                                                             <  2:4197443/93‑1 (MQ=255)
  tgttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑Cgcgcgcgg                                                                                           <  2:3747792/93‑3 (MQ=255)
  tgttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑Cgcgcgcgg                                                                                           <  2:5426041/93‑3 (MQ=255)
    ttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑Cgcgcgcggca                                                                                         <  1:4019283/93‑5 (MQ=255)
    ttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATACCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑Cgcgggcggca                                                                                         >  2:5317762/1‑89 (MQ=255)
           tGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACtt                                                                                  >  2:5500364/1‑93 (MQ=255)
             aTTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTtg                                                                                <  1:1011489/93‑1 (MQ=255)
              tttAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTtgt                                                                               <  1:2605031/93‑1 (MQ=255)
              tttAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTtgt                                                                               <  2:8152604/93‑1 (MQ=255)
                 aGATGGCGCGTTACTTTCTCTCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGtt                                                                            <  2:7021859/93‑1 (MQ=255)
                 aGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGtt                                                                            >  2:1529165/1‑93 (MQ=255)
                  gATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGttt                                                                           >  1:1348135/1‑93 (MQ=255)
                  gATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGttt                                                                           >  1:7759627/1‑93 (MQ=255)
                    tGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGc                                                                         <  2:3233289/93‑1 (MQ=255)
                    tGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGc                                                                         <  2:8790846/93‑1 (MQ=255)
                        aCGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGAcc                                                                     >  2:8355481/1‑93 (MQ=255)
                         cGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCg                                                                    <  1:2463946/93‑1 (MQ=255)
                           ttACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGgt                                                                  <  2:3075359/93‑1 (MQ=255)
                           ttACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGgt                                                                  >  2:6463319/1‑93 (MQ=255)
                             tcTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGgtcg                                                                <  1:7329643/92‑1 (MQ=255)
                              cTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGgtcgt                                                               <  1:3392760/93‑1 (MQ=255)
                               tttCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGgtcgtc                                                              <  2:1352000/93‑1 (MQ=255)
                               tttCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGgtcgtc                                                              >  1:6166076/1‑93 (MQ=255)
                                   tcCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcacca                                                          <  1:8061239/93‑1 (MQ=255)
                                   tcCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATAAACTTTGTGTTTGCGACCGGTCGTcaccc                                                          >  1:2748392/1‑92 (MQ=255)
                                    cccGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcaccac                                                         >  2:4575490/1‑93 (MQ=255)
                                      cGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGt                                                       <  1:966024/93‑1 (MQ=255)
                                        aCCATACGTTATCCCCTTACGCCAAAGAAGCTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTg                                                     <  2:2768799/93‑1 (MQ=255)
                                         ccATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGa                                                    >  1:4922987/1‑93 (MQ=255)
                                          cATCCGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTCTGCGACCGGTCGTCACCACGTTGAt                                                   <  2:6508531/93‑1 (MQ=255)
                                          cATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAt                                                   >  2:4103000/1‑93 (MQ=255)
                                            tACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtgt                                                 <  1:1669674/93‑1 (MQ=255)
                                             aCGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtgtg                                                <  1:4437590/93‑1 (MQ=255)
                                               gTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTggg                                              <  1:3833992/93‑1 (MQ=255)
                                               gTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTggg                                              <  1:6270518/93‑1 (MQ=255)
                                               cttATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTggg                                              <  1:1834632/92‑1 (MQ=255)
                                                  aTCCCCTGACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGGTGATGTGGGGCa                                           >  2:860159/1‑93 (MQ=255)
                                                    ccccTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCaaa                                         <  2:8417860/93‑1 (MQ=255)
                                                    ccccTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCaaa                                         <  2:6076998/93‑1 (MQ=255)
                                                    ccccTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCAGCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCaaa                                         >  1:4340705/1‑93 (MQ=255)
                                                     cccTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAt                                        <  1:3265994/93‑1 (MQ=255)
                                                     cccTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATAAACTTTGTGTTTGCGACCGTTCGTCACCACGTTGATGTGGGGCAAAt                                        >  1:3133050/1‑93 (MQ=255)
                                                      ccTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAtt                                       <  2:4697637/93‑1 (MQ=255)
                                                      ccTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAtt                                       <  2:3486527/93‑1 (MQ=255)
                                                      ccTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAtt                                       >  1:2921082/1‑93 (MQ=255)
                                                       cTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGC‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAGATTc                                      >  1:7561434/1‑93 (MQ=255)
                                                       cTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTc                                      <  1:3085971/93‑1 (MQ=255)
                                                        ttACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTcg                                     >  1:4615924/1‑93 (MQ=255)
                                                         tACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTcgc                                    >  2:2479675/1‑93 (MQ=255)
                                                         tACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTcgc                                    >  2:7772787/1‑93 (MQ=255)
                                                            gCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGAt                                 <  1:6502601/93‑1 (MQ=255)
                                                             ccAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGCTa                                >  1:3663250/1‑93 (MQ=255)
                                                             ccAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATa                                <  2:7240640/93‑1 (MQ=255)
                                                             ccAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATa                                <  1:4082065/93‑1 (MQ=255)
                                                              cAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATaa                               <  2:7946363/93‑1 (MQ=255)
                                                               aaaGCAACTCTGAAGCTGCGCAC‑‑TGCGCGCTGCATCAGCTTTGTGTTTCAGACCGGACGTCACCACGTTGAAGTGGGGTAAATTCGCGATagt                              >  2:5087814/1‑91 (MQ=255)
                                                               aaaGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCACATTCGCGATAAt                              >  1:7874401/1‑93 (MQ=255)
                                                               aaaGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAAt                              >  1:4225316/1‑93 (MQ=255)
                                                                acgAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATc                             <  1:6839566/91‑1 (MQ=255)
                                                                aaGAAACTCTCAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATc                             <  1:8498011/93‑1 (MQ=255)
                                                                 aGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCt                            >  2:4468458/1‑93 (MQ=255)
                                                                 aGAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCt                            <  2:5627498/93‑1 (MQ=255)
                                                                  gAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTg                           >  1:5949026/1‑93 (MQ=255)
                                                                  gAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTg                           >  2:1047394/1‑93 (MQ=255)
                                                                  gAAACTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTg                           >  1:5581596/1‑93 (MQ=255)
                                                                    aaCTCTGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGga                         <  1:3459565/93‑1 (MQ=255)
                                                                      ctctGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGgaga                       >  1:3722363/1‑93 (MQ=255)
                                                                      ctctGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGgaga                       <  1:6441385/93‑1 (MQ=255)
                                                                      ctctGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGgaga                       >  1:6953407/1‑93 (MQ=255)
                                                                      ctctGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGgaga                       <  1:8213431/93‑1 (MQ=255)
                                                                       tctGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGAt                      <  2:6722846/93‑1 (MQ=255)
                                                                         tGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTa                    >  2:6356553/1‑93 (MQ=255)
                                                                         tGAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTa                    >  1:2329729/1‑93 (MQ=255)
                                                                          gAAGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTaa                   <  1:8691869/93‑1 (MQ=255)
                                                                           aaGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAg                  <  1:1172437/93‑1 (MQ=255)
                                                                            aGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGt                 <  1:8774253/93‑1 (MQ=255)
                                                                            aGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGt                 <  2:2144042/93‑1 (MQ=255)
                                                                            aGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGt                 <  2:257316/93‑1 (MQ=255)
                                                                            aGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGt                 >  1:6645910/1‑93 (MQ=255)
                                                                            aGCTGCTCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGt                 >  2:7679276/1‑93 (MQ=255)
                                                                             gctgctCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTc                <  1:5238976/93‑1 (MQ=255)
                                                                              ctgctCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCt               <  2:4002840/93‑1 (MQ=255)
                                                                               tgctCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCtt              >  1:5042610/1‑93 (MQ=255)
                                                                                gctCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTa             >  2:4970700/1‑93 (MQ=255)
                                                                                gctCAC‑‑CGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTa             >  1:5329444/1‑93 (MQ=255)
                                                                                    cccgcgcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATg         <  1:207734/89‑1 (MQ=255)
                                                                                    accgcgcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGGGATTAAGTCTTACATg         >  1:57602/5‑93 (MQ=255)
                                                                                    accgcgcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATg         >  2:1268089/5‑93 (MQ=255)
                                                                                    accgcgcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATg         >  1:7078986/5‑93 (MQ=255)
                                                                                     ccgcgcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGa        >  2:3892053/4‑93 (MQ=255)
                                                                                      cgcgcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGAt       <  1:1877675/91‑1 (MQ=255)
                                                                                      cgcgcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGAt       <  1:5086501/91‑1 (MQ=255)
                                                                                       gcgcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGAtt      <  1:396042/92‑1 (MQ=255)
                                                                                          cgcgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTAcc   <  2:2256997/93‑1 (MQ=255)
                                                                                          cgcgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTAcc   <  1:19732/93‑1 (MQ=255)
                                                                                          cgcgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGAGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTAcc   <  2:5865301/93‑1 (MQ=255)
                                                                                           gcgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTACCt  <  2:1928218/93‑1 (MQ=255)
                                                                                               ||                                                                                         
GTTGTTGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCAC‑‑CGCGCGC‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTACCT  >  NC_000913/4007815‑4007996

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