Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,949,410 2 bp→TC coding (588‑589/600 nt) yecD → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,949,4100CT75.3% 1935.5 / inf 936L196L (CTC→CTTyecDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (129/102);  new base T (345/360);  total (474/462)
Fisher's exact test for biased strand distribution p-value = 6.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01
*NC_0009131,949,4110AC75.6% 2359.3 / inf 931N197H (AAC→CAC) yecDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (125/102);  new base C (344/360);  total (469/462)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

Only 100 of 965 total aligned reads displayed.

GCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACT  >  NC_000913/1949318‑1949503
                                                                                            ||                                                                                            
gCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCAGATCTACCCGCGCATCGCACGTGTGCGTAGCGGCGAAGAGGTCCTc                                                                                               >  2:2637820/1‑93 (MQ=255)
 ccTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCa                                                                                              >  2:4401728/1‑93 (MQ=255)
   tGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGAGGAAGAGATCCTCAAc                                                                                            >  1:2466212/1‑93 (MQ=255)
    gTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAcg                                                                                           >  2:1409311/1‑93 (MQ=255)
     tAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCAcgc                                                                                          >  1:3314500/1‑93 (MQ=255)
       gTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGt                                                                                        >  1:5309802/1‑93 (MQ=255)
       gTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGt                                                                                        >  2:1358173/1‑93 (MQ=255)
        tGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGtt                                                                                       <  2:4681348/93‑1 (MQ=255)
        tGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGtt                                                                                       >  1:6009064/1‑93 (MQ=255)
        tGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGtt                                                                                       >  1:988814/1‑93 (MQ=255)
        tGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGtt                                                                                       <  2:3792162/93‑1 (MQ=255)
            gccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGa                                                                                   <  1:3731880/93‑1 (MQ=255)
            gccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGa                                                                                   <  1:3403970/93‑1 (MQ=255)
            gccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGa                                                                                   <  2:192540/93‑1 (MQ=255)
             ccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGAt                                                                                  >  1:2694504/1‑93 (MQ=255)
             ccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAt                                                                                  <  1:1646920/93‑1 (MQ=255)
                gcgcCGAGCAACACAATAACAGCATTAATCATAGCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTa                                                                               <  2:5937633/93‑1 (MQ=255)
                   ccGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACAt                                                                            >  1:363980/1‑93 (MQ=255)
                     gAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCg                                                                          >  1:4523542/1‑93 (MQ=255)
                     gAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCg                                                                          <  1:506008/93‑1 (MQ=255)
                     gAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCg                                                                          >  2:3295748/1‑93 (MQ=255)
                       gcaacaCAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGCGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGt                                                                        <  2:78238/93‑1 (MQ=255)
                        caacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTc                                                                       <  2:3617424/90‑1 (MQ=255)
                        caacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTc                                                                       <  1:1001413/90‑1 (MQ=255)
                        caacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTc                                                                       <  2:2080578/90‑1 (MQ=255)
                         aacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCt                                                                      >  1:2526304/3‑93 (MQ=255)
                         aacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCt                                                                      <  2:1101445/91‑1 (MQ=255)
                         aacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCt                                                                      <  2:64237/91‑1 (MQ=255)
                          acacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCtg                                                                     >  1:4106912/2‑92 (MQ=255)
                           cacaATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTAc                                                                    <  1:4976843/93‑1 (MQ=255)
                            acaATAACAGCATTACTCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGCTTTACATCGGTCTAcc                                                                   <  1:5481195/93‑1 (MQ=255)
                              aaTAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGc                                                                 >  1:562527/1‑93 (MQ=255)
                               aTAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCa                                                                >  1:2276012/1‑93 (MQ=255)
                                tAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCaa                                                               >  2:2831670/1‑93 (MQ=255)
                                 aaCAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAAt                                                              >  1:5200052/1‑93 (MQ=255)
                                   cAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGCAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATgg                                                            <  1:4184678/93‑1 (MQ=255)
                                       aTCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGc                                                        >  2:2161611/1‑93 (MQ=255)
                                        tCAGCCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCa                                                       <  1:5837806/93‑1 (MQ=255)
                                         taatcATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAt                                                      >  2:73990/5‑93 (MQ=255)
                                         taatcATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAt                                                      <  1:2870243/89‑1 (MQ=255)
                                          aatcATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATc                                                     <  2:5404844/90‑1 (MQ=255)
                                           aCCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATcc                                                    >  1:4400584/1‑93 (MQ=255)
                                           aCCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATcc                                                    <  1:4267920/93‑1 (MQ=255)
                                            ccATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCt                                                   >  2:166393/1‑93 (MQ=255)
                                             cATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTa                                                  <  2:5181746/93‑1 (MQ=255)
                                             cATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTa                                                  >  1:1089758/1‑93 (MQ=255)
                                             cATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTa                                                  >  1:3191016/1‑93 (MQ=255)
                                               tatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaaa                                                >  2:2529704/1‑93 (MQ=255)
                                               tatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaaa                                                <  2:4249160/93‑1 (MQ=255)
                                               tatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaaa                                                <  1:2354942/93‑1 (MQ=255)
                                               tatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaaa                                                >  2:1277546/1‑93 (MQ=255)
                                                   tACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATggg                                            >  2:229231/1‑93 (MQ=255)
                                                   tACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATggg                                            <  1:1860266/93‑1 (MQ=255)
                                                    aCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATgggg                                           >  1:1201859/1‑92 (MQ=255)
                                                     cccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTg                                          >  1:3545024/1‑93 (MQ=255)
                                                     cccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTg                                          <  2:3489987/93‑1 (MQ=255)
                                                      ccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGc                                         <  2:5279783/93‑1 (MQ=255)
                                                      ccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGc                                         >  2:3960977/1‑93 (MQ=255)
                                                      ccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGc                                         >  1:3492415/1‑93 (MQ=255)
                                                       cgcgCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCg                                        >  2:5622812/1‑93 (MQ=255)
                                                       cgcgCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCg                                        >  2:3381483/1‑93 (MQ=255)
                                                       cgcgCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCg                                        <  2:3190331/93‑1 (MQ=255)
                                                        gcgcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCgt                                       >  1:3114568/1‑92 (MQ=255)
                                                         cgcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGt                                      <  2:895219/93‑1 (MQ=255)
                                                         cgcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGt                                      <  2:1542170/93‑1 (MQ=255)
                                                         cgcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGt                                      >  1:6092082/1‑93 (MQ=255)
                                                          gcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGtt                                     >  2:3092592/1‑93 (MQ=255)
                                                          gcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGtt                                     >  1:1803562/1‑93 (MQ=255)
                                                          gcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGtt                                     >  1:3824324/1‑93 (MQ=255)
                                                          gcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGtt                                     >  2:175179/1‑93 (MQ=255)
                                                            aTCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTgg                                   >  1:1442945/1‑93 (MQ=255)
                                                             tCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg                                  <  1:1563947/93‑1 (MQ=255)
                                                             tCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg                                  <  2:5670398/93‑1 (MQ=255)
                                                                 ccGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATc                              >  1:2798357/1‑93 (MQ=255)
                                                                   gtgtGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCAc                            >  2:4903417/1‑93 (MQ=255)
                                                                    tgtgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCAcc                           >  1:1948767/1‑93 (MQ=255)
                                                                    tgtgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCAcc                           >  2:5014586/1‑93 (MQ=255)
                                                                    ggtgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCAcc                           <  1:1356310/92‑1 (MQ=255)
                                                                      tgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAg                         <  2:6071380/93‑1 (MQ=255)
                                                                         gTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCt                      >  2:3700711/1‑93 (MQ=255)
                                                                         gTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCt                      >  1:206038/1‑93 (MQ=255)
                                                                              gTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAag                 >  2:4541065/1‑93 (MQ=255)
                                                                               tGGAAGAGATCCTTCACGCTTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAaga                <  1:2993388/93‑1 (MQ=255)
                                                                               tGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAaga                <  2:4807115/93‑1 (MQ=255)
                                                                               tGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTCAATGGGTGCGGTTGGGGATCACCAGCCTTGAaga                >  2:2977387/1‑93 (MQ=255)
                                                                               tGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAaga                <  1:4834548/93‑1 (MQ=255)
                                                                                ggAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAagag               >  2:5789025/1‑93 (MQ=255)
                                                                                ggAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAagag               <  2:4148018/93‑1 (MQ=255)
                                                                                  aaGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTa             <  1:4696007/93‑1 (MQ=255)
                                                                                  aaGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTa             <  2:549488/93‑1 (MQ=255)
                                                                                  aaGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTa             <  1:690361/93‑1 (MQ=255)
                                                                                  aaGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAACCTTGAAGAGTa             <  1:5396749/93‑1 (MQ=255)
                                                                                  aaGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTa             <  1:2165895/93‑1 (MQ=255)
                                                                                   agatATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTAt            <  2:1835789/93‑1 (MQ=255)
                                                                                   agagATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTAt            <  1:911419/93‑1 (MQ=255)
                                                                                    gagaTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATg           <  1:5721804/93‑1 (MQ=255)
                                                                                    gagaTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTCAAGAGTATg           >  2:1196560/1‑93 (MQ=255)
                                                                                     agaTCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGc          >  1:3926228/1‑93 (MQ=255)
                                                                                      gaTCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGcc         <  2:5110330/93‑1 (MQ=255)
                                                                                          cTCCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGCCCGcc     <  1:2632954/93‑1 (MQ=255)
                                                                                            ||                                                                                            
GCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACT  >  NC_000913/1949318‑1949503

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