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.0% 1063.6 / 85.8 363intergenic (+109/‑11)pldB/b4555lysophospholipase L(2)/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (8/21);  new base G (179/153);  total (188/175)
Fisher's exact test for biased strand distribution p-value = 1.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01
*NC_0009134,007,9072.G92.0% 1051.6 / 85.8 363intergenic (+109/‑11)pldB/b4555lysophospholipase L(2)/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (8/21);  new base G (180/152);  total (188/175)
Fisher's exact test for biased strand distribution p-value = 6.53e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

Only 100 of 406 total aligned reads displayed.

GTTGTTGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGC‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTACCT  >  NC_000913/4007815‑4007996
                                                                                               ||                                                                                         
gttgttGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑gcgcgc                                                                                             >  1:5811177/1‑93 (MQ=255)
gttgttGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑gcgcgc                                                                                             >  2:1372485/1‑93 (MQ=255)
 ttgttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑gcgcgcg                                                                                            >  1:7740754/1‑92 (MQ=255)
 ttgttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑gcgcgcg                                                                                            <  2:2156028/93‑2 (MQ=255)
  tgttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑gcgcgcgg                                                                                           >  1:4636122/1‑91 (MQ=255)
  tgttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑gcgcgcgg                                                                                           >  2:3127180/1‑91 (MQ=255)
   gttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑gcgcgcggc                                                                                          >  2:149315/1‑90 (MQ=255)
   gttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑gcgcgcggc                                                                                          >  2:4992364/1‑90 (MQ=255)
    ttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑gcgcgcggca                                                                                         <  2:6239479/93‑5 (MQ=255)
           tGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACtt                                                                                  <  1:2551115/93‑1 (MQ=255)
            gATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACttt                                                                                 >  1:1566620/1‑93 (MQ=255)
            gATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACttt                                                                                 >  2:3833345/1‑93 (MQ=255)
             aTTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTtg                                                                                >  2:1088504/1‑93 (MQ=255)
              tttAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTtgt                                                                               >  2:5869611/1‑93 (MQ=255)
              tttAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTtgt                                                                               >  2:2292434/1‑93 (MQ=255)
               ttAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTtgtg                                                                              >  2:2553653/1‑93 (MQ=255)
                 aGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGtt                                                                            <  1:3993307/93‑1 (MQ=255)
                    tgtcACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAATAAACTCTGAAGCTGCTGACC‑‑GCGCGCGGCATCAACTTTGTCCTTGc                                                                         >  2:6760973/4‑93 (MQ=255)
                      gCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGa                                                                       <  2:5310494/93‑1 (MQ=255)
                       cACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGAc                                                                      >  1:7631314/1‑93 (MQ=255)
                        aCGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGAcc                                                                     >  1:4297375/1‑93 (MQ=255)
                         cGTTACTTTCTCCCGACCATACGTTATCTCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCg                                                                    <  1:4908404/93‑1 (MQ=255)
                         cGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCg                                                                    >  1:4422845/1‑93 (MQ=255)
                         cGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCg                                                                    <  2:2877792/93‑1 (MQ=255)
                          gTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCgg                                                                   >  2:5193617/1‑93 (MQ=255)
                           ttACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGgt                                                                  >  2:3980476/1‑93 (MQ=255)
                           ttACTTCCTCCCGACCAAACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGgt                                                                  <  2:5519271/93‑1 (MQ=255)
                            tACTTTCTCCCGACCATACGTTATCCTCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGgtc                                                                 >  2:7028716/1‑93 (MQ=255)
                              cTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGgtcgt                                                               <  2:7767711/93‑1 (MQ=255)
                               tttCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGgtcgtc                                                              >  2:772347/1‑93 (MQ=255)
                                 tctcCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcac                                                            >  2:328008/1‑93 (MQ=255)
                                 tctcCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcac                                                            >  2:1737714/1‑93 (MQ=255)
                                  ctcCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcacc                                                           >  2:2330889/1‑93 (MQ=255)
                                  ctcCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcacc                                                           >  1:1327617/1‑93 (MQ=255)
                                   tcCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcacca                                                          <  1:3412540/93‑1 (MQ=255)
                                   tcCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcacca                                                          >  1:2829260/1‑93 (MQ=255)
                                   tcCCGACCATACGTTATCACCTTACGCCAAAGAAACTCTGAAGCAACTAACC‑‑GCGCGCAACATCAAAAAAGTGTTTGCAACCGGTCGTcacca                                                          <  2:7243593/93‑1 (MQ=255)
                                   ccccGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcacca                                                          <  1:7073589/92‑1 (MQ=255)
                                     ccGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACg                                                        <  2:178779/93‑1 (MQ=255)
                                      cGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGt                                                       <  1:352636/93‑1 (MQ=255)
                                      cGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGt                                                       >  1:527696/1‑93 (MQ=255)
                                       gACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGtt                                                      <  1:4239098/93‑1 (MQ=255)
                                         ccATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGa                                                    <  2:6532634/93‑1 (MQ=255)
                                          cATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAt                                                   >  1:4585609/1‑93 (MQ=255)
                                            tACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGTGCGCGG‑‑CATTAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtgt                                                 <  1:6387172/93‑1 (MQ=255)
                                            tACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtgt                                                 >  2:1598729/1‑93 (MQ=255)
                                            tACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtgt                                                 <  1:17594/93‑1 (MQ=255)
                                            tACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtgg                                                 >  1:716918/1‑92 (MQ=255)
                                             aCGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtgtg                                                <  1:2025620/93‑1 (MQ=255)
                                             aCGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtgtg                                                <  1:6785052/93‑1 (MQ=255)
                                             aCGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtgtg                                                <  1:6482503/93‑1 (MQ=255)
                                              cGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTgg                                               <  1:6198000/93‑1 (MQ=255)
                                               gTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTggg                                              >  2:1148797/1‑93 (MQ=255)
                                               cttATCCCCCTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTggg                                              <  1:1860702/92‑1 (MQ=255)
                                                ttATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTgggg                                             >  1:2416596/1‑93 (MQ=255)
                                                  aTCCCCTTACGCCAAAGAAACTCTGAATCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCa                                           <  1:7668963/93‑1 (MQ=255)
                                                    ccccTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCaaa                                         <  2:7532522/93‑1 (MQ=255)
                                                      ctttACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAtt                                       <  2:7437710/91‑1 (MQ=255)
                                                      ccTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAtt                                       >  1:3047116/1‑93 (MQ=255)
                                                      ccTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAtt                                       >  2:358744/1‑93 (MQ=255)
                                                       cTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTc                                      >  1:1001671/1‑93 (MQ=255)
                                                       cTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTc                                      >  1:6675451/1‑93 (MQ=255)
                                                          aCGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTcgcg                                   >  2:4852883/1‑93 (MQ=255)
                                                           cGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGa                                  <  1:7348329/93‑1 (MQ=255)
                                                           cGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGa                                  <  1:6963509/93‑1 (MQ=255)
                                                            gCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGAt                                 >  2:2985876/1‑93 (MQ=255)
                                                             ccAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATa                                <  2:5445009/93‑1 (MQ=255)
                                                             ccAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATa                                >  1:7319086/1‑93 (MQ=255)
                                                              cAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATac                               >  2:3370487/1‑92 (MQ=255)
                                                               aaaGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAAt                              >  2:7151503/1‑93 (MQ=255)
                                                               aaaGAAACTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAAt                              >  1:452360/1‑93 (MQ=255)
                                                                   aaaCTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGGTGTGGGGCAAATTCGCGATCATCTgg                          >  1:469421/1‑93 (MQ=255)
                                                                    aaCTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCACATTCGCGATAATCTGga                         >  2:7934631/1‑93 (MQ=255)
                                                                    aaCTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGga                         <  1:4367854/93‑1 (MQ=255)
                                                                     aCTCTGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGgag                        >  1:815224/1‑93 (MQ=255)
                                                                       tctGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGAt                      >  1:5042827/1‑93 (MQ=255)
                                                                       tctGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGAt                      >  1:5979722/1‑93 (MQ=255)
                                                                        ctGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGAtt                     <  2:4592392/93‑1 (MQ=255)
                                                                         tGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTa                    <  1:5411367/93‑1 (MQ=255)
                                                                         tGAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTa                    <  1:7568998/93‑1 (MQ=255)
                                                                          gAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTaa                   >  1:5608898/1‑93 (MQ=255)
                                                                          gAAGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTaa                   <  2:6917084/93‑1 (MQ=255)
                                                                           aaGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAg                  >  1:3883251/1‑93 (MQ=255)
                                                                           aaGCTGCTCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAg                  >  2:4002189/1‑93 (MQ=255)
                                                                             gctgctCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTc                <  1:3298344/93‑1 (MQ=255)
                                                                               tgctCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCtt              >  2:598236/1‑93 (MQ=255)
                                                                                gctCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTa             <  1:5845904/93‑1 (MQ=255)
                                                                                gctCACC‑‑GCGCGCGGCATCAACTTCGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTa             <  2:5719786/93‑1 (MQ=255)
                                                                                 ctCACC‑‑GCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTAc            >  1:7451805/1‑93 (MQ=255)
                                                                                  cccgc‑‑gcgcgG‑‑CATCAACTTTGTGGTTGCGACCGGTCGGCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATg         <  2:6665115/89‑1 (MQ=255)
                                                                                  accgc‑‑gcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATg         <  2:2479113/89‑1 (MQ=255)
                                                                                  accgc‑‑gcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATg         >  2:4253665/5‑93 (MQ=255)
                                                                                  accgc‑‑gcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATg         >  2:4337479/5‑93 (MQ=255)
                                                                                  accgc‑‑gcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATg         <  2:6463337/89‑1 (MQ=255)
                                                                                   ccgc‑‑gcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGa        <  1:3190090/90‑1 (MQ=255)
                                                                                      c‑‑gcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTa     <  2:1842897/93‑1 (MQ=255)
                                                                                      c‑‑gcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTa     <  1:4108220/93‑1 (MQ=255)
                                                                                         gcgcgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCCCCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTAc    <  1:945156/93‑1 (MQ=255)
                                                                                         gcgcgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTAc    <  2:4742289/93‑1 (MQ=255)
                                                                                           gcgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTACCt  >  2:4473920/1‑93 (MQ=255)
                                                                                               ||                                                                                         
GTTGTTGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACC‑‑GCGCGC‑‑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 < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: