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,9530AG98.9% 574.8 / ‑3.2 182G61G (GGT→GGC) murAUDP‑N‑glucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (94/86);  total (96/86)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 189 total aligned reads displayed.

ACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGC  >  NC_000913/3333863‑3334042
                                                                                          |                                                                                         
aCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGcc                                                                                         >  1:4301018/1‑93 (MQ=255)
aCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGcc                                                                                         <  2:3968179/93‑1 (MQ=255)
 cGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGccc                                                                                        >  1:2303412/1‑93 (MQ=255)
   cATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAg                                                                                      >  2:2176173/1‑93 (MQ=255)
    aTGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGc                                                                                     >  1:4379901/1‑93 (MQ=255)
     tGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCt                                                                                    >  1:1633705/1‑93 (MQ=255)
     tGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCt                                                                                    >  2:134868/1‑93 (MQ=255)
     tGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCt                                                                                    <  1:1042880/93‑1 (MQ=255)
      ggTTTTAACCAGATCGTAAGGTGCTCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTCCGTTCTACTTTCGCGCCCAGCTg                                                                                   <  2:725306/93‑1 (MQ=255)
      ggTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTg                                                                                   <  1:602997/93‑1 (MQ=255)
      ggTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTg                                                                                   >  1:1135905/1‑93 (MQ=255)
       gTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTgg                                                                                  <  1:1560064/93‑1 (MQ=255)
         tttAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCt                                                                                >  2:4552103/1‑93 (MQ=255)
          ttAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCtt                                                                               >  1:1176345/1‑93 (MQ=255)
          ttAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCtt                                                                               <  1:3492638/93‑1 (MQ=255)
           tAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTa                                                                              >  1:2220734/1‑93 (MQ=255)
            aaCCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTag                                                                             <  2:48754/93‑1 (MQ=255)
             aCCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTagc                                                                            >  2:1504148/1‑93 (MQ=255)
              ccAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTagca                                                                           >  1:1880652/1‑93 (MQ=255)
              ccAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTagca                                                                           <  2:1246396/93‑1 (MQ=255)
               cAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTagcag                                                                          >  1:3681364/1‑93 (MQ=255)
                 gATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCt                                                                        >  1:2122277/1‑93 (MQ=255)
                  aTCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCtt                                                                       <  2:2800912/93‑1 (MQ=255)
                  aTCGTAAGCTGCGCAGAATTCATTAACGTCGCGTGCATGAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCtt                                                                       <  1:2521772/93‑1 (MQ=255)
                    cGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCa                                                                     >  1:653215/1‑93 (MQ=255)
                     gTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCAt                                                                    <  2:1931444/93‑1 (MQ=255)
                     gTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCAt                                                                    <  1:3925419/93‑1 (MQ=255)
                     gTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCAt                                                                    >  1:3819455/1‑93 (MQ=255)
                     gTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCAt                                                                    >  2:1733006/1‑93 (MQ=255)
                      tAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCAtt                                                                   <  2:428512/93‑1 (MQ=255)
                      tAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCAtt                                                                   <  2:1424048/93‑1 (MQ=255)
                        aGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGa                                                                 <  2:2357514/93‑1 (MQ=255)
                         ggTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGAt                                                                >  1:2314276/1‑93 (MQ=255)
                          gTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATg                                                               <  1:1986064/93‑1 (MQ=255)
                           tGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGt                                                              <  2:447264/93‑1 (MQ=255)
                            gcgcAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTa                                                             <  1:1003797/93‑1 (MQ=255)
                            gcgcAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTa                                                             <  1:4525737/93‑1 (MQ=255)
                             cgcAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTAt                                                            <  2:3541284/93‑1 (MQ=255)
                              gcAGAATACATTAACGTCGGGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATc                                                           <  1:920740/93‑1 (MQ=255)
                              gcAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATc                                                           >  1:1918421/1‑93 (MQ=255)
                              gcAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATc                                                           <  2:2726480/93‑1 (MQ=255)
                                aGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGa                                                         <  1:2867957/93‑1 (MQ=255)
                                aGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGa                                                         <  1:2141531/93‑1 (MQ=255)
                                aGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGa                                                         >  2:4162043/1‑93 (MQ=255)
                                  aaTACATTAACGTCGCGCGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCTAGCTGGCTTAGCAGCTTCATTGATGTATCGACg                                                       <  1:1081467/93‑1 (MQ=255)
                                   aTACATTAACGTCGCGGGCAGCAATATGCACAGAACCATTATGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGt                                                      >  2:3037743/1‑93 (MQ=255)
                                    tACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTc                                                     >  2:785100/1‑93 (MQ=255)
                                    tACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTc                                                     >  1:201233/1‑93 (MQ=255)
                                     aCATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCt                                                    <  1:3528630/93‑1 (MQ=255)
                                     aCATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCt                                                    <  1:2969965/93‑1 (MQ=255)
                                     aCATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCt                                                    >  2:3356517/1‑93 (MQ=255)
                                     aCATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCt                                                    <  1:4027490/93‑1 (MQ=255)
                                     aCATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCt                                                    >  2:3078422/1‑93 (MQ=255)
                                      cATTAACGTCGCAGGCATCGATATGCACAGAACCATTACGTTCAACTATCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTcat                                                   <  2:4225667/93‑3 (MQ=255)
                                       aTTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCttt                                                  <  2:557987/93‑1 (MQ=255)
                                        ttAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTc                                                 <  2:1622610/93‑1 (MQ=255)
                                        ttAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTc                                                 >  2:3267008/1‑93 (MQ=255)
                                          aaCGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAg                                               <  2:3675094/93‑1 (MQ=255)
                                           aCGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTTTTTCCGt                                              >  2:985763/1‑93 (MQ=255)
                                           aCGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGt                                              >  1:1448983/1‑93 (MQ=255)
                                            cGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGtt                                             <  1:767897/93‑1 (MQ=255)
                                            cGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGtt                                             <  1:2384790/93‑1 (MQ=255)
                                            cGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGtt                                             >  2:444334/1‑93 (MQ=255)
                                             gTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGttt                                            >  1:3256629/1‑93 (MQ=255)
                                             gTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGttt                                            >  1:2674514/1‑93 (MQ=255)
                                             gTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGtgt                                            >  1:3596195/1‑91 (MQ=255)
                                              tCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTc                                           >  1:1005299/1‑93 (MQ=255)
                                              tCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTc                                           >  2:1887612/1‑93 (MQ=255)
                                              gcgcgGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTc                                           <  1:1278487/92‑1 (MQ=255)
                                                 cgGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCggg                                        <  2:2773129/93‑1 (MQ=255)
                                                  gggCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGa                                       >  1:1730218/1‑93 (MQ=255)
                                                    gCATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACg                                     <  1:1465462/93‑1 (MQ=255)
                                                     cATCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGt                                    <  1:4356008/93‑1 (MQ=255)
                                                       tCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTc                                  <  2:2276224/93‑1 (MQ=255)
                                                       tCAATATGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTc                                  >  1:3105884/1‑93 (MQ=255)
                                                          atatGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTgg                               >  2:1974947/1‑93 (MQ=255)
                                                           tatGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGa                              >  1:867985/1‑93 (MQ=255)
                                                           tatGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGa                              >  2:1921643/1‑93 (MQ=255)
                                                            atGCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGAt                             >  2:4091787/1‑93 (MQ=255)
                                                              gCACAGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGAtct                           <  1:300466/93‑1 (MQ=255)
                                                                acaGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGGCTTTCAGTTTCGGGACGTTCTGGATctct                         >  2:3443659/1‑93 (MQ=255)
                                                                 caGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTa                        >  1:2792191/1‑93 (MQ=255)
                                                                  aGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTAc                       <  2:4067303/93‑1 (MQ=255)
                                                                  aGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTAc                       <  1:523805/93‑1 (MQ=255)
                                                                  aGAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTAc                       >  2:4425691/1‑93 (MQ=255)
                                                                   gAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTAcc                      >  2:3796360/1‑93 (MQ=255)
                                                                   gAACCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGAGTTATCGACGTCTTTAAGTTTCGGGAGGTTTTGGATCTCTAcc                      >  1:2765485/1‑93 (MQ=255)
                                                                   gAACCATTACGTTCTAATTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTAcc                      >  1:738471/1‑93 (MQ=255)
                                                                    aaCCATTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCg                     >  1:2062538/1‑93 (MQ=255)
                                                                        aTTACGTTCTACTTTCGCGCCTAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGttc                 <  2:2574396/93‑1 (MQ=255)
                                                                        aTTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGttt                 >  1:3332653/1‑92 (MQ=255)
                                                                        aTTACGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGttc                 >  1:1781932/1‑93 (MQ=255)
                                                                           aCGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGttcttc              <  1:565208/93‑1 (MQ=255)
                                                                            cGTTCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTcc             <  1:2770889/93‑1 (MQ=255)
                                                                              ttCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGc           >  2:4241517/1‑93 (MQ=255)
                                                                                cTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCAATGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCa         <  1:4572000/93‑1 (MQ=255)
                                                                                 tACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCag        <  2:2879406/93‑1 (MQ=255)
                                                                                 tACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCag        >  2:4583605/1‑93 (MQ=255)
                                                                                   ctctCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCagta      <  2:104357/90‑1 (MQ=255)
                                                                                   cTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCagta      <  2:2108054/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: ‑