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.9% 1620.3 / 107.1 506intergenic (+109/‑11)pldB/b4555lysophospholipase L(2)/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (21/15);  new base G (245/225);  total (266/240)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-01
*NC_0009134,007,9072.G92.9% 1620.0 / 107.1 506intergenic (+109/‑11)pldB/b4555lysophospholipase L(2)/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (21/15);  new base G (245/225);  total (266/240)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-02

Only 100 of 544 total aligned reads displayed.

GTTGTTGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGC‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTACCT  >  NC_000913/4007815‑4007996
                                                                                             ||                                                                                         
gttgttGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCgcgcgc                                                                                             <  1:2017977/93‑1 (MQ=255)
 ttgttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCgcgcgcg                                                                                            <  2:6970413/93‑2 (MQ=255)
  tgttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCgcgcgcgg                                                                                           >  1:807806/1‑91 (MQ=255)
   gttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCgcgcgcggc                                                                                          >  2:7288327/1‑90 (MQ=255)
           tGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACtt                                                                                  >  1:6747736/1‑93 (MQ=255)
            gATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACttt                                                                                 >  1:2813422/1‑93 (MQ=255)
             aTTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTtg                                                                                <  2:4609783/93‑1 (MQ=255)
              tttAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTtgt                                                                               >  2:4711996/1‑93 (MQ=255)
                tAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTgtgt                                                                             >  2:5039626/1‑93 (MQ=255)
                 aGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGtt                                                                            <  1:6498693/93‑1 (MQ=255)
                  gATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGttt                                                                           >  2:4854651/1‑93 (MQ=255)
                  gATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGttt                                                                           <  1:396926/93‑1 (MQ=255)
                  gATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGttt                                                                           <  1:5793612/93‑1 (MQ=255)
                    tGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGc                                                                         <  2:5241722/93‑1 (MQ=255)
                    tGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGc                                                                         >  2:4527880/1‑93 (MQ=255)
                     ggCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCg                                                                        <  2:4429868/93‑1 (MQ=255)
                      gCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGa                                                                       <  1:7157888/93‑1 (MQ=255)
                       cACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGAc                                                                      <  2:3406068/93‑1 (MQ=255)
                       cACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGAc                                                                      <  2:7043524/93‑1 (MQ=255)
                        aCGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGAcc                                                                     <  2:163631/93‑1 (MQ=255)
                         cGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCg                                                                    <  1:519565/93‑1 (MQ=255)
                         cGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCg                                                                    >  2:1313634/1‑93 (MQ=255)
                           ttACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGgt                                                                  <  2:6583604/93‑1 (MQ=255)
                           ttACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGgt                                                                  <  2:956609/93‑1 (MQ=255)
                            tACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGgtc                                                                 >  1:6960091/1‑93 (MQ=255)
                             aCTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGgtcg                                                                <  1:3386232/93‑1 (MQ=255)
                             aCTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGgtcg                                                                <  2:1502640/93‑1 (MQ=255)
                              cTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGgtcgt                                                               <  2:2468824/93‑1 (MQ=255)
                              cTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGgtcgt                                                               <  2:1209532/93‑1 (MQ=255)
                              cTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGgtcgt                                                               >  1:1861644/1‑93 (MQ=255)
                               tttCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGgtcgtc                                                              >  2:6842907/1‑93 (MQ=255)
                               tttCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGgtcgtc                                                              >  1:448767/1‑93 (MQ=255)
                                ttCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTca                                                             <  2:669204/93‑1 (MQ=255)
                                ttCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACGCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTca                                                             >  1:6217778/1‑93 (MQ=255)
                                 tctcCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcac                                                            <  1:2566057/93‑1 (MQ=255)
                                  ctcCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcacc                                                           >  2:4206429/1‑93 (MQ=255)
                                   tcCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcacca                                                          >  1:2685147/1‑93 (MQ=255)
                                    cccGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcaccac                                                         >  2:4021280/1‑93 (MQ=255)
                                    cccGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcaccac                                                         <  1:6650267/93‑1 (MQ=255)
                                      cGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGt                                                       <  1:3720825/93‑1 (MQ=255)
                                      cGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGt                                                       <  2:3812354/93‑1 (MQ=255)
                                         ccATACGTCATCCCCTTACTCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCACCTTTGTGTTTGCGACCGCTGTTCACCACGTTgt                                                    >  1:7041890/1‑92 (MQ=255)
                                           aTACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtg                                                  >  2:3879063/1‑93 (MQ=255)
                                           aTACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtg                                                  >  1:1058281/1‑93 (MQ=255)
                                            tACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtgt                                                 >  1:5721527/1‑93 (MQ=255)
                                             aCGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtgtg                                                <  2:2837977/93‑1 (MQ=255)
                                               gTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTggg                                              >  2:6473092/1‑93 (MQ=255)
                                               gTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTggg                                              <  2:2179138/93‑1 (MQ=255)
                                                  aTCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCa                                           >  2:3682623/1‑93 (MQ=255)
                                                  aTCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCa                                           <  1:951554/93‑1 (MQ=255)
                                                    ccccTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCaaa                                         <  1:5108248/93‑1 (MQ=255)
                                                    ccccTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCaaa                                         <  2:2083087/93‑1 (MQ=255)
                                                    ccccTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCaaa                                         <  1:5935591/93‑1 (MQ=255)
                                                     cccTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAt                                        >  1:2181967/1‑93 (MQ=255)
                                                      ccTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAtt                                       >  2:2690715/1‑93 (MQ=255)
                                                       cTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTc                                      >  1:6155180/1‑93 (MQ=255)
                                                        ttACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTcg                                     >  1:1247725/1‑93 (MQ=255)
                                                        ttACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTcg                                     >  2:1917817/1‑93 (MQ=255)
                                                        ttACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTcg                                     >  2:75463/1‑93 (MQ=255)
                                                           cGCCAAAGAGACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGa                                  <  1:1007851/93‑1 (MQ=255)
                                                           cGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGa                                  <  1:1368981/93‑1 (MQ=255)
                                                            gCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGAt                                 >  1:3081487/1‑93 (MQ=255)
                                                            gCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGAt                                 >  2:5426830/1‑93 (MQ=255)
                                                            gCCAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGAt                                 >  1:2333553/1‑93 (MQ=255)
                                                             ccAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATa                                >  2:2344294/1‑93 (MQ=255)
                                                             ccAAAGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATa                                >  1:2469024/1‑93 (MQ=255)
                                                               aaaGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAAt                              >  2:314821/1‑93 (MQ=255)
                                                               aaaGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAAt                              <  1:155648/93‑1 (MQ=255)
                                                                 aGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCt                            >  1:1948257/1‑93 (MQ=255)
                                                                 aGAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCt                            <  1:2949058/93‑1 (MQ=255)
                                                                  gAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTg                           <  2:5854693/93‑1 (MQ=255)
                                                                  gAAACTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTg                           >  1:1098324/1‑93 (MQ=255)
                                                                   aaaCTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTgg                          >  1:5637021/1‑93 (MQ=255)
                                                                    aaCTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGga                         >  1:6573949/1‑93 (MQ=255)
                                                                    aaCTCTGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGga                         <  2:3284494/93‑1 (MQ=255)
                                                                      ctctGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGgaga                       >  2:6931058/1‑93 (MQ=255)
                                                                       tctGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGAt                      <  1:736040/93‑1 (MQ=255)
                                                                       tctGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGAt                      >  2:6386156/1‑93 (MQ=255)
                                                                        ctGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGAtt                     >  1:6350828/1‑93 (MQ=255)
                                                                        ctGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGAtt                     >  2:3537426/1‑93 (MQ=255)
                                                                         tGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTa                    >  1:1675250/1‑93 (MQ=255)
                                                                         tGAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTa                    <  1:3213359/93‑1 (MQ=255)
                                                                          gAAGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTaa                   >  1:3582579/1‑93 (MQ=255)
                                                                           aaGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCCCCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAg                  <  1:3472769/93‑1 (MQ=255)
                                                                            aGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGt                 >  1:4776290/1‑93 (MQ=255)
                                                                            aGCTGCTCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGt                 <  2:5925676/93‑1 (MQ=255)
                                                                              ctgctCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCt               >  2:2917284/1‑93 (MQ=255)
                                                                               tgctCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCtt              >  2:4972745/1‑93 (MQ=255)
                                                                                gctCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTa             <  2:5693211/93‑1 (MQ=255)
                                                                                gctCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTa             <  2:5626111/93‑1 (MQ=255)
                                                                                gctCACCGCGCGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTa             <  2:4117181/93‑1 (MQ=255)
                                                                                   ccgcgcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGa        <  1:4210282/90‑1 (MQ=255)
                                                                                    cgcgcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGAt       >  1:4969283/3‑93 (MQ=255)
                                                                                    cgcgcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGAt       <  2:1837918/91‑1 (MQ=255)
                                                                                     gcgcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGAtt      >  2:6162928/2‑93 (MQ=255)
                                                                                     gcgcgcgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGAtt      >  1:5364982/2‑93 (MQ=255)
                                                                                       gcgcgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTAc    <  2:3045329/93‑1 (MQ=255)
                                                                                        cgcgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTAcc   <  2:2290733/93‑1 (MQ=255)
                                                                                         gcgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCCCCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTACCt  <  1:4562380/93‑1 (MQ=255)
                                                                                         gcgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTACCt  >  1:7457942/1‑93 (MQ=255)
                                                                                             ||                                                                                         
GTTGTTGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCGCGC‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTACCT  >  NC_000913/4007815‑4007996

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