Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,333,953 A→G G61G (GGT→GGC murA ← UDP‑N‑glucosamine 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,333,9530AG99.6% 859.8 / NA 268G61G (GGT→GGCmurAUDP‑N‑glucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (152/115);  minor base C (0/1);  total (152/116)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 274 total aligned reads displayed.

ACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGC  >  NC_000913/3333863‑3334042
                                                                                          |                                                                                         
 cGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGccc                                                                                        <  1:2368586/93‑1 (MQ=255)
    aTGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGc                                                                                     <  1:5550938/93‑1 (MQ=255)
    aTGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGc                                                                                     >  2:1873127/1‑93 (MQ=255)
     tGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCt                                                                                    >  2:2443932/1‑93 (MQ=255)
        ttttAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGc                                                                                 >  1:4276759/1‑93 (MQ=255)
          ttAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCtt                                                                               >  1:3547951/1‑93 (MQ=255)
            aaCCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTag                                                                             >  1:5210972/1‑93 (MQ=255)
              ccAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTagca                                                                           >  2:249435/1‑93 (MQ=255)
                aGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTagcagc                                                                         >  2:5442586/1‑93 (MQ=255)
                 gATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCt                                                                        <  2:3790946/93‑1 (MQ=255)
                 gATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCt                                                                        >  1:2473825/1‑93 (MQ=255)
                   tCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTc                                                                      >  1:2807018/1‑93 (MQ=255)
                   tCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTc                                                                      >  1:5372193/1‑93 (MQ=255)
                   tCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTc                                                                      >  1:4036076/1‑93 (MQ=255)
                    cGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCa                                                                     >  1:1070775/1‑93 (MQ=255)
                     gTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCAt                                                                    >  1:3916579/1‑93 (MQ=255)
                     gTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCAt                                                                    <  1:1747881/93‑1 (MQ=255)
                     gTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCAt                                                                    >  2:1165636/1‑93 (MQ=255)
                      tAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCAtt                                                                   >  2:1699734/1‑93 (MQ=255)
                      tAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCAtt                                                                   <  2:1324946/93‑1 (MQ=255)
                      tAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCAtt                                                                   >  1:5657613/1‑93 (MQ=255)
                      tAAGGGGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCAtt                                                                   >  2:5412865/1‑93 (MQ=255)
                       aaGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTg                                                                  >  2:3231612/1‑93 (MQ=255)
                        aGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGa                                                                 >  1:2422717/1‑93 (MQ=255)
                         ggTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGAt                                                                <  1:2769250/93‑1 (MQ=255)
                             cgcAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTAt                                                            <  1:3876382/93‑1 (MQ=255)
                              gcAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATc                                                           >  1:6091855/1‑93 (MQ=255)
                               cAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCg                                                          >  2:5515200/1‑93 (MQ=255)
                                aGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGa                                                         >  1:5769154/1‑93 (MQ=255)
                                aGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGACGTATCGa                                                         <  2:5807939/93‑1 (MQ=255)
                                 gAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGAc                                                        >  2:4848038/1‑93 (MQ=255)
                                  aaTACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACg                                                       >  1:844038/1‑93 (MQ=255)
                                   aTACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGt                                                      >  1:6165523/1‑93 (MQ=255)
                                   aTACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGt                                                      >  2:477544/1‑93 (MQ=255)
                                     aCATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGtgt                                                    <  2:3296342/93‑3 (MQ=255)
                                      cATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCtt                                                   >  2:3823771/1‑93 (MQ=255)
                                      cATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATGGATGTATCGACGTCtt                                                   <  2:3392057/93‑1 (MQ=255)
                                       aTTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCttt                                                  <  1:5422618/93‑1 (MQ=255)
                                       aTTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCttt                                                  >  2:3029697/1‑93 (MQ=255)
                                        ttAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTc                                                 >  2:2961534/1‑93 (MQ=255)
                                        ttAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTc                                                 <  2:432498/93‑1 (MQ=255)
                                        ttAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGCGTTTc                                                 >  1:4527509/1‑93 (MQ=255)
                                         tAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCGGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCa                                                <  2:3518959/93‑1 (MQ=255)
                                            cGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGtt                                             <  1:379988/93‑1 (MQ=255)
                                            cGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGtt                                             <  1:5957271/93‑1 (MQ=255)
                                             gTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGttt                                            <  2:1097415/93‑1 (MQ=255)
                                              tCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTc                                           <  2:4388721/93‑1 (MQ=255)
                                              tCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTc                                           >  1:1887679/1‑93 (MQ=255)
                                                gcgGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCgg                                         <  2:686894/93‑1 (MQ=255)
                                                 cgGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCggg                                        <  1:4467782/93‑1 (MQ=255)
                                                  gggCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGa                                       <  2:4189335/93‑1 (MQ=255)
                                                  gggCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGa                                       >  2:625551/1‑93 (MQ=255)
                                                   gacATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCGTAGCAGCTTCATTTATGTATCGACGTCTTTCAGTTTCGGGAc                                      <  2:3730298/91‑1 (MQ=255)
                                                    gCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACg                                     >  2:4870236/1‑93 (MQ=255)
                                                      aTCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGtt                                   >  2:5318793/1‑93 (MQ=255)
                                                      aTCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGtt                                   >  1:6227297/1‑93 (MQ=255)
                                                       tCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTc                                  >  1:6054536/1‑93 (MQ=255)
                                                        cAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCt                                 <  1:5914386/93‑1 (MQ=255)
                                                        cAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCt                                 <  1:1001790/93‑1 (MQ=255)
                                                         aaTATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTg                                <  2:565370/93‑1 (MQ=255)
                                                         aaTATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTg                                >  2:2063312/1‑93 (MQ=255)
                                                         aaTATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTg                                <  2:258824/93‑1 (MQ=255)
                                                         aaTATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGAATCGACGTCTTTCAGTTTCGGGACGTTCTg                                >  1:919599/1‑93 (MQ=255)
                                                           tatGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGa                              >  1:3308878/1‑93 (MQ=255)
                                                            atGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGAt                             >  2:1204085/1‑93 (MQ=255)
                                                             tGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGAtc                            >  2:893122/1‑93 (MQ=255)
                                                             tGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTAGGGACGTTCTGGAtc                            >  1:699835/1‑93 (MQ=255)
                                                                acaGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATctct                         <  1:5086090/93‑1 (MQ=255)
                                                                acaGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATctct                         >  2:3881888/1‑93 (MQ=255)
                                                                 caGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTa                        >  1:2243608/1‑93 (MQ=255)
                                                                 caGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTa                        <  2:4423507/93‑1 (MQ=255)
                                                                  aGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTAc                       <  1:5307852/93‑1 (MQ=255)
                                                                     aCCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCgg                    <  2:5543174/93‑1 (MQ=255)
                                                                      ccATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGt                   <  1:1651638/93‑1 (MQ=255)
                                                                       cATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGtt                  >  2:6308396/1‑93 (MQ=255)
                                                                         ttACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCCTTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGTGTCTCTACCGGttct                >  1:1367989/1‑93 (MQ=255)
                                                                         ttACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTTTACCGGttct                <  2:4693633/93‑1 (MQ=255)
                                                                         ttACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGttct                >  2:2573581/1‑93 (MQ=255)
                                                                           aCGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGttcttc              <  1:2993327/93‑1 (MQ=255)
                                                                           aCGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGttcttc              <  2:3312511/93‑1 (MQ=255)
                                                                            cGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTcc             <  2:4042861/93‑1 (MQ=255)
                                                                            cGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTcc             <  2:2239679/93‑1 (MQ=255)
                                                                             gTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCg            >  1:5001029/1‑93 (MQ=255)
                                                                              ttCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGc           <  2:1992929/93‑1 (MQ=255)
                                                                               tCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGcc          >  1:5728635/1‑93 (MQ=255)
                                                                               tCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGcc          >  2:817451/1‑93 (MQ=255)
                                                                               tCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGcc          >  2:2610377/1‑93 (MQ=255)
                                                                                 tACTTTCGCCCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCag        <  1:890750/93‑1 (MQ=255)
                                                                                  aCTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCagt       <  1:589943/93‑1 (MQ=255)
                                                                                  aCTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCagt       >  2:6197652/1‑93 (MQ=255)
                                                                                   cTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCagta      <  1:3046563/93‑1 (MQ=255)
                                                                                   cTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCagta      >  2:1439716/1‑93 (MQ=255)
                                                                                   cTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCagta      <  1:4804748/93‑1 (MQ=255)
                                                                                   cTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCagta      <  2:3166900/93‑1 (MQ=255)
                                                                                    tttCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCagtag     >  1:19874/1‑93 (MQ=255)
                                                                                    ttccgcgCCCAGCTGGCTTAGCAGCCTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCagtag     <  1:4183820/90‑1 (MQ=255)
                                                                                     ttCGCGCCCATCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCagtagt    <  1:3251783/93‑1 (MQ=255)
                                                                                     ttCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCagtagt    <  1:3108196/93‑1 (MQ=255)
                                                                                      tCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTg   >  2:5281705/1‑93 (MQ=255)
                                                                                       cgcgCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGc  <  1:970253/93‑1 (MQ=255)
                                                                                          |                                                                                         
ACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGC  >  NC_000913/3333863‑3334042

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: