Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,4410AC51.4% 81.2 / 130.3 72G384G (GGT→GGGyqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (31/4);  new base C (0/37);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 2.35e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.17e-14

GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACC  >  NC_000913/2984292‑2984585
                                                                                                                                                       |                                                                                                                                                
gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCacac                                                                                                                                                 <  1:1034859/151‑1 (MQ=255)
gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCacac                                                                                                                                                 >  1:727009/1‑151 (MQ=255)
 aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAcc                                                                                                                                                >  1:882537/1‑151 (MQ=255)
  atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAccc                                                                                                                                               >  2:819680/1‑151 (MQ=255)
    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTg                                                                                                                                             <  1:1040390/151‑1 (MQ=255)
    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTg                                                                                                                                             <  1:899424/151‑1 (MQ=255)
    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTg                                                                                                                                             <  1:217276/151‑1 (MQ=255)
          tGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGccccc                                                                                                                                       >  2:560354/1‑151 (MQ=255)
            agaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgc                                                                                                                                     >  1:21918/1‑151 (MQ=255)
               aTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgccg                                                                                                                                   >  1:809898/1‑150 (MQ=255)
                     aaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACaa                                                                                                                             >  2:944756/1‑150 (MQ=255)
                      aaTGGGGGAGCCCAAAAAATTTTTGTGGACTTGCCACCAATAATTATTCAACAGATAAAGGCCGGCAGATATTTATCGGTAAGGAGAAAATGCAGGGCTTTTCGTCAGTTCAAATGGTCAATGCCCCCCCCTGGCCCCCGCCGATACAAAg                                                                                                                           <  1:196696/151‑1 (MQ=255)
                          aGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGtt                                                                                                                       >  2:888725/1‑151 (MQ=255)
                           gggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTg                                                                                                                      >  1:369966/1‑151 (MQ=255)
                           gggAGCCCAAAAAATTTATGCCGACTTGCCAGCAATCATTATTCAACAGATTAAGGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACTTTCAAGGTCCAATGCCACCCCCTGGCCCCCGCCGATACAAAGc                                                                                                                          <  1:1116418/147‑1 (MQ=255)
                           gggAGCCCAAAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGCCGGCAAATTTTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCCCTTTAAATGTTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGtt                                                                                                                       <  2:45955/150‑1 (MQ=255)
                                     aaaaTGTATGGCGACTGCCCCCCCATATTTATTAAACAGATTAATGACGTAAAATATTTCTCGGAAGGGAGAAAATGCAGACTTTTTCGTCCCTTTCAATGGTCACTCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                            <  1:586653/151‑1 (MQ=255)
                                       attgtatgtCGTCTGGCCACCCATTTTTATACAACAGATAAAGGACGGCAGATTTTTATCGGTAAGGAAAAAATGAAGATCTTTTCGTCACTTTAAATGGCCAACCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                                                          <  1:892628/149‑1 (MQ=255)
                                       aaTGTATGTCGCCTTGCCAGCCATCATTATTAAAAAGATTAAGGACGGCAGTTATTTATCGGTAAGGAGAAAATGCAGGGCTTTTCGTCACGTTAAATGGTCAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                                                          <  2:1030133/151‑1 (MQ=255)
                                        atgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                         >  1:74325/1‑151 (MQ=255)
                                        atgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                         >  2:552286/1‑151 (MQ=255)
                                             tgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGCCGGCAGATTTTTATCGGTAAGGAGAAAATGCAGAGTTTTTCTTCACTTTCAATGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCa                                                                                                     <  1:1102613/150‑1 (MQ=255)
                                             tggctaCTTCCCACCCATCATTTTTCAACAGATTAATGCCGGAAAATATTTATCGGTAAGGAGAAAAGGCAGATTTTTTCGTCACGTTAAATGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                                                          <  1:411551/140‑1 (MQ=255)
                                              gtcgtcTTGCCACCCATCATTATTAAACAGATAAAGGACGGCAAATATTTATGGGGAAGGAGAAAATGCAGAGCTTTTCGCCACGTTAAAGGGCCATTGCCCCACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCag                                                                                                    <  1:430701/150‑2 (MQ=255)
                                                cgccTTGCCAGCCATCATTATTAAAAAGATTAATGACGGCAGATATTTATGGGTAGGGAGAAAATGAAGGCTTATTCGCCCCTTTAAATGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTa                                                                                                  <  1:1104058/147‑1 (MQ=255)
                                                cGACTTGCCAGCAATCATTATTCAACAGATTAATGACGGCAAATTTTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACTTTAAATGTTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                                                                   <  1:686324/149‑1 (MQ=255)
                                                 gACTTGCCAGCCATCTTTATCAAACAGATAAAGGCCGCCAGTTATTTATGGGAAAGGAGGAAATGCAGAGCTATTCGCCACGTTAAAGGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc                                                                                                 <  1:422728/150‑1 (MQ=255)
                                                   cTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                                                                              >  1:378063/1‑151 (MQ=255)
                                                      cccacccaTCATTATCCAACAGATAAAGGACGGCGGATTTTTAGGGGAAAGGAGAAAAGGGGGATTTTTTCGTCACTTCCAAGGGTCATTGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                            <  2:89303/149‑1 (MQ=255)
                                                       ccggccaTCATTATTCAACAGATTAATGCGGGCAAATTTTTTTGGGAAAGGGGAAAATGCAGAGTTTTTCGCCACTTTCAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt                                                                                             <  2:409265/145‑1 (MQ=255)
                                                       ccggccaTCATTATTAAACAGATTAATGACGGCAGATTTTTTTGGGAAAGGAGAAAATGCAGAGTTTTTCGCCCCGTTAAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                            <  1:925233/146‑1 (MQ=255)
                                                       ccagccaTCATTATTAAAAAGATAAATGACGGAAGTTATTTACCGGTAAGGAGAAAATGCAGAGTTTTTCGTCACGTTAAATGGTCAATGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                            <  2:768367/149‑1 (MQ=255)
                                                          gccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAt                                                                                       >  1:696477/1‑151 (MQ=255)
                                                           caataattataCAACAAATTAATGACGGCAGATTTTTAGGGGAAGGGAGAAAATGCAGAGCTTTTCGTCCCGTTAAATGGTAAATCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCa                                                                                        <  1:373344/144‑1 (MQ=255)
                                                               catTATTCTACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                  >  2:865757/1‑151 (MQ=255)
                                                               catTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                  >  2:360834/1‑151 (MQ=255)
                                                                 ttttccAACGGATTAAGGCCGGAAGATATTTATCGGTAGGGGGGAAAGCAGGTGTTTTTCCTCAGGTAAAAGGGTCATTCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCa                                                                                        <  2:144455/138‑1 (MQ=255)
                                                                 ttattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                 >  1:459024/1‑150 (MQ=255)
                                                                 ttattCAACAGATTAATGACGGCAAATTTTTATCGGAAAGGAGAAAATGCAGAGTTTTTCGTCACTTTAAATGGTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                                                                   <  1:19934/148‑1 (MQ=255)
                                                                 ttattCAACAGATTAATGACGGAAGATATTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTCAATGGTAAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                  <  1:100537/149‑1 (MQ=255)
                                                                   attCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                              >  2:467627/1‑151 (MQ=255)
                                                                           gATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaa                                                                       >  1:477004/1‑150 (MQ=255)
                                                                           gATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaa                                                                       >  2:1010000/1‑150 (MQ=255)
                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                                                                     >  2:670351/1‑151 (MQ=255)
                                                                               aaggACGGCAGATATTTACCGGAAAGGAGGAAATGCAGAGCTTTTCGTCCCGTTCATTGGCCAACCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                   <  1:986325/147‑1 (MQ=255)
                                                                               aaTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                    >  2:756852/1‑149 (MQ=255)
                                                                                aTGACGGAAGATATTTAGCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGACAACGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                 <  2:791789/151‑1 (MQ=255)
                                                                                   aCGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGACCCAGAGAAACCAGGAtt                                                              >  2:790547/1‑151 (MQ=255)
                                                                                               tATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                   >  1:683231/1‑150 (MQ=255)
                                                                                                  cGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                               >  2:149045/1‑151 (MQ=255)
                                                                                                   ggTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc                                              >  2:677042/1‑151 (MQ=255)
                                                                                                   ggTAAGGAGAAAATGCAGAGTTTTTCGTCCCGTTAAAGGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                               <  2:1051260/150‑1 (MQ=255)
                                                                                                      aaGGAGAAATTGCAGACCTTTTCGTCACGTTAAATGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTa                                           <  2:497739/151‑1 (MQ=255)
                                                                                                       aggaggAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTAt                                          >  1:603663/1‑151 (MQ=255)
                                                                                                        ggaggaAATGCAGACCTTTTCGCCATGTTAAATGGTAAATCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                            <  2:537524/148‑1 (MQ=255)
                                                                                                          aggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                       >  2:121564/1‑151 (MQ=255)
                                                                                                             aaaTGCAGTGCTATTCGTGCCGTTAAGTGGTTAAACCCTCCCCCCGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                       <  2:503099/148‑1 (MQ=255)
                                                                                                             aaaTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                    >  1:657923/1‑151 (MQ=255)
                                                                                                             aaaTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                    >  1:1075473/1‑151 (MQ=255)
                                                                                                             aaaTGAAGAGCTTTTCGTCACGTTCAATGGCCAAGCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg                                        <  1:442353/147‑1 (MQ=255)
                                                                                                                     gagtttttCGCCCCGTTAAAGGGTCAACGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAcc                                 <  1:278639/141‑1 (MQ=255)
                                                                                                                     gagCTTTTCGTCTTGTTCATTGGTAAAGCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGa                               <  2:970359/149‑1 (MQ=255)
                                                                                                                      agCTATTCGTCCCGTTCATGGGTCAACCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg                            <  1:644719/151‑1 (MQ=255)
                                                                                                                        cTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg                          >  1:500418/1‑151 (MQ=255)
                                                                                                                         ttttCGTCACTTTAAATGGTCAAGGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg                            <  2:1008956/146‑1 (MQ=255)
                                                                                                                           tttgtcggtttCAACGGTAATTGCCTCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt                             <  2:704959/136‑1 (MQ=255)
                                                                                                                             cGTCACTTTAAAGGGTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                       <  1:819639/149‑1 (MQ=255)
                                                                                                                              gTCACTTTAAATGGTCAAGCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc                     <  2:962623/150‑1 (MQ=255)
                                                                                                                              gTCACTTTAAAGGGTCAACGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc                     <  1:750184/150‑1 (MQ=255)
                                                                                                                               tctcGTTAAATGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgcg                    <  2:1046196/147‑1 (MQ=255)
                                                                                                                               tCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCgcgc                   >  2:524922/1‑151 (MQ=255)
                                                                                                                                               aaTGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacatt      <  2:432239/149‑1 (MQ=255)
                                                                                                                                               aaTGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg     >  1:941556/1‑150 (MQ=255)
                                                                                                                                                   ccACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGAcc  <  1:986041/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                
GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACC  >  NC_000913/2984292‑2984585

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