Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,4410AC46.6% 115.2 / 25.5 73G384G (GGT→GGGyqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (39/0);  new base C (0/34);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.35e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTG  >  NC_000913/2984292‑2984582
                                                                                                                                                      |                                                                                                                                             
gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCacac                                                                                                                                              >  2:321724/1‑151 (MQ=255)
    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTg                                                                                                                                          >  1:579051/1‑151 (MQ=255)
                   aaaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATAc                                                                                                                            >  1:86344/1‑150 (MQ=255)
                     aaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACaa                                                                                                                          >  1:615753/1‑150 (MQ=255)
                     aaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGCCACTTTAAATGGTCAATGCCCCCCCCTGCCCCCCGCCGATACaaa                                                                                                                         <  2:125500/151‑1 (MQ=255)
                     aaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGAAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACa                                                                                                                           >  2:91442/1‑149 (MQ=255)
                       aTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGc                                                                                                                       >  2:742585/1‑151 (MQ=255)
                           gggAGCCCAAAAAATGTATGTCGACTTCCCAGCCATTATTATTCAACAGATTAATGACGGCAGTTATTTTTCGGAAAGGAGAAACTGCAGGGTTTTTCGTCCCTTTCATTGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGtt                                                                                                                    <  1:282179/150‑1 (MQ=255)
                           gggAGCCAAAAAAATGTATGTCGATTTGCCACCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGAAGGGAGAAAATGCAGAGCTATTCGTCACGTTAAATGGTCAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTg                                                                                                                   <  2:390141/151‑1 (MQ=255)
                            ggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGc                                                                                                                  >  2:186491/1‑151 (MQ=255)
                                cccAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGt                                                                                                              >  2:601063/1‑151 (MQ=255)
                                ccaaaaaaaTGTATGTCGACTTGCCACCCATAATTATTAAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGAAGAGCTTTTCGTCACGTTAAATGGTCAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGcc                                                                                                                 <  1:788003/143‑1 (MQ=255)
                                   aaaaatTGTATGTCGATTTGCCAGCCATCATTATCCAACAAATTAAGGCCGGAAAATATTTATCGGAAAGGAGAAAATGCAGAGCTTTTCGTCACGTAAAATGGTAAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAg                                                                                                               <  1:551267/145‑1 (MQ=255)
                                   aaaaaaTGTATGTCTACTTGCCAGCCATCTTTATTAAACAGATTAATGACGGAAGATTTTTATCGGAAAGGGGAAAATGCAGAGTTTTTCGTCACTTTCAAGGGTCAATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                           <  1:251363/149‑1 (MQ=255)
                                     aaaaTTTATCTATACTTCCCAGCAATAATTATTCAAAAAATTAATGACGGCAGTTTTTTAGCGGTAAGGAGAAAACACAGATCTTTTCGTCACTTTAAAGGGTAACCCCCACCCCCAGCCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                         <  2:397637/151‑1 (MQ=255)
                                           tatgtCGACTTGCCACCCATAATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACTTTCAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCggg                                                                                                    <  1:540491/150‑1 (MQ=255)
                                            atgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGc                                                                                                   >  2:202056/1‑150 (MQ=255)
                                             tgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAt                                                                                                 >  1:14902/1‑151 (MQ=255)
                                              gtCGACTTGCCAGCCATCTTTATTCAACAGATTAATGACGGCAGATATTTATCGGAAAGGAGAAAATGCAGAGCTTTTCGTCCCGTTCAATGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAt                                                                                                 <  1:347984/150‑1 (MQ=255)
                                                cgctttGCCAGCCATCATTATTCAACAAATAAATGACGGCAGATATTTATGGGTAAGGAGAAAAGGCAGAGCTTTTCGTCACTTTAAAGGTCCAACCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc                                                                                              <  1:410668/147‑1 (MQ=255)
                                                cGACTTGCCAGCCATCATTATTAAACAGATAAAGGACGGCATATATTTATCGGAAAGGAGAAAATGCAGAGCTATTCGTCACGTTAAAGGGTCAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                                                                <  2:724860/149‑1 (MQ=255)
                                                cGACTTGCCACCCATCATTATTCAACAGATTAAGGCCGGCAGATTTTTATCGGTAAGGAGAAAAGGCAGACCTATTCGTCCCGTTAAATGGCCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTa                                                                                               <  2:119358/150‑1 (MQ=255)
                                                    ttGCCAGCCATCATTATTCAAAAGATAAATGACGGCAGATATTTATGGGAAAGGAGAAAATGCAGAGCTTTTCGTCACGTTCAATGGTAAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGt                                                                                            <  2:264882/149‑1 (MQ=255)
                                                     tGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                         >  1:409162/1‑151 (MQ=255)
                                                       ccagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAc                                                                                       >  2:740326/1‑151 (MQ=255)
                                                       ccagccaTAATTAATCAACATATTAATGACGGAAGATATTTACCGGTAAGGAAAAAATACAAAGCTCTCGGCCTCTTAGGAAGCGAAAATCAACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt                                                                                          <  1:196703/148‑1 (MQ=255)
                                                        cagccaTTTTTATCCAACAGATTAATGACGGAAAATATTTATCGGAAAGGAGGAAATGCAGAGCTTTTCTTCCCTTTCAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                         <  2:202379/148‑1 (MQ=255)
                                                           ccaTCTTTTTACAACAGATAAATGACGGCAGATTTTTATCGGAAAGGAGAAAATGCAGAGCTTTTCGTCCCTTTAAAGGTTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAt                                                                                    <  1:523380/150‑1 (MQ=255)
                                                           ccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt                                                                                   >  2:397185/1‑151 (MQ=255)
                                                            catcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt                                                                                  >  1:263146/1‑151 (MQ=255)
                                                                 ttattCAACAGATTAAGGACGGCAGATTTTTAGGGGAAAGGAGAAAAGGCAGAGCTATTCTCCACGTTAAAGGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                                                                <  1:202210/148‑1 (MQ=255)
                                                                  tattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa                                                                             >  2:399036/1‑150 (MQ=255)
                                                                  tattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGTTTTTCGTCACTTTCAATGGTCAATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa                                                                             <  1:453191/150‑1 (MQ=255)
                                                                   attCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                           >  1:601369/1‑151 (MQ=255)
                                                                   attCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                           >  1:455214/1‑151 (MQ=255)
                                                                      cAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                                                                         >  1:36378/1‑150 (MQ=255)
                                                                       aaCAGATTAATGACGGCAGATATTTATCGGAAAGGAGAAAAGGCAGAGTTATTCGTCACTTTCAATGGTCAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                                                                         <  2:36361/149‑1 (MQ=255)
                                                                         cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                        >  1:495587/1‑148 (MQ=255)
                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                                                                  >  2:628550/1‑151 (MQ=255)
                                                                             taaaTGACGGAAGATATTTAGCGGAAAGGAGGAAAGGAAGAGTTTTTCGTCCCGTTAAAGGGTAAAGGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaa                                                                    <  2:760546/146‑1 (MQ=255)
                                                                               aagaaCGGAAAATATTTACCGGTAAGGAGAAAATGCAGAGTTTTTCGTCACTTTAAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAg                                                               <  2:815308/147‑1 (MQ=255)
                                                                               aaTGGCGGAAAATATTTATCGGAAAGGAGAAAATGCAGAGCTTTTCGCCAGTTTAAATGGCCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                           <  1:97190/139‑1 (MQ=255)
                                                                                aTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                              >  1:100039/1‑151 (MQ=255)
                                                                                  gACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt                                                            >  1:241412/1‑151 (MQ=255)
                                                                                     ggaaGATATTTATCGGAAAGGAGAAAATGAAGAGTTTTTCGTCCCGTTCAATGGTCAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                         <  1:628599/148‑1 (MQ=255)
                                                                                     ggCAGATATTTATCGGTAAGGAGGAAAGGCAGATCTTTTCGTCACTTTAAATGGTCAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc                                                          <  1:612236/150‑1 (MQ=255)
                                                                                     ggCAAATTTTTATGGGAAAGGAGAAAAGGCAGAGCTTTTCGTCACTTTAAATGTTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                 <  1:647751/143‑1 (MQ=255)
                                                                                       cAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                       >  2:362649/1‑151 (MQ=255)
                                                                                             tgtATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTGTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCg                                                  >  2:159648/3‑150 (MQ=255)
                                                                                                aTCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                >  1:591065/1‑149 (MQ=255)
                                                                                                 tCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc                                             >  1:673519/1‑151 (MQ=255)
                                                                                                  cGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                            >  2:225132/1‑151 (MQ=255)
                                                                                                  cGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                            >  1:229470/1‑151 (MQ=255)
                                                                                                   ggTAAGGAGAAAATGCAGCGTTTTTCGTCACGTTAAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg                                              <  2:204277/148‑1 (MQ=255)
                                                                                                      aaGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                         >  2:647435/1‑150 (MQ=255)
                                                                                                           ggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg                                   >  2:61497/1‑151 (MQ=255)
                                                                                                             aaaTGCAGAGTTTTTCGCCACGTTCAATGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                    <  2:232197/148‑1 (MQ=255)
                                                                                                             aaaTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                 >  1:180658/1‑151 (MQ=255)
                                                                                                                  cAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                             >  1:292378/1‑150 (MQ=255)
                                                                                                                      gtttttcttcACGTTAATGGTTAACTCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt                          <  1:739465/145‑1 (MQ=255)
                                                                                                                      gtttttCGTCACGTTAAATGTTCAAGCCCCCCCCGTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt                          <  2:816576/145‑1 (MQ=255)
                                                                                                                        ttttCGTCACGTTAAAGGGCCAACCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc                        <  1:376532/147‑1 (MQ=255)
                                                                                                                         aTTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGAt                     >  1:419815/1‑151 (MQ=255)
                                                                                                                           tCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                   >  2:351278/1‑151 (MQ=255)
                                                                                                                            cGTCACTTTCAATGGAAAACCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                   <  1:496942/150‑1 (MQ=255)
                                                                                                                            cGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc                  >  1:680112/1‑151 (MQ=255)
                                                                                                                             ggcacgtaaactgctCAATGCCCCCCCCTGCCACCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTGCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                    <  2:88074/134‑1 (MQ=255)
                                                                                                                              tCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCgcgc                >  2:524660/1‑151 (MQ=255)
                                                                                                                                aCGTTAAATGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcca              <  1:129602/151‑1 (MQ=255)
                                                                                                                                    tCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcat          >  2:677120/1‑151 (MQ=255)
                                                                                                                                      aaTGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg        >  2:812229/1‑151 (MQ=255)
                                                                                                                                      aaTGGTAAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacca           <  2:371229/148‑1 (MQ=255)
                                                                                                                                       aTGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg        >  1:512843/1‑150 (MQ=255)
                                                                                                                                            cAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg  >  2:664334/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                             
GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTG  >  NC_000913/2984292‑2984582

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