Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,669,340 A→G 100% 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,3400AG98.7% 324.8 / ‑5.4 104intergenic (‑152/‑252)yhjA/treFputative cytochrome C peroxidase/cytoplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base A (4/0);  new base G (59/41);  total (63/41)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

GCAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCGGTA  >  NC_000913/3669244‑3669438
                                                                                                |                                                                                                  
gcAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTa                                                                                                  >  2:2111421/1‑99 (MQ=255)
gcAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAg                                                                                                 >  2:2162667/1‑100 (MQ=255)
gcAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAg                                                                                                 >  2:1403571/1‑100 (MQ=255)
gcAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATAAATATAAACTATCAATAAGGTa                                                                                                  >  1:2807541/1‑99 (MQ=255)
 cAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGc                                                                                                >  2:2996398/1‑100 (MQ=255)
   aatgtttctTCGAGACAAAGCCATCATAAAAATTAGGGAAATATGCTTATTGATAGGAATCAAACAATTACATGAATATAAACTATCAATAAGGTAGc                                                                                                <  1:379796/92‑1 (MQ=255)
     cGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATAACATAAATATAAACTATCAATAAGGTAGCCtgt                                                                                            >  1:77068/1‑100 (MQ=255)
       tctcTGCGAGACAAAGCCATTATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGc                                                                                          <  1:2837046/100‑1 (MQ=255)
        ctctGCGAGACAAAGCCATCATAAAAAACAGGGAGCTACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGcc                                                                                         >  1:1104494/1‑100 (MQ=255)
        ctctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGcc                                                                                         <  2:53104/100‑1 (MQ=255)
         tctGCGATACAAAGCCATCATAAAAAACAGGTAGATACGCTTTTTGATAGGAATCCATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCt                                                                                        <  2:2419797/100‑1 (MQ=255)
          ctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTa                                                                                       >  2:1383503/1‑100 (MQ=255)
          ctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTa                                                                                       >  2:1208100/1‑100 (MQ=255)
          ctGCGAGACAAAGACATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTa                                                                                       <  2:1695680/100‑1 (MQ=255)
            gCGAGATAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTAtt                                                                                     <  2:2650890/100‑1 (MQ=255)
            gCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTAtt                                                                                     >  1:2301336/1‑100 (MQ=255)
            gCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAAATATCAATAAGGTAGCCTGTGCATAtt                                                                                     >  2:1915061/1‑100 (MQ=255)
                gaCAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTAccc                                                                                 <  1:961023/100‑1 (MQ=255)
                gaCAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTAccc                                                                                 <  1:1796561/100‑1 (MQ=255)
                 aCAAAGCCATCATAAAAAACAGTGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCa                                                                                <  1:2905718/100‑1 (MQ=255)
                 aCAAAGCCATCATAAAAAACAGTGAGAGACGCTTCATGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTAcc                                                                                  <  2:1351898/98‑1 (MQ=255)
                 aCAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCa                                                                                <  2:754488/100‑1 (MQ=255)
                  cAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACTCAg                                                                               >  2:694062/1‑100 (MQ=255)
                   aaaGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGa                                                                              <  1:377551/100‑1 (MQ=255)
                   aaaGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATAAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGa                                                                              >  2:258126/1‑100 (MQ=255)
                    aaatcatcatATTAAACATGGATAAACGCTTTTTTATAGGTATCAATCAATTATATGAATATAAACTTTCAATGAGGTAGCCTGTGCCTATTACCTAg                                                                               <  2:1387662/94‑1 (MQ=255)
                      gCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCa                                                                           >  1:547747/1‑100 (MQ=255)
                        catcatAAAAAACAGGGAGATACGCTTTTTGATAGTAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAgg                                                                         <  1:2523942/100‑1 (MQ=255)
                         aaaatAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAg                                                                          <  2:101913/95‑1 (MQ=255)
                          tcatAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGaa                                                                       <  1:2871911/100‑1 (MQ=255)
                            atAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAAt                                                                      >  1:2897557/1‑99 (MQ=255)
                            atAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATg                                                                     >  2:1700736/1‑100 (MQ=255)
                            atAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAACCTATGCCTATTACCCAGATCAGGAATg                                                                     >  1:418171/1‑100 (MQ=255)
                            atAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGaa                                                                       >  1:557949/1‑98 (MQ=255)
                             tAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATgg                                                                    >  1:2863494/1‑100 (MQ=255)
                              aaaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGc                                                                   >  1:550209/1‑100 (MQ=255)
                              aaaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCAAGATCGGGaaatg                                                                    >  2:2311501/1‑96 (MQ=255)
                               aaaaaCAGGGAGATATGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCa                                                                  <  1:2996431/100‑1 (MQ=255)
                               aaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCa                                                                  <  2:550206/100‑1 (MQ=255)
                               aaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCa                                                                  <  1:1700756/100‑1 (MQ=255)
                                aaatCAGGGAGATACGTTTTTTGATATGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCa                                                                  <  1:1496956/99‑1 (MQ=255)
                                aaaaCAGGGAGATACGCTTTTTGATAGGACTCAATCAATTACATGAATATAAACAATCAATAAAATAACATGTACCTATTACACAGATCAGCAATGTCat                                                                 >  1:581483/1‑100 (MQ=255)
                                aaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCat                                                                 >  1:2857079/1‑100 (MQ=255)
                                 aaaCAGGGAGATACGCTTTTTGATAGGAATCAATCGATTACATGAATATAAACTACCAATAAGGTAGCCTGTGCCTATTTCCCAGATCAGGAATGGCa                                                                  <  2:1463914/98‑1 (MQ=255)
                                   aCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCatat                                                               >  1:3009163/1‑99 (MQ=255)
                                          gaTACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTg                                                       >  1:2195046/1‑100 (MQ=255)
                                            tACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcg                                                     >  2:342783/1‑100 (MQ=255)
                                             aCGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgc                                                      >  1:1731065/1‑98 (MQ=255)
                                              cGCTTTTTGATAGGAATCAATCAATTACATGCATATAATCTATCAATAAGATAGCCTATACCTATTACAAAAATCAGAAATGGCATATCATGATTgcagc                                                   >  1:1109797/1‑100 (MQ=255)
                                               gCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcggcg                                                  <  2:559792/100‑1 (MQ=255)
                                                    ttGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGc                                             <  2:498881/100‑1 (MQ=255)
                                                    ttGATAGGAATCAATCAATTACATGAATATAAACGATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGc                                             <  2:77068/100‑1 (MQ=255)
                                                     tGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCa                                            >  2:61601/1‑100 (MQ=255)
                                                     tGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCa                                            <  1:2478659/100‑1 (MQ=255)
                                                     tGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAAATCAGGAATGGCATATCATGATAGCGGCGAGAGCa                                            >  1:1620976/1‑100 (MQ=255)
                                                       aTAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAgg                                          <  2:2301310/100‑1 (MQ=255)
                                                       aTAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAgg                                          >  2:1553049/1‑100 (MQ=255)
                                                       aTAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAgg                                          >  2:1553293/1‑100 (MQ=255)
                                                          ggAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGAtt                                       >  1:1293363/1‑100 (MQ=255)
                                                          ggAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGAtt                                       >  2:1579188/1‑100 (MQ=255)
                                                            aatcaatcAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCATATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGa                                     >  2:2125592/1‑100 (MQ=255)
                                                               caatcaatTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATg                                  >  1:1884401/1‑100 (MQ=255)
                                                               caatcaatTACATAAATATAAACTATCAATACGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATg                                  >  1:2041869/1‑100 (MQ=255)
                                                                aatcaatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACACAGATCAGGAATGGCATATCATGATTGCGGCGATAGCAGGATTGaa                                    >  2:134297/1‑97 (MQ=255)
                                                                aatcaatTACATGAATATAAACTATCAATAAGGTAAACTGTGCATATTACCCAGATCAGGAATGGCATATCATGATTACGGCGAGAGAAGGATTGAATGt                                 >  2:1376932/1‑100 (MQ=255)
                                                                 atcaatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGtt                                >  2:1317316/1‑100 (MQ=255)
                                                                   caatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCt                              >  2:2313352/1‑100 (MQ=255)
                                                                    aatTACATGCATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGGTCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCt                              <  2:2897525/99‑1 (MQ=255)
                                                                    aatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtg                             <  1:2701016/100‑1 (MQ=255)
                                                                     atTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtga                            <  1:1583159/100‑1 (MQ=255)
                                                                      ttACATAAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGACATATCATGATTGCGGCGAGAGCACGATTGAATGTTCt                              <  2:1871851/97‑1 (MQ=255)
                                                                       tACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaa                          <  1:1244148/100‑1 (MQ=255)
                                                                       tACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaa                          >  2:420921/1‑100 (MQ=255)
                                                                        aCATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaa                         >  1:1848893/1‑100 (MQ=255)
                                                                         cttGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTAACCAGATCAGGAATGACATATCATGATTGTGTCGAAAGCAGGAGTGAATGTTGtgaaaa                         >  2:2183592/3‑99 (MQ=255)
                                                                         cATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                        <  2:2358180/100‑1 (MQ=255)
                                                                           tGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatga                      >  2:1782778/1‑100 (MQ=255)
                                                                           tGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatga                      >  2:2324608/1‑100 (MQ=255)
                                                                           tGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatga                      >  2:2736915/1‑100 (MQ=255)
                                                                            ggatatatACTATCAATAAGGTAACCTGTCTCTATTACCCAGATCAGGAATTGCATATCATGATTGCTGCGAGTGCAGGATTGAATGTTCtgataat                        >  2:685451/3‑97 (MQ=255)
                                                                            gAATATAAACTATCAATAAGGTAACCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGAGTGCTGCGAGCGCATGATTGAATGTTCtgaaaatgaa                     <  2:2800081/100‑1 (MQ=255)
                                                                              atatAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatgaaaa                   <  1:1136223/100‑1 (MQ=255)
                                                                              atatAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatgaaaa                   <  1:1015388/100‑1 (MQ=255)
                                                                              atatAAACTATCAATAAGGTAGCCTGTACCTATTACCAAGATAAGGAATGGCATATCATGATTGCGGCGATAGCAGGATTGTATGTTCtgaaaatgaaaa                   >  2:1412870/1‑100 (MQ=255)
                                                                                 taAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGTAGGATTGAATGTTCTGAAAATGAAAaaca                <  2:1607406/100‑1 (MQ=255)
                                                                                 taAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaaca                <  1:2367416/100‑1 (MQ=255)
                                                                                  aaaCTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaacaa               >  1:1899623/1‑100 (MQ=255)
                                                                                  aaaCTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaacaa               >  1:1784090/1‑100 (MQ=255)
                                                                                  aaaCTATCAATAAGGTAGCCTGTGCCTAGTACCCAGATCAGGAATGGCATATCATGATTGCGGCGGGAGCAGGATTGAATGTTCTGAAAATGAAAaaca                <  2:1816481/99‑1 (MQ=255)
                                                                                   aaCTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaacaac              >  1:828823/1‑100 (MQ=255)
                                                                                     cTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAacac            >  1:1423612/1‑100 (MQ=255)
                                                                                     cTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGAGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAacac            >  1:475527/1‑100 (MQ=255)
                                                                                      tATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACAcc           <  2:1293351/100‑1 (MQ=255)
                                                                                         cAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCaaa        >  1:1292917/1‑100 (MQ=255)
                                                                                         cAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCaaa        >  2:1071305/1‑100 (MQ=255)
                                                                                         cAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCaaa        >  1:2889150/1‑100 (MQ=255)
                                                                                          aaTAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAAt       >  2:2459199/1‑100 (MQ=255)
                                                                                          aaTAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAAt       >  1:1087514/1‑100 (MQ=255)
                                                                                           aTAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATc      <  1:1579206/100‑1 (MQ=255)
                                                                                                aTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACAAAATCGGTa  >  2:2553985/1‑99 (MQ=255)
                                                                                                |                                                                                                  
GCAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCGGTA  >  NC_000913/3669244‑3669438

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