Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,859,4320CT79.7% 179.7 / ‑3.2 64G780G (GGC→GGTmutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (16/35);  minor base G (13/0);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 5.17e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCT  >  NC_000913/2859288‑2859578
                                                                                                                                                |                                                                                                                                                  
gACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCACCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCa                                                                                                                                              <  2:48623/151‑1 (MQ=255)
  cccAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCaaa                                                                                                                                            <  2:120476/151‑1 (MQ=255)
      gTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAagag                                                                                                                                        <  2:607916/151‑1 (MQ=255)
        tACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTTCAGCTCTTGCAGGTGTGCCGAAAGAgg                                                                                                                                       >  2:430675/1‑150 (MQ=255)
         aCCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGGAGCTCTGGGAGGGGGGGCAAAAAAggg                                                                                                                                      >  1:626570/1‑149 (MQ=255)
             gagaAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGttatta                                                                                                                                 <  2:232861/151‑1 (MQ=255)
              agaAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTaa                                                                                                                                <  1:533943/151‑1 (MQ=255)
                     ggAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCa                                                                                                                         <  1:146190/151‑1 (MQ=255)
                       aaGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGGGTGCCCAAAGAGGTTATTAAGCGCGCAc                                                                                                                        >  1:969107/1‑150 (MQ=255)
                       aaGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGGGGGCCAAAAGAGGGTATTAAGCGCGCa                                                                                                                         >  2:186620/1‑149 (MQ=255)
                          gCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGGTATTAAGCGCGCACGGc                                                                                                                     >  2:223711/1‑150 (MQ=255)
                              cGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAg                                                                                                                <  1:869469/151‑1 (MQ=255)
                                  aaCGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGGGGGCCCAAAGAGGGTATTAAGCGCGCCCGGCAAAAGCTg                                                                                                             >  1:1051296/1‑150 (MQ=255)
                                        cATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGc                                                                                                      <  1:92795/151‑1 (MQ=255)
                                         aTCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGGGTGCCAAAAGAGGGTATTAAGCGCGCACGGCAAAAGCTGCGTGAg                                                                                                       >  2:298686/1‑149 (MQ=255)
                                             cGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTgg                                                                                                   >  2:168986/1‑149 (MQ=255)
                                             cGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGa                                                                                                  <  1:604940/150‑1 (MQ=255)
                                             cGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGGTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTgg                                                                                                   >  1:327783/1‑149 (MQ=255)
                                               aTGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAg                                                                                               <  2:106887/151‑1 (MQ=255)
                                                tGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGGTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGaaa                                                                                                >  2:516684/1‑149 (MQ=255)
                                                 gCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGGTATTAAGCGCGCACGGCCAAAGCTGCGTGAGCTGGaaa                                                                                                >  1:735822/1‑148 (MQ=255)
                                                  cACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGGGGGCCAAAAGAGGGTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGaaa                                                                                                >  2:999571/1‑147 (MQ=255)
                                                    cTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGGGTGCCAAAAGAGGGTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCAttt                                                                                          >  1:290618/1‑151 (MQ=255)
                                                      ggAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTc                                                                                         <  2:38764/150‑1 (MQ=255)
                                                       gAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGc                                                                                       <  1:852877/151‑1 (MQ=255)
                                                        aGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGGTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGcc                                                                                      >  2:1021751/1‑151 (MQ=255)
                                                                gACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAAcgccg                                                                              <  2:460327/151‑1 (MQ=255)
                                                                      aTTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTa                                                                        <  2:961766/151‑1 (MQ=255)
                                                                          ccTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGGGTGCCCAAAGAGGGTATTAAGCGCGCACGGCCAAAGCTGCGTGAGCTGGAAAACATTTCGCCGAACGCCGCCCCCACGc                                                                     >  1:644484/1‑150 (MQ=255)
                                                                              tATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAg                                                                  <  2:806348/149‑1 (MQ=255)
                                                                               aTGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTgg                                                               <  1:888933/151‑1 (MQ=255)
                                                                               aTGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTg                                                                <  2:375537/150‑1 (MQ=255)
                                                                               aTGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGGGTGCCAAAAGAGGGTATTAAACGCGCACGGCAAAAGCTGCGGGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAgggg                                                               >  1:1023417/1‑151 (MQ=255)
                                                                                 gCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGa                                                              <  2:274942/150‑1 (MQ=255)
                                                                                            aGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGGACGCaaa                                                   >  2:923538/1‑150 (MQ=255)
                                                                                             ggATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCCGGGGTGCCAAAAGAGGATTTTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCCAGTGGATGGTACGCAAAt                                                  >  1:572190/1‑150 (MQ=255)
                                                                                             ggATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAAt                                                  <  1:68138/150‑1 (MQ=255)
                                                                                                tGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCCGCTTTGGCAGGTGGGGCCCAAGAGGGTATTAAGAGCGCACGGCGAAAACATCGTGAGGTGGAAAGAATTTTGGCGAACGCCGCCGCTACCCAAGTGGCTgggaaatcaacagc                                               >  1:958035/1‑138 (MQ=255)
                                                                                                 ggCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCtt                                             <  1:669257/151‑1 (MQ=255)
                                                                                                   cgcgGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGGGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCCACGCAAGTGGATGGTACGCAAAAGTCttt                                            >  1:481847/1‑150 (MQ=255)
                                                                                                    gcggcgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGGAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTTGATGGTACGCAAATGTCTTTGc                                          <  1:56723/151‑1 (MQ=255)
                                                                                                     cggcgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGGAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAACTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGAAGGGACGCAAATGTCTTTTCt                                         >  2:925932/1‑151 (MQ=255)
                                                                                                      ggcgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGGGTGCCCAAAGAGGGTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTTGCCGAACGCCGCCGCTACGCAAGTGGGTGGGACGCAAATGTCttt                                            >  1:744634/1‑147 (MQ=255)
                                                                                                       gcgAGCAAAAGCTACGGCCTGGCGGTTGCCGCTCTGGCAGGGGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGGACGCAAATGTCTTTGCTGt                                       >  2:836685/1‑151 (MQ=255)
                                                                                                       gcgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGt                                       <  2:27815/151‑1 (MQ=255)
                                                                                                       gcgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGt                                       <  1:282261/151‑1 (MQ=255)
                                                                                                       gcgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGGTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCttt                                            >  2:284164/1‑146 (MQ=255)
                                                                                                          aGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAg                                    <  1:483802/151‑1 (MQ=255)
                                                                                                          aGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAg                                    <  1:932154/151‑1 (MQ=255)
                                                                                                            cAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGGGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGGACGCAAATGTCTTTGCTGTCa                                     >  2:847663/1‑148 (MQ=255)
                                                                                                              aaaGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTAcc                                <  1:878238/151‑1 (MQ=255)
                                                                                                                 gCTACGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCaga                             <  1:663840/151‑1 (MQ=255)
                                                                                                                    aCGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCagaag                           <  2:968014/150‑1 (MQ=255)
                                                                                                                    aCGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCagaag                           <  1:350365/150‑1 (MQ=255)
                                                                                                                    aCGGCCTGGCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCagaa                            <  1:337643/149‑1 (MQ=255)
                                                                                                                           ggCGGTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACtt                     <  2:553950/149‑1 (MQ=255)
                                                                                                                               gTTGCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTg               >  2:477644/1‑151 (MQ=255)
                                                                                                                                  gCAGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCg             <  1:183373/150‑1 (MQ=255)
                                                                                                                                    aGCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTc          >  1:930509/1‑151 (MQ=255)
                                                                                                                                     gCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAATCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTc          <  1:43158/150‑1 (MQ=255)
                                                                                                                                     gCTCTGGCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCg         <  1:748832/151‑1 (MQ=255)
                                                                                                                                           gCAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGctc   >  1:741718/1‑151 (MQ=255)
                                                                                                                                            cAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGctct  <  1:76781/151‑1 (MQ=255)
                                                                                                                                            cAGGTGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGctct  >  2:1066270/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCT  >  NC_000913/2859288‑2859578

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