Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,4410AC41.4% 161.7 / 34.3 87G384G (GGT→GGGyqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (46/5);  new base C (0/36);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 2.05e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAC  >  NC_000913/2984292‑2984589
                                                                                                                                                     |                                                                                                                                                    
gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCacac                                                                                                                                                     <  3:1299941/151‑1 (MQ=255)
  atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAccc                                                                                                                                                   <  4:303003/151‑1 (MQ=255)
   tatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCt                                                                                                                                                  >  2:81793/1‑151 (MQ=255)
    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACATATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCt                                                                                                                                                  <  2:55560/150‑1 (MQ=255)
    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCt                                                                                                                                                  <  3:688133/150‑1 (MQ=255)
        cATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGccc                                                                                                                                             >  3:526591/1‑151 (MQ=255)
         aTGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcccc                                                                                                                                            >  3:491419/1‑151 (MQ=255)
            agaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgc                                                                                                                                         <  2:877350/151‑1 (MQ=255)
                   aaaaaTGAGGGAGCCAAAAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACATATTAATGACGGAAATTTTTTATCGGTAAGGAGAAAATGAAGAGCTTTTCGTCCCTTTCAATGGTAAATCCCCCCCCCTGGCCCCCGCCGATACa                                                                                                                                  <  3:997101/151‑1 (MQ=255)
                       aTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGc                                                                                                                              >  3:1307162/1‑151 (MQ=255)
                          aGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGtt                                                                                                                           >  2:518343/1‑151 (MQ=255)
                           gggAGCCTAAAAAATGTTGGCCGACTCCCCACCCATAATTATCCAACAGATTAATGCCGGGAAATTTTTATCGGTAAGGAGAAAAGCCAAGTCTATTCGCCACTTTCAATGATCAATGCCACCCCTTGGCCCCCGCCGATACAAAGCGTTg                                                                                                                          <  4:325909/151‑1 (MQ=255)
                           gggACCCAAAAAAATGTATGTCGCCTTCCCAGCCATCTTTATAAAACAAATAAATGACGGCATATTTTTATGGGTAAGGAGAAAATGCAGAGTTTTTCGTCACGTTAAATGGTAAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTg                                                                                                                          <  3:238904/151‑1 (MQ=255)
                            ggAGCCCAAAAAATGTATGTCGACTTGCCACCCATCATTATTCAACAGATTAATGACGGCAGATATTTAGCGGAAAGGAGAAAATGCAGATCTTTTCGCCACGTTAAATGGTCAAGCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGc                                                                                                                         <  3:178922/151‑1 (MQ=255)
                             gAGCCCAAAAAATTTATGTCGACTTCCCAGCCATTATTATTCAACAGATTAATGACGGCAGATTTTTATGGGAAAGGAGAAAATGCAGATTTTTTCGCCACTTTCAAGGGACAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGc                                                                                                                         <  3:1185818/150‑1 (MQ=255)
                                cccAAAAAGTGTATGTCGCTTTCCCAGCCATCATTATTAAAAAGTTTAAGGCCGGCATATTTTTATGGTAAAGGAAAAAATGAAGAGCTTTTCGTCACGTTAAATGGAAATTGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCa                                                                                                                       <  4:303159/149‑1 (MQ=255)
                                 ccAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTc                                                                                                                    >  2:532337/1‑151 (MQ=255)
                                  cAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTcc                                                                                                                   >  4:418587/1‑151 (MQ=255)
                                   aaaaaaTGATTGTCGTCTTCCCAGCAAACATTATCCAACAAATAAATGGCGGAAGATATTTATCGGAAGGGAGAAAAGGCAGGTCTTTTCGCCCATTTAAAGGATCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTcc                                                                                                                   <  2:367381/148‑1 (MQ=255)
                                     aaaaTGTTTGTCGACTTGCCACCCATCATTATTCAACAGATTAAGGCCGGCAGATATTTAGCGGTAAGGAGAAAATGCAGAGCTTTTCGTCCCTTAAAATGGTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                                <  3:694747/151‑1 (MQ=255)
                                      aaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                                                               >  4:112185/1‑151 (MQ=255)
                                        atgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                             >  4:467803/1‑151 (MQ=255)
                                        atgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                             >  1:87170/1‑151 (MQ=255)
                                         tgtatgcCGACTTGCCAGCCATAATTATTCAAAAGATTAAGGAGGGCAGATATTTAGCGGTAAGGAGAAAATGCAGAGTTATTCGTCACTTTAAATGGTAAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                             <  1:89580/150‑1 (MQ=255)
                                                cGACTTGCCAGCAATCTTTATTCAACAGATTAATGACGGAAGATATTTAGCGGAAAGGAGGAAATGCAGAGTTTTTCGTCACTTTAAAGGGTAAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTa                                                                                                      <  2:999232/150‑1 (MQ=255)
                                                 gACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                                    >  3:228671/1‑151 (MQ=255)
                                                    ttGCCACCAATTATTATTCAAAATATTAATGACGGTAAATATTTTTGGGAAAGGGGAAAAGGCGGAGTTTTTCGTCAGGTCCAATGGAAAAGGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                             <  3:1291588/139‑1 (MQ=255)
                                                     tGCCACCAATTATTATACAAAAAATAAAGGACGGAAAATTTTTATCGGAAAGGAGAAAAGGCAAAGCTTTTCGTAACTTTCAATGGTAAATCCCCCCCCCCGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                                                                                  <  3:1090276/149‑1 (MQ=255)
                                                       ccagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAc                                                                                              >  3:244124/1‑151 (MQ=255)
                                                              tcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                        >  2:627731/1‑150 (MQ=255)
                                                              tcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGATCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                        >  4:255855/1‑150 (MQ=255)
                                                               catTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                      >  3:641676/1‑151 (MQ=255)
                                                                 ttataCAAAAGATAAAGGAGGGAAGATTTTTATGGGAAAGGAGGAAAGGCAGAGTTTTTCGTCACTTTAAATGGTAAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                        <  2:491385/147‑1 (MQ=255)
                                                                  tattCAACAGATAAATGACGGTATATTTGTACCGGTAAGCAACAAATGGCGGCCTATTCGCTTTGTTAAGGGGTCAAGCCCCCACCCAGACCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                                                                       <  3:255177/147‑1 (MQ=255)
                                                                   attCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                                  >  2:413811/1‑151 (MQ=255)
                                                                     tCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                                 >  4:395351/1‑150 (MQ=255)
                                                                      cAACAGATTAAGGACGCCATATCTTTACGGGTAAGAAGGAAATGCAGAGTTTTTCTTCCCTTTCATTGGTCAACGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                                                                                <  2:813565/150‑1 (MQ=255)
                                                                       aaaaGATTAATGACGGCAGATTTTTGGCGGAAGGGAGAAAACGCAGATCTATTCGTCAGTTAAAATGGTCATCCCCCCCCCGTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                               <  3:956308/147‑1 (MQ=255)
                                                                       aaCAGATTAATGACGGAAGATATTTATCGGAAAGGAGAAAATGCAGAGCTTTTCGTCACTTACAATGGTAAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGaa                                                                                   <  2:1367105/146‑1 (MQ=255)
                                                                       aaCAAATTAATGACGGAAGATATTTATCGGTAAGGAGAAAATGCAGAGTTTTTCGTCACTTTAAATGGTAAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                                                                              <  3:1017229/151‑1 (MQ=255)
                                                                         cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                               >  2:867359/1‑148 (MQ=255)
                                                                         aaGATTAATGACGGCAGATATTTAGCGGAAAGGAGGAAAGGCAGAGTTTTTCGTCACTTTAAATGGTCAACGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCagag                                                                             <  3:1064633/149‑1 (MQ=255)
                                                                          aGATTAATGACGGAAGATATTTATGGGTAAGGAGGAAATGCAGGGTTTTTCGTCCCGTTAAATGGTCAATCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                                                                                <  2:228644/146‑1 (MQ=255)
                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                                                                         >  3:248533/1‑151 (MQ=255)
                                                                             ttAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                        >  2:1039686/1‑151 (MQ=255)
                                                                              tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                        >  4:245457/1‑150 (MQ=255)
                                                                               aaggACGGCAGTTTTTTATGGGTAAGGAGAAAAGGCAGAGCTTTTCGTCACGTTAAATGGTCAACCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                        <  3:593586/146‑1 (MQ=255)
                                                                               aaTGACGGAAGATATTTTTCGGTAAGGAGAAAATGAAGATTTTTTCGTCCCTTTAAATGTTCAATGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                       <  2:1248954/150‑1 (MQ=255)
                                                                                 tGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                     >  1:157020/1‑150 (MQ=255)
                                                                                  gACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt                                                                   >  4:343891/1‑151 (MQ=255)
                                                                                    cggaaattttttATGGGAAGGGAGGAAATGCAGATTTTTTCGCCACGTTAAATGGCCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                     <  3:532370/138‑1 (MQ=255)
                                                                                     ggCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACTTTAAATGGTCAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAtt                                                                  <  1:34542/149‑1 (MQ=255)
                                                                                     ggCAAATATTTAGCGGTAGGGAGAAAATGAAGACCTATTCGTCACGTTCAATGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa                                                                    <  3:518377/147‑1 (MQ=255)
                                                                                        aGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                              >  1:159745/1‑150 (MQ=255)
                                                                                            atTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCg                                                         >  3:121097/1‑151 (MQ=255)
                                                                                             tttATGGGAAAGGAGAAAATGCAGAGCTATTCGTCACTTTAAATGTTCAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcct                                                         <  3:695134/150‑2 (MQ=255)
                                                                                             tttATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTCAATGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTCGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                              <  4:60846/145‑1 (MQ=255)
                                                                                                 tgggaaaGGAGGAAATGCAGTGCTATTCCTCTCTTTAAATGGTAAAGCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                        <  2:1098551/142‑1 (MQ=255)
                                                                                                   ggaaagaagAAAATGCAGAGTTATTCGTCACTTTAAATGGTCAATCCCCCCCCTGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc                                                    <  1:63390/146‑1 (MQ=255)
                                                                                                     tAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                                >  3:938182/1‑151 (MQ=255)
                                                                                                      aaGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                                >  2:573997/1‑150 (MQ=255)
                                                                                                       aggaggAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTAt                                              >  1:354826/1‑151 (MQ=255)
                                                                                                        ggaggaAATGCAGAGCTATTCGTCACTTTAAATGTTCAAGCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTAt                                              <  2:1218090/150‑1 (MQ=255)
                                                                                                          aggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                           >  3:989381/1‑151 (MQ=255)
                                                                                                          aggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                           >  3:837003/1‑151 (MQ=255)
                                                                                                           gaaaaTGCAGAGCTATTCGCCCCGTTAAGTGGTCAACGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg                                            <  3:867396/147‑1 (MQ=255)
                                                                                                             aaaTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                        >  2:320113/1‑151 (MQ=255)
                                                                                                             aaaTGCAGAGCTATTCGTCACGTTAAATGGTCAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                           <  1:1897/148‑1 (MQ=255)
                                                                                                                   agagCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                    >  2:1132044/1‑149 (MQ=255)
                                                                                                                   agagCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                    >  3:1140947/1‑149 (MQ=255)
                                                                                                                      gCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc                               >  3:840441/1‑151 (MQ=255)
                                                                                                                       cTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc                               >  2:1362251/1‑150 (MQ=255)
                                                                                                                        ttttCGTCACTTTAAATGGTCAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc                               <  2:70007/147‑1 (MQ=255)
                                                                                                                        tATTCGCCACGTTAAATGGTCAACCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc                               <  4:104613/149‑1 (MQ=255)
                                                                                                                          ttCGTCACGTTAAATGGTCAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGAt                            <  2:387162/150‑1 (MQ=255)
                                                                                                                           tCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                          >  1:434873/1‑151 (MQ=255)
                                                                                                                           tCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                          >  2:721349/1‑151 (MQ=255)
                                                                                                                             gTCACGTCCAATGTTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                          <  3:948428/149‑1 (MQ=255)
                                                                                                                              tCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCgcgc                       >  2:167997/1‑151 (MQ=255)
                                                                                                                              tCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCgcgc                       >  2:165537/1‑151 (MQ=255)
                                                                                                                              tCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCgcgc                       >  3:163417/1‑151 (MQ=255)
                                                                                                                               cACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc                      >  2:1020843/1‑151 (MQ=255)
                                                                                                                               cACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc                      >  3:291193/1‑151 (MQ=255)
                                                                                                                                  gTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc                   >  4:43580/1‑151 (MQ=255)
                                                                                                                                    tCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcat                 >  3:865434/1‑151 (MQ=255)
                                                                                                                                       aTGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg               >  2:964609/1‑150 (MQ=255)
                                                                                                                                         ggTCAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattga              <  2:526557/149‑1 (MQ=255)
                                                                                                                                           tCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacatt          >  1:480343/1‑151 (MQ=255)
                                                                                                                                                   acacCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAc  >  2:1215997/1‑151 (MQ=255)
                                                                                                                                                   acacCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAc  >  2:648699/1‑151 (MQ=255)
                                                                                                                                                    cacCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAc  >  4:17203/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAC  >  NC_000913/2984292‑2984589

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