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.2% 733.2 / 49.2 218intergenic (+109/‑11)pldB/b4555lysophospholipase L(2)/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (7/10);  new base G (98/103);  total (105/113)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
*NC_0009134,007,9072.G92.2% 720.3 / 49.2 218intergenic (+109/‑11)pldB/b4555lysophospholipase L(2)/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (7/10);  new base G (97/103);  total (105/113)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 241 total aligned reads displayed.

GTTGTTGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GC‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTACC  >  NC_000913/4007815‑4007995
                                                                                                     ||                                                                                        
gttgttGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCg‑‑c‑‑g‑‑c‑‑gc                                                                                            >  2:2852694/1‑93 (MQ=255)
 ttgttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCg‑‑c‑‑g‑‑c‑‑gcg                                                                                           >  2:5355088/1‑92 (MQ=255)
 ttgttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCg‑‑c‑‑g‑‑c‑‑gcg                                                                                           <  2:5382382/93‑2 (MQ=255)
    ttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCg‑‑c‑‑g‑‑c‑‑gcggca                                                                                        <  1:2409349/93‑5 (MQ=255)
    ttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCg‑‑c‑‑g‑‑c‑‑gcggca                                                                                        >  2:5628778/1‑89 (MQ=255)
          cTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACt                                                                                  >  1:5650964/1‑93 (MQ=255)
          cTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACt                                                                                  >  1:2729426/1‑93 (MQ=255)
           tGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACtt                                                                                 >  2:4326864/1‑93 (MQ=255)
           tGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACtt                                                                                 >  1:4770032/1‑93 (MQ=255)
            gATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACttt                                                                                <  2:1008892/93‑1 (MQ=255)
              tttAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTtgt                                                                              >  1:2507896/1‑93 (MQ=255)
              tttAGATGGCACGTTAATTTCTCCCGAGCATACGTTATCACCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTtgg                                                                              >  1:1517712/1‑92 (MQ=255)
                tAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTgtgt                                                                            <  2:2680729/93‑1 (MQ=255)
                tAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTgtgt                                                                            >  2:4567741/1‑93 (MQ=255)
                tAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTgtgt                                                                            <  2:3092772/93‑1 (MQ=255)
                 aGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGtt                                                                           >  2:705389/1‑93 (MQ=255)
                  gATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGttt                                                                          >  1:1930036/1‑93 (MQ=255)
                   aTGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTg                                                                         <  1:687378/93‑1 (MQ=255)
                   aTGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTg                                                                         <  1:5307627/93‑1 (MQ=255)
                    tGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGc                                                                        >  2:5198404/1‑93 (MQ=255)
                      gCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGa                                                                      <  1:5043115/93‑1 (MQ=255)
                       cACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGAc                                                                     <  2:569028/93‑1 (MQ=255)
                         cGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCg                                                                   >  1:5443982/1‑93 (MQ=255)
                         cGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCg                                                                   <  2:545551/93‑1 (MQ=255)
                          gTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCgg                                                                  >  1:345489/1‑93 (MQ=255)
                          gTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCgg                                                                  >  1:1433611/1‑93 (MQ=255)
                           ttACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGgt                                                                 <  1:6062132/93‑1 (MQ=255)
                            tACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGgtc                                                                >  1:5031950/1‑93 (MQ=255)
                             aCTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGgtcg                                                               >  1:1013704/1‑93 (MQ=255)
                              cTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGgtcgt                                                              <  2:2237971/93‑1 (MQ=255)
                              cTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGgtcgt                                                              <  1:3768692/93‑1 (MQ=255)
                               tttCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGgtcgtc                                                             >  2:3369756/1‑93 (MQ=255)
                               tttCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGACATCAACTTTGTGTTTGCGACCGgtcgtc                                                             >  1:5392803/1‑93 (MQ=255)
                                ttCTCCCGACCATACGTTTTCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTca                                                            >  1:2887989/1‑93 (MQ=255)
                                 tctcCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTcac                                                           <  2:2296746/93‑1 (MQ=255)
                                   tcCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCGCCC‑‑G‑‑C‑‑GG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTcaccg                                                         >  2:1491889/1‑92 (MQ=255)
                                    cccGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTcaccac                                                        <  2:2111790/93‑1 (MQ=255)
                                     ccGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACg                                                       <  2:4767832/93‑1 (MQ=255)
                                      cGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTAGTGTTTGCGACCGGGCGTCACCACGt                                                      >  1:5489590/1‑93 (MQ=255)
                                       gACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGtt                                                     <  1:4930589/93‑1 (MQ=255)
                                       gACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGtt                                                     >  1:1204093/1‑93 (MQ=255)
                                        aCCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTg                                                    >  2:1653441/1‑93 (MQ=255)
                                         ccATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGa                                                   <  1:4343567/93‑1 (MQ=255)
                                         ccATACGTTATCCACTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑CGGGG‑‑CATCAACTTTGTGTCTGCGACCGGTCGTCACCACGTTGa                                                   <  1:460414/93‑1 (MQ=255)
                                              cGTTATCCCATTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTgg                                              <  2:4175267/93‑1 (MQ=255)
                                               gTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTggg                                             <  1:2992981/93‑1 (MQ=255)
                                                ttATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTgggg                                            >  2:2582432/1‑93 (MQ=255)
                                                ttATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTgggg                                            >  2:3864478/1‑93 (MQ=255)
                                                 tATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGc                                           >  1:1161352/1‑93 (MQ=255)
                                                  aTCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCa                                          >  2:4783699/1‑93 (MQ=255)
                                                   tCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCaa                                         >  2:409587/1‑93 (MQ=255)
                                                    ccccTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGaaaa                                        >  1:3222974/1‑93 (MQ=255)
                                                    ccccTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCaaa                                        <  2:5557182/93‑1 (MQ=255)
                                                    ccccTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCaaa                                        <  2:4941312/93‑1 (MQ=255)
                                                    ccccTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCaaa                                        <  2:836463/93‑1 (MQ=255)
                                                     cccTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAt                                       <  1:2100006/93‑1 (MQ=255)
                                                       cTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTc                                     >  1:3671966/1‑93 (MQ=255)
                                                         tACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTcgc                                   <  1:4652842/93‑1 (MQ=255)
                                                         tACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTcgc                                   <  1:5968217/93‑1 (MQ=255)
                                                           cGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGa                                 <  2:160920/93‑1 (MQ=255)
                                                            gCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGAt                                <  2:3740822/93‑1 (MQ=255)
                                                             ccAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATa                               >  1:3477861/1‑93 (MQ=255)
                                                              cAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATaa                              <  2:1078968/93‑1 (MQ=255)
                                                              cAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATaa                              >  1:788575/1‑93 (MQ=255)
                                                              cAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATaa                              <  2:382801/93‑1 (MQ=255)
                                                               aaaGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAAt                             >  1:583053/1‑93 (MQ=255)
                                                                aaGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATc                            >  2:1706846/1‑93 (MQ=255)
                                                                 aGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCt                           <  2:3933317/93‑1 (MQ=255)
                                                                  gAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTg                          >  1:5229002/1‑93 (MQ=255)
                                                                  gAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTg                          <  1:5907996/93‑1 (MQ=255)
                                                                   aaaCTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTgg                         <  1:1549342/93‑1 (MQ=255)
                                                                    aaCTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGga                        >  2:3284101/1‑93 (MQ=255)
                                                                    aaCTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGga                        <  2:2044645/93‑1 (MQ=255)
                                                                    aaCTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGga                        <  1:4536281/93‑1 (MQ=255)
                                                                    aaCTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGga                        <  2:105218/93‑1 (MQ=255)
                                                                    aaCTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGga                        >  1:5347294/1‑93 (MQ=255)
                                                                    aaCTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGga                        >  2:3302979/1‑93 (MQ=255)
                                                                       tctGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGAt                     >  1:5081484/1‑93 (MQ=255)
                                                                         tGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTa                   <  1:4083431/93‑1 (MQ=255)
                                                                          gAAGCTGCTCACCG‑‑C‑‑GCAC‑‑GG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTaa                  <  1:5475936/93‑1 (MQ=255)
                                                                          gAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTaa                  <  2:3234756/93‑1 (MQ=255)
                                                                           aaGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAg                 <  1:2554444/93‑1 (MQ=255)
                                                                            aGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGCTAATCTGGAGATTAAGt                >  2:1310580/1‑93 (MQ=255)
                                                                            aGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGt                <  2:4639733/93‑1 (MQ=255)
                                                                            aGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAATCGCGATAATCTGGAGATTAAGt                >  1:4301358/1‑93 (MQ=255)
                                                                            aGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTGTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGt                <  2:5950835/93‑1 (MQ=255)
                                                                             gctgctCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTc               >  2:769262/1‑93 (MQ=255)
                                                                             gctgctCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTc               >  2:5270429/1‑93 (MQ=255)
                                                                             gctgctCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTc               <  2:143704/93‑1 (MQ=255)
                                                                              ctgctCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGCTTAAGTCt              >  2:1563447/1‑93 (MQ=255)
                                                                              ctgctCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCt              <  1:1001325/93‑1 (MQ=255)
                                                                               tgctCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCtt             >  2:4821017/1‑93 (MQ=255)
                                                                                gctCACCG‑‑CGGG‑‑C‑‑GG‑‑CATCAAATTTGTGTTTGCGACCGGTAGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTa            >  1:1086467/1‑93 (MQ=255)
                                                                                gctCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTa            <  2:5422841/93‑1 (MQ=255)
                                                                                gctCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTa            <  1:3326616/93‑1 (MQ=255)
                                                                                gctCACCG‑‑C‑‑G‑‑C‑‑GCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTa            <  2:2599361/93‑1 (MQ=255)
                                                                                  accgcg‑‑c‑‑g‑‑c‑‑gG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATg        <  1:3083372/89‑1 (MQ=255)
                                                                                   ccgcg‑‑c‑‑g‑‑c‑‑gG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGa       <  2:697684/90‑1 (MQ=255)
                                                                                      cg‑‑c‑‑g‑‑c‑‑gG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTa    <  1:2816241/93‑1 (MQ=255)
                                                                                       g‑‑c‑‑g‑‑c‑‑gg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGTTTAc   >  2:1880644/1‑93 (MQ=255)
                                                                                                     ||                                                                                        
GTTGTTGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑C‑‑G‑‑C‑‑GC‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTACC  >  NC_000913/4007815‑4007995

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