Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,4410AC47.1% 97.1 / 22.1 68G384G (GGT→GGGyqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (33/3);  new base C (0/32);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 2.60e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984292‑2984591
                                                                                                                                                     |                                                                                                                                                      
gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCacac                                                                                                                                                       >  2:77666/1‑151 (MQ=255)
gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCacac                                                                                                                                                       <  1:1383447/151‑1 (MQ=255)
   tatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCt                                                                                                                                                    >  2:1008864/1‑151 (MQ=255)
     tGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTgg                                                                                                                                                  >  1:782143/1‑151 (MQ=255)
          tGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATCTTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGccccc                                                                                                                                             >  2:626696/1‑151 (MQ=255)
             gaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgcc                                                                                                                                          <  2:80500/151‑1 (MQ=255)
               aTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgccg                                                                                                                                         >  1:355748/1‑150 (MQ=255)
                       aTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGc                                                                                                                                >  2:857965/1‑151 (MQ=255)
                         gggggAGCCCAAAAAATGTATGTCCACTTGCCAGCCATAATTATTCAACAGATTAAGGACGGCAGATATTTATGGGTAGGGAGAAAAGGAAGGGCTTTTCGCCACTTTAAATGGAAAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGt                                                                                                                              <  2:1420066/149‑1 (MQ=255)
                             gAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGcc                                                                                                                          >  2:941086/1‑151 (MQ=255)
                                     aaaaTGTATGTCGCCTTCCCAGCCATCATTATTCAACATATTAATGACGGCAGTTTTTTATGGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGCCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                                  <  2:713040/151‑1 (MQ=255)
                                     aaaaTGTATGTCGCCTGGCCACCAATCATTATAAAACAGATAAAGGACGGCAGATATTTATGGGTAAGGAGAAAATGAAGACCTTTTCGTCCCGTAAAATGGTCAACCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                                  <  1:384771/151‑1 (MQ=255)
                                     aaaaTGTATGTCGACTTCCCAGCCATCATTATTAAACAGATAAATGACGGAAGATATTTAGCGGAAGGGAGAAAATGCAGAGTTTTTCGTCACGTTAAAGGGTCAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                                  <  1:7324/151‑1 (MQ=255)
                                         tgtttgtCGATTCCCCAGCCATCTTTATAAAAAAGATTAATGCCGGCAGATTTTTATCGGTAAGGAGAAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                               <  2:897520/150‑1 (MQ=255)
                                           ttgttCGTCTCGCCACCCAATATTATCCAACATATTAATGACGGCAAATTTTTATCGGTAGGGGAGAAATGCAGAGCTTTTCGTCACGTACAGTGGAAAGTGCCCCACCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCggg                                                                                                             <  1:1374086/146‑1 (MQ=255)
                                                 gcttggccacccaATATTATCCAACAAATAAACGCCGGAAGATATATATCGAAAAGGAGAAAATGCAGAGCTTTTCGTCACGTAAAAGGGTCCATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                                                                         <  2:478831/143‑1 (MQ=255)
                                                   cTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                                                                                    >  1:582183/1‑151 (MQ=255)
                                                    ttGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt                                                                                                   >  1:856527/1‑151 (MQ=255)
                                                     tGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                                  >  2:446110/1‑151 (MQ=255)
                                                       ccagctaTCATTATTCAACAGATTAATGCCGGCAGTTATTTATCGGTAAGGAGAAAATGCAGAGCTATTCGTCCCGTACAATGTTAAAGGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                                  <  1:568335/149‑1 (MQ=255)
                                                              tcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                          >  2:586468/1‑150 (MQ=255)
                                                               catTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                        >  1:1408075/1‑151 (MQ=255)
                                                                attattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                       >  2:1440312/1‑151 (MQ=255)
                                                                   attCAACAGATTAATGACGGCAGATTTTTATGGGTAAGGAGAAAATGCAGAGCTTTTCGTCACTTTCAATGGTAAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGaa                                                                                     <  1:641514/150‑1 (MQ=255)
                                                                   attCAAAATATTAATGCCGGCAGATATTTCTCGGTAAGGAGGAAATGCAGAGCTTTTCGTCTCTTTCAATGGTAAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGaa                                                                                     <  2:186482/150‑1 (MQ=255)
                                                                    ttCAACAGATTAATGACGGCAGTTATTTATCGGTAAGGAGAAAATGCAGAGCTATTCGTCACTTTCAATGGTCAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                                                                         <  1:991506/145‑1 (MQ=255)
                                                                       aaaaGATTAATGACGGAAAATTTTTATCGGAAGGGAGAAAATGAAGAGCTATTCGTCACTTTAATGGGTAAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                          <  1:137917/138‑1 (MQ=255)
                                                                       aaCAGATTAATGCCGGCATATTTTTATCGTTAAGGAGGAAATGCAGAGCTTTTCGTCACTTTCATTGTTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                                 <  1:202384/150‑1 (MQ=255)
                                                                          aaaTTAATGACGGCAGATATTTACCGGTAAGAAGAAAATGCAGAGCTTTTCGTCCCGTTCAATGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAgaga                                                                              <  1:80630/148‑1 (MQ=255)
                                                                           gATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaa                                                                             >  2:1477176/1‑150 (MQ=255)
                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                                                                           >  1:1029514/1‑151 (MQ=255)
                                                                              tAATGACGGCAGATTTTTATGGGTAAGGAGGAAATGCAGAGTTTTTCGTCCCGTTCAATGGTAAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                          <  1:670681/150‑1 (MQ=255)
                                                                              tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                          >  1:423293/1‑150 (MQ=255)
                                                                               aaggACGGCAGATATTTACCGGTAAGGAGAAAAGGCAGAGCTTTTCTTCCCTTTAAAGGGTAAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAg                                                                        <  1:1495654/148‑1 (MQ=255)
                                                                                     ggaaaatatTTATCGGAAAGGAGGAAATGCAGTGTTTTTCGTCACGTTAATGGGTAATTCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt                                                                     <  2:984216/143‑1 (MQ=255)
                                                                                      gCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg                                                                 >  1:1205027/1‑151 (MQ=255)
                                                                                        aGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                                >  1:409483/1‑150 (MQ=255)
                                                                                          atatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAc                                                             >  2:669499/1‑151 (MQ=255)
                                                                                          atatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAc                                                             >  2:975556/1‑151 (MQ=255)
                                                                                           tatTTATGGGTAAGGAGAAAATGCGGTCCTTTTGTTCCCGTCCAATGGCAAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt                                                                     <  1:593408/142‑1 (MQ=255)
                                                                                           tatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                            >  2:441867/1‑151 (MQ=255)
                                                                                            atTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATTGTCAATGCCACACGCTGGCCCCCGCTGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                            >  2:669020/1‑150 (MQ=255)
                                                                                             tttATCGGTAGGGAGAAAAGGCAGATCTTTTCGTCACTTTCAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCg                                                           <  2:838450/150‑1 (MQ=255)
                                                                                             tttATCGGAAAGAAGAAAAGCCAGAGTTTTTCGCCACTTTAAATGGCCAACGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                                <  1:1135464/145‑1 (MQ=255)
                                                                                              ttATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                          >  1:243962/1‑150 (MQ=255)
                                                                                                 tCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc                                                      >  2:1465104/1‑151 (MQ=255)
                                                                                                  cGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                                     >  2:400822/1‑151 (MQ=255)
                                                                                                   ggaaaGGAGAAAATGCAGGGCTTTTCGTCCCGTCCAATGGTCAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc                                                    <  2:378914/148‑1 (MQ=255)
                                                                                                    gTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGcc                                                   >  1:585353/1‑151 (MQ=255)
                                                                                                    gTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGcc                                                   >  2:1092246/1‑151 (MQ=255)
                                                                                                    gTAAGGAGAAAAGGCAGAGTTATTCGTCTCGTTCATTGGTCAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                                     <  1:727930/149‑1 (MQ=255)
                                                                                                         gaggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg                                              >  2:1324648/1‑151 (MQ=255)
                                                                                                            gaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg                                            >  2:1432224/1‑150 (MQ=255)
                                                                                                             aaaGGCAGAGCTATTCGTCTGTTTAAATGGTAAACCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg                                            <  2:1303619/149‑1 (MQ=255)
                                                                                                                   agagCTATTCGTCACGTTCAATGGTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                      <  1:822718/149‑1 (MQ=255)
                                                                                                                    gagCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt                                   >  2:71620/1‑151 (MQ=255)
                                                                                                                         aTTCGCCCCGTTCAATGGTCAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc                                 <  2:1240040/148‑1 (MQ=255)
                                                                                                                          ttcttcACTTTCAATGGTCAACGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGAt                              <  2:251261/150‑1 (MQ=255)
                                                                                                                          ttCGTCACTTTAAATGGTCAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg                                  <  1:843441/146‑1 (MQ=255)
                                                                                                                          ttCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                             >  2:1036327/1‑151 (MQ=255)
                                                                                                                          ttCGCCACTTTAAATGGTAAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt                                   <  2:766405/145‑1 (MQ=255)
                                                                                                                            cGTCCCGTTCAATGGTCAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                            <  1:1524829/150‑1 (MQ=255)
                                                                                                                            cGTCACGTTAAAGGATCACGCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc                           <  1:1463382/151‑1 (MQ=255)
                                                                                                                                  gTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc                     >  2:769849/1‑151 (MQ=255)
                                                                                                                                    tCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcat                   >  1:1156069/1‑151 (MQ=255)
                                                                                                                                      aagggTCAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacca                    <  1:2724/145‑1 (MQ=255)
                                                                                                                                      aaTGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg                 >  2:1155642/1‑151 (MQ=255)
                                                                                                                                         gttaaaTGCCCCCCCCCGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAGCCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattga                <  1:206356/145‑1 (MQ=255)
                                                                                                                                         ggagaacgCCACACCCTGCCCGCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg                 <  2:1191262/141‑1 (MQ=255)
                                                                                                                                                 ccACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCg      <  1:109810/151‑1 (MQ=255)
                                                                                                                                                   acacCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAc    >  2:1453009/1‑151 (MQ=255)
                                                                                                                                                   acacCCTGGCCCCCGACGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAc    >  2:1323319/1‑151 (MQ=255)
                                                                                                                                                    cacCCCGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAc    <  2:465084/150‑1 (MQ=255)
                                                                                                                                                     acCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGc  >  2:1131176/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984292‑2984591

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