Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,4470GC44.4% 115.8 / 25.8 72G382G (GGC→GGGyqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (31/9);  new base C (0/32);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 1.25e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCATCC  >  NC_000913/2984297‑2984597
                                                                                                                                                      |                                                                                                                                                      
tGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTgg                                                                                                                                                        >  2:776659/1‑151 (MQ=255)
 gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTgg                                                                                                                                                        <  1:658160/150‑1 (MQ=255)
 gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGc                                                                                                                                                       >  1:286384/1‑151 (MQ=255)
      gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccg                                                                                                                                                  >  1:1212794/1‑151 (MQ=255)
       agaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgc                                                                                                                                                 <  2:878679/151‑1 (MQ=255)
        gaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgcc                                                                                                                                                <  1:786924/151‑1 (MQ=255)
         aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgccg                                                                                                                                               >  2:764971/1‑151 (MQ=255)
              aaaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATAc                                                                                                                                           >  2:191507/1‑150 (MQ=255)
                    gAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGt                                                                                                                                    >  2:1114721/1‑151 (MQ=255)
                      gggAGCCCAAAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATGGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGTAAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTg                                                                                                                                  <  1:693921/151‑1 (MQ=255)
                            ccAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTc                                                                                                                            >  1:959645/1‑151 (MQ=255)
                              aGAAAATGTATGTCGCCTTGCCACCCATCATTATTCAACAGATTAATGCCGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTc                                                                                                                            <  2:770266/149‑1 (MQ=255)
                                aaaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                                        >  2:951703/1‑151 (MQ=255)
                                  aaTGTATGTCGACTTGCCAGCCATCATTATCAAACAAATTAAGGCCGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTCAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                                                                      <  1:471986/151‑1 (MQ=255)
                                    tgtatgtCGCCTTGCCCGCCATCATTATCAAAAAGATTAAGGCCGGCAGATTTTTATGGATAAGGAGAAAAGGCAGATCTTTTCTCCCCTTTAAAGGGACAACCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                                     <  1:1266769/150‑1 (MQ=255)
                                      tatgtCGTCTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGAAGAGCTTTTCGTCACTTTAAATGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg                                                                                                                    <  2:353026/149‑1 (MQ=255)
                                      tatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTATTCGTCACTTTAAATGGTCAATCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg                                                                                                                    <  1:179675/149‑1 (MQ=255)
                                               ttcccacccaTCTTTATTAAAAAAATTAATGCCGGAAAATTTTAAGGGGTAGGGAGAAAAGGCAGAGTTTTTCGTCCCTTTAAATAGTCAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                                                                                          <  1:435004/147‑1 (MQ=255)
                                               ttcccacccaTATTTATTAAACAGATTAATGCCGGAATATATTTATCGGTAAGGAGAAAATGCAGAGTTTTTCGTCCCTTTAAAGGGTAAGCCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                                            <  1:558957/145‑1 (MQ=255)
                                               ttGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGTTTTTCGTCACGTTAAATGTTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                                            <  1:666555/148‑1 (MQ=255)
                                               ttGCCAGCCATAATTATTAAACAGATAAATGACGAGATATTTTTATGGGAAAGGAGAAAATGCAGATCTTTTCGTCACTTTAAAGGGTAAAGCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                                            <  1:463874/148‑1 (MQ=255)
                                                 cccacccaTCATTATTAAAAAGATAAATGACGGCAGATTTTTATCGGAAAGGAGAAAATGCAGAGTTTTTCGTCCCGTTCAAGGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                                                                                          <  1:951693/147‑1 (MQ=255)
                                                  ccaccctTCTTTATTCAACAGTTAAATGCCGGAAAATTTTTATGGGAAAGGAGAAAATGCAGATCTTTTCGTCCCGTTAATTGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt                                                                                                         <  2:479127/148‑1 (MQ=255)
                                                  ccaaccaTAATTATTCAACAGATTAATGACGGAATATATTTATCGGAAAGGAGGAAAGGAAGGGTTTTTCGTCACGTTAAATGGTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa                                                                                                       <  1:1010087/150‑1 (MQ=255)
                                                   cagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                                     >  2:1288431/1‑151 (MQ=255)
                                                    agccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                                     >  1:222368/1‑150 (MQ=255)
                                                    agccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                                     >  1:1066732/1‑150 (MQ=255)
                                                      ccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt                                                                                                  >  2:836303/1‑151 (MQ=255)
                                                         tcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                                >  2:1305959/1‑150 (MQ=255)
                                                            ttattCAACAGATAAATGCCGGAAGATTTTTATCGGAAAGGAGAAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt                                                                                                 <  1:302022/146‑1 (MQ=255)
                                                             ttttCAACAGTTTAATGGCGGCAGATATTTATGGGAAAGGAGAAATGGCAGAGCTTTTCGTCACTTTCAAGGGTAAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                              <  2:837116/146‑1 (MQ=255)
                                                             tattCAACAGATTAATGACGGAAGATATTTATCGGTAAGGAGAAAATGCAGAGTTTTTCGTCACTTTAAATGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                                                                               <  1:1277685/147‑1 (MQ=255)
                                                              attCAACAGATTAATGACGGCAGATATTTATCGTTAAGGAGAAAATGCAGAGCTTTTCGCCACGTTCAATGGTCAATGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                                          <  2:1055152/151‑1 (MQ=255)
                                                                  aaaaaaTTAATGACGGCAAATTTTTATCGGAAAGGAGAAAATGCAGAGTTATTCGTCACGTTAAGGGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                                                                                        <  1:292141/144‑1 (MQ=255)
                                                                  aaCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                                         >  2:202976/1‑148 (MQ=255)
                                                                  aaCAGATTAAGGCCGGCAGATATTTATCGGTAAGGAGAAAATGCAGGGTTATTCGTCACTTTAAATGGTCAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                                          <  2:349342/147‑1 (MQ=255)
                                                                          aaggACGGAAGATTTTTAGCGGTAAGGAGGAAAGGCAGAGCTATTCGTCACGTTAAATGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                               <  1:1127808/147‑1 (MQ=255)
                                                                          aaTGCCGGCAGATTTTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCCCTTTAAATGGTCAATGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                               <  2:1279641/150‑1 (MQ=255)
                                                                          aaTGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGTTATTCGTCCCTTTAAATGTTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                               <  2:281051/150‑1 (MQ=255)
                                                                           aTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                             >  1:575150/1‑151 (MQ=255)
                                                                            tGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                             >  1:200809/1‑150 (MQ=255)
                                                                              aCGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAtt                                                                          >  2:1167820/1‑151 (MQ=255)
                                                                                   aGATATTTATCGGAAAGGAGAAAATGCAGAGCTTTTCTCCACTTTAAAGGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg                                                                       <  1:1219583/149‑1 (MQ=255)
                                                                                      tatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                                  >  2:127105/1‑151 (MQ=255)
                                                                                      tatTTATCGGAAAGGAGAAAATGAAGAGTTTTTCTTCACTTTAAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCa                                                                     <  1:937664/148‑1 (MQ=255)
                                                                                         ttgtgggTAAGGGGGAAAGGCAGAGCTTTTCGCCACGTTAAAGGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                                <  1:295661/145‑1 (MQ=255)
                                                                                           accGGAAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAAGGGTCAACCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGaa                                                              <  1:459429/148‑1 (MQ=255)
                                                                                           aTCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                               >  2:314711/1‑149 (MQ=255)
                                                                                             cGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                                           >  1:852085/1‑151 (MQ=255)
                                                                                                   ggaggaAATGAAGAGCTTTTCGTCACGTAAAAGGATCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc                                                            <  1:1003481/144‑1 (MQ=255)
                                                                                                     agaaaaTGCAGGGCTTTTCGCCACGTTAAATGGTCAAGCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATc                                                     <  1:939508/146‑1 (MQ=255)
                                                                                                      ggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg                                                  >  1:365931/1‑151 (MQ=255)
                                                                                                         aaTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                                >  1:969042/1‑150 (MQ=255)
                                                                                                           tGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAcc                                             >  2:825054/1‑151 (MQ=255)
                                                                                                                agCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg                                        >  2:232057/1‑151 (MQ=255)
                                                                                                                  cTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg                                      >  1:1220483/1‑151 (MQ=255)
                                                                                                                     ttcgccagtttaaaTGGTAAAGCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAcc                                             <  1:67466/129‑1 (MQ=255)
                                                                                                                     ttCGTCACGTTAAAGGGTAAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGa                                     <  2:652393/149‑1 (MQ=255)
                                                                                                                        gTCACTTTAAATGGTCAATGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                                   <  2:418971/148‑1 (MQ=255)
                                                                                                                           aCGTTAAAGGGACAAGGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                                   <  2:222371/145‑1 (MQ=255)
                                                                                                                             tttCAAGGGTCAATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcca                             <  1:660069/148‑1 (MQ=255)
                                                                                                                              ttaaaTGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCgcgc                               <  1:397120/143‑1 (MQ=255)
                                                                                                                              ttCAATGGTAAATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGctggtttctctgg                           <  1:202973/150‑13 (MQ=255)
                                                                                                                               tCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcat                         >  2:502161/1‑151 (MQ=255)
                                                                                                                                     gTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacat                   >  2:359697/1‑151 (MQ=255)
                                                                                                                                        aatcccccccccTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattga                <  2:946420/142‑1 (MQ=255)
                                                                                                                                        aaTGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg                 <  2:1211418/150‑1 (MQ=255)
                                                                                                                                        aaTGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg                 >  1:304533/1‑150 (MQ=255)
                                                                                                                                        aaTGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg                 >  2:793853/1‑150 (MQ=255)
                                                                                                                                            ccACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCg            <  1:86159/151‑1 (MQ=255)
                                                                                                                                             caccCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGa           <  1:610032/151‑1 (MQ=255)
                                                                                                                                              acacCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAc          >  1:972756/1‑151 (MQ=255)
                                                                                                                                                  ccTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTc      >  1:440191/1‑151 (MQ=255)
                                                                                                                                                   cTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTc      >  1:497864/1‑150 (MQ=255)
                                                                                                                                                      gCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCATcc  >  1:62217/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
TGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCATCC  >  NC_000913/2984297‑2984597

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: