Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,669,340 A→G 60.0% intergenic (‑152/‑252) yhjA ← / → treF putative cytochrome C peroxidase/cytoplasmic trehalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,669,3400AG60.0% 93.2 / 150.7 138intergenic (‑152/‑252)yhjA/treFputative cytochrome C peroxidase/cytoplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base A (33/22);  new base G (45/37);  total (78/60)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 139 total aligned reads displayed.

CTGCAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCGGTAT  >  NC_000913/3669242‑3669439
                                                                                                  |                                                                                                   
   cAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGTAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGc                                                                                                 >  1:3240866/1‑100 (MQ=255)
    aaaCGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGcc                                                                                                <  1:3129886/100‑1 (MQ=255)
        gTCTCTGTGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCtgtg                                                                                            <  2:2344893/100‑1 (MQ=255)
        gTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCtgtg                                                                                            <  1:1108326/100‑1 (MQ=255)
         tctcTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGc                                                                                           >  2:558029/1‑100 (MQ=255)
          ctctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGcc                                                                                          <  2:1658064/100‑1 (MQ=255)
           tctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCt                                                                                         <  1:426594/100‑1 (MQ=255)
            ctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTa                                                                                        >  2:925180/1‑100 (MQ=255)
            ctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTa                                                                                        >  2:3015920/1‑100 (MQ=255)
            ctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTa                                                                                        <  1:1750366/100‑1 (MQ=255)
            ctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTa                                                                                        >  1:19136/1‑100 (MQ=255)
            ctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCt                                                                                         <  1:2632831/99‑1 (MQ=255)
            ctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTa                                                                                        >  2:1398125/1‑100 (MQ=255)
             tGCGAGACAAAGCCATCATAAAAAACAGGGGGATTCGCTTTTTGATAGGAAGCAATCAATTGCATGAATATAAACTATCAATAAGATAGCCTGTGCCTa                                                                                        <  2:1505467/99‑1 (MQ=255)
             tGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATCGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTAt                                                                                       >  2:2890312/1‑100 (MQ=255)
             tGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGATTATAAACTATCAATAAGATAGCCTGTGCCTAt                                                                                       <  2:3213656/100‑1 (MQ=255)
             ngCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTAt                                                                                       >  2:1942955/2‑100 (MQ=255)
              gCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTAtt                                                                                      >  1:518179/1‑100 (MQ=255)
              gCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTAtt                                                                                      >  2:476497/1‑100 (MQ=255)
              gCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTAtt                                                                                      >  2:672548/1‑100 (MQ=255)
              gCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTAtt                                                                                      <  2:2217807/100‑1 (MQ=255)
               cGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTa                                                                                     <  1:1428668/100‑1 (MQ=255)
               cGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTa                                                                                     <  1:23769/100‑1 (MQ=255)
               cGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTa                                                                                     >  2:779015/1‑100 (MQ=255)
               cGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTa                                                                                     >  1:2974675/1‑100 (MQ=255)
                   tcAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGTTAGCCTGTGCCTATTACCCa                                                                                 <  1:2071097/99‑1 (MQ=255)
                    cAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTTTTTTTTTTTAATCAATTACATGAATATAAACTATGAATAAGGTAGCCTGTGCCTATTACCCa                                                                                 >  1:411566/1‑99 (MQ=255)
                    cAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCaa                                                                                >  2:2071110/1‑99 (MQ=255)
                    cAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAg                                                                                >  1:1640188/1‑100 (MQ=255)
                    cAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACACAg                                                                                >  1:1537338/1‑100 (MQ=255)
                      aaGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACAAGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGAt                                                                              >  1:643369/1‑100 (MQ=255)
                       aGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATc                                                                             >  2:1175056/1‑100 (MQ=255)
                       aGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATc                                                                             <  2:2419812/100‑1 (MQ=255)
                       aGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATc                                                                             >  2:1210473/1‑100 (MQ=255)
                        gCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCa                                                                            >  1:3161297/1‑100 (MQ=255)
                          catcatAAAAAACAGGGAGATACGATTTATAATAGGAATAAATCAATTACATGAATATACACTATCAATAAGATAGCATGTGCCTATTACCCAGATCAg                                                                           <  1:628601/99‑1 (MQ=255)
                           aaaatAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAg                                                                           <  1:703236/95‑1 (MQ=255)
                            tcatAAAAAATAGGGAGATACGCTTTTTGATAGGAATCAATAAATTACATGAATATAAACTATAAATAAGATAGCATGTGCCTATTACCCAGATAAGGaa                                                                        >  1:1707795/1‑100 (MQ=255)
                            tcatAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGaa                                                                        >  1:2606254/1‑100 (MQ=255)
                              atAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCATGTGCCTATTACCCAGATAAgca                                                                         >  1:3201216/1‑95 (MQ=255)
                              atAAAAAACAGGGAGATACGCTTTTTGACAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATg                                                                      <  2:1060800/100‑1 (MQ=255)
                                aaaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGc                                                                    >  2:956548/1‑100 (MQ=255)
                                aaaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGc                                                                    >  1:1752253/1‑100 (MQ=255)
                                aaaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGc                                                                    <  2:476622/100‑1 (MQ=255)
                                 aaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCa                                                                   <  2:2512607/100‑1 (MQ=255)
                                 aaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCa                                                                   >  2:694497/1‑100 (MQ=255)
                                  aaaaCAGGGAGATACTCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCat                                                                  <  1:2466356/100‑1 (MQ=255)
                                   aaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATGAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCata                                                                 >  1:2166253/1‑100 (MQ=255)
                                    aaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCatat                                                                <  2:1938456/100‑1 (MQ=255)
                                      cAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCa                                                              <  2:855425/100‑1 (MQ=255)
                                      cAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATc                                                               >  1:2595950/1‑99 (MQ=255)
                                        gggAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATg                                                            >  2:2313587/1‑100 (MQ=255)
                                         ggAGATACGCTTTTTTATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGa                                                           <  1:2194802/100‑1 (MQ=255)
                                          gagaTACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGAt                                                          >  1:1606147/1‑100 (MQ=255)
                                              tACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcg                                                      <  2:540251/100‑1 (MQ=255)
                                              tACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcg                                                      <  2:2531277/100‑1 (MQ=255)
                                               aCGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATCGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcgg                                                     >  1:1475952/1‑100 (MQ=255)
                                                 gCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcggcg                                                   <  1:1229158/100‑1 (MQ=255)
                                                  cTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcggcg                                                   >  2:3081483/1‑99 (MQ=255)
                                                  cTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCga                                                  >  2:3080543/1‑100 (MQ=255)
                                                      ttGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGc                                              <  1:1220101/100‑1 (MQ=255)
                                                      ttGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGc                                              <  1:1619225/100‑1 (MQ=255)
                                                       tGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCa                                             >  1:3043213/1‑100 (MQ=255)
                                                        gATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAg                                            <  2:2560959/100‑1 (MQ=255)
                                                         aTAGGAATCAATCAATTACATGAATATAAACTATAAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCaag                                           >  2:2273949/1‑98 (MQ=255)
                                                           aGGAATCAATCAATTACATGAATATAAACTATCGATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAgg                                           >  1:2241603/1‑98 (MQ=255)
                                                           aGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGAt                                         >  1:3317430/1‑100 (MQ=255)
                                                            ggAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGAtt                                        >  2:418311/1‑100 (MQ=255)
                                                            ggAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGAtt                                        >  1:3352037/1‑100 (MQ=255)
                                                             gAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTg                                       <  1:2402210/100‑1 (MQ=255)
                                                             gAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTg                                       <  1:3015894/100‑1 (MQ=255)
                                                              aatcaatcAATTACATGAATATAAACTATCAATAAGGTAGACTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGa                                      <  1:1394741/100‑1 (MQ=255)
                                                               atcaatcaATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGaa                                     <  2:2794229/100‑1 (MQ=255)
                                                                tcaatctaTTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAAt                                    <  1:2903592/100‑1 (MQ=255)
                                                                 caatcaatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAAt                                    <  1:856847/99‑1 (MQ=255)
                                                                     caatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCt                               <  1:2789072/100‑1 (MQ=255)
                                                                     caatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCt                               >  2:2482965/1‑100 (MQ=255)
                                                                        ttACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaa                            <  2:2602446/100‑1 (MQ=255)
                                                                        ttACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaa                            <  1:542809/100‑1 (MQ=255)
                                                                           cATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                         >  2:375434/1‑100 (MQ=255)
                                                                           cATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                         >  2:1458007/1‑100 (MQ=255)
                                                                           cATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                         >  2:2655040/1‑100 (MQ=255)
                                                                             tGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatga                       >  2:2785831/1‑100 (MQ=255)
                                                                             tGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatga                       >  2:75798/1‑100 (MQ=255)
                                                                               aaTATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatgaaa                     >  2:3199464/1‑100 (MQ=255)
                                                                                atatAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatgaaaa                    <  1:1535939/100‑1 (MQ=255)
                                                                                 tataAACTATCAATAAGGTAGACTGTGCCTATTATCCAGATCAGGAATAGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAGTGaaaaa                   <  1:1782078/100‑1 (MQ=255)
                                                                                 tataAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGaaaaa                   >  2:2599842/1‑100 (MQ=255)
                                                                                 tataAACCATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGaaaaa                   >  1:1610828/1‑100 (MQ=255)
                                                                                  ataAGCGATCAATAAGGTAGCCTGTGCCGATTACCCGGAGCAGGAATGGCATAGCAGGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaac                  <  1:2601777/100‑1 (MQ=255)
                                                                                     aaCTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaacaac               >  1:3193343/1‑100 (MQ=255)
                                                                                     aaCTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaacaac               <  1:502442/100‑1 (MQ=255)
                                                                                       cTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAacac             <  2:2365971/100‑1 (MQ=255)
                                                                                          tCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCaa          >  2:2756451/1‑100 (MQ=255)
                                                                                           cAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCaaa         >  1:1871489/1‑100 (MQ=255)
                                                                                           cAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCaaa         >  2:2807583/1‑100 (MQ=255)
                                                                                            aaTAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAAt        <  2:6153/100‑1 (MQ=255)
                                                                                            aaTAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAAt        >  2:1681229/1‑100 (MQ=255)
                                                                                                 gATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCGGTa   >  2:3373372/1‑100 (MQ=255)
                                                                                                  aTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCGGTAt  <  1:2577258/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
CTGCAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCGGTAT  >  NC_000913/3669242‑3669439

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