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

GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAC  >  NC_000913/2984292‑2984589
                                                                                                                                                       |                                                                                                                                                    
gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAACGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCacac                                                                                                                                                     <  2:224382/151‑1 (MQ=255)
  atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAccc                                                                                                                                                   >  2:662124/1‑151 (MQ=255)
    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTg                                                                                                                                                 >  1:1005158/1‑151 (MQ=255)
    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTg                                                                                                                                                 <  1:607691/151‑1 (MQ=255)
    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTg                                                                                                                                                 >  1:364614/1‑151 (MQ=255)
      gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGc                                                                                                                                               <  1:1069358/151‑1 (MQ=255)
           gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccg                                                                                                                                          <  1:410810/151‑1 (MQ=255)
                  gAAACATGAGGGAGCCCAGAAAATGTATGTCTACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGTAGGAAATGCAGTGCTATTCGTCTCGTTCAATTGTCTATGCCTCACCCTGGCCCccg                                                                                                                                          <  1:1022384/144‑1 (MQ=255)
                         gAGGGAGCCAAGAAAATGTATGTCGACTTGCCAGCCATAATTATTCAACAGATTAATGACGGCAGATATTTATGGGTAAGGAGAAAATGCAGAGCTTTTCGTCACTTTAAATGGTCAATCCCACCCCCTGGCCCCCGCCGATACAAAGc                                                                                                                              <  2:99819/149‑1 (MQ=255)
                             gAGCCCAAAAAATGTATGTCGCTTGCCCGGCCATCATTATTCAACAGATTAATGCCGGCAGATATTTATCGGTAAGGAGAAAAGGAGGAGTTTTTCTTCACTTTCAAGGGTAAAGCCCCCACCCGGGCCCCCGCCGATACAAAGCGTTGcc                                                                                                                        <  2:692043/151‑1 (MQ=255)
                                  cAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTcc                                                                                                                   >  2:75165/1‑151 (MQ=255)
                                   aaaaaaTGTATGTCGACTTGCCACCCATCATTATACAACAGATTAAGGACGGCAGATTTTTATGGGTAAGGAGAAAATGCAGAGCTATTCGTCCCGTTCAATGGTAAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTc                                                                                                                    <  1:389364/147‑1 (MQ=255)
                                   aaaaaaTGTATGTCGACTGGCCAGCCATCATTATCAAACAGATAAAGGAGGGAAATTATTTAGCGGGAGGGAGAAAATGAAGAGCTTTTCGTCACGTCCAATGGTAAAGGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAg                                                                                                                      <  1:61610/145‑1 (MQ=255)
                                   aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                                  >  2:39391/1‑151 (MQ=255)
                                     aaatTGTATGTCGACTTGCCAGCCATAATTATTAAAAAGATAAATGACGGAAGATATTTATGGGAAGGGAGAAAATGCAGAGCTTTTCTCCACGTACAATGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                                  <  2:1028793/149‑1 (MQ=255)
                                      aaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                                                               >  1:850335/1‑151 (MQ=255)
                                      aaaTGTATGCCGACTTGCCAGCCTTCATTATTCAAAAAATTAATGAGGGAAGATTTTTAGGGGTAAGGAGAAAAGGAAGAGCTTTTCGTCACTTTAAATGGCAAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                                                               <  1:1136285/151‑1 (MQ=255)
                                           tatgtCGACTTGCCAGCCATCATTATTCAAAAGATTAATGACGGCAGATTTTTAGCGGAAAGGAGAAAAGGCAGACCTTTTCGTCACGTTAAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                             <  1:270397/148‑1 (MQ=255)
                                              gtCGACTTGCCACCCATCTTTATTAAACAGATTAATGACGGCAGATATTTATCGGAAAGGAGGAAATGCAGAGCTTTTCGTCACTTTAAAGGGTCAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAt                                                                                                        <  2:766226/150‑1 (MQ=255)
                                               tCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                                                                       >  1:1048382/1‑150 (MQ=255)
                                               tCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                                                                       >  2:1120557/1‑150 (MQ=255)
                                                cGACTTGCCATCAAGAATTATACAACAAATAAATGGCGGAAAATATTTAGAGGGAGGGAAAAAATCCAGATGTATTCGTCATTTTCAATGTTCACCCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc                                                                                                     <  2:62285/151‑1 (MQ=255)
                                                   atttgccgccaatcatTATTAAAAAAATTAATGACGGCAAATATTTAGCGGAAAGGAGAAAAGGAAGGGTTTTTGGTCACGTTAAAGGGTAAAGGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGc                                                                                                          <  2:90998/133‑1 (MQ=255)
                                                    tggccagcaaATATTATTAAACAGATAAATGCCGGAAAATATTTAGCGGAAAGGAGAAAATGAAGGCCTTTTCGTCACGTTAAAGGGTAAACCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg                                                                                                            <  1:130564/138‑1 (MQ=255)
                                                     tGCCAGCCATCATTATTCAACAGATTAATGACGGCAGTTTTTTATCGGAAAGGAGAAAATGCAGAGCTTTTCGTCCCTTTCAATGGCCAATCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                                                                                  <  1:850485/149‑1 (MQ=255)
                                                     tGCCAGCCATCATTATTCAACAGATTAAGGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTAAATGGTCAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                                                                                  <  2:225655/149‑1 (MQ=255)
                                                       ccaagccaTATTTATACAACAAATTAATGACGAAAGATATTTAGCGGAAAGGAGGAAATGAAGATCTTTTCGTCCCGTAAAATGGTAAACCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa                                                                                               <  1:426346/148‑1 (MQ=255)
                                                         cagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                             >  1:694015/1‑151 (MQ=255)
                                                            ccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt                                                                                          >  2:656074/1‑151 (MQ=255)
                                                              atcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCCCACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt                                                                                         <  1:121940/150‑1 (MQ=255)
                                                                catTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                      >  2:194365/1‑151 (MQ=255)
                                                                catTATCCAACAGATGATTGTCGAAATATTTTTATCGGAAAGGAGGATTTTCAGGCTTCTTCATAACTTCCAATGTTAATTACCACCCCCGGGCCACCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAt                                                                                           <  2:509830/146‑1 (MQ=255)
                                                                 attattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                     >  1:1093310/1‑151 (MQ=255)
                                                                  ttattaAAAAGATTAATGACGGCAATTTTTTATGGGAAAGGAGAAAATGCAGAGTTTTTCGTCACGTTAAAGGGTCAATCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                      <  1:1034157/149‑1 (MQ=255)
                                                                    attCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                                  >  2:610223/1‑151 (MQ=255)
                                                                     ttCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                                 >  1:64792/1‑151 (MQ=255)
                                                                       cAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                                 >  2:175571/1‑149 (MQ=255)
                                                                        aaCAGTTTAAGGACGGCAGATTTTTATCGGTAAGGAGAAAAGGCAGAGTTTTTCTTCCCGTTCAAGGTTAAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                                                                                <  2:952554/149‑1 (MQ=255)
                                                                         aCATATTATTGACGGCAAATATTTAGCGGAAAGGAGAAAATGCAGAGTTTTTCTTCACTTAAAGTGGACAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCagag                                                                             <  1:563402/151‑1 (MQ=255)
                                                                                     cGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc                                                                 >  2:383910/1‑151 (MQ=255)
                                                                                      ggaaGATATTTATCGGTAAGAAGAAAATGCAGAGTTATTCGCCACTTTCAATGTTAAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa                                                                    <  2:355448/144‑1 (MQ=255)
                                                                                         aGATATTTATCGGTAGGGAGGAAATGCAGAGTTTTTCGTACCTTTCAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg                                                               <  1:655521/149‑1 (MQ=255)
                                                                                              tttAGCGGAAGGAAGAAAATGGAGAGCTTTTCCTCTTTTCAACTGCTCAACCCCACACCCCGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCa                                                             <  2:1099083/146‑1 (MQ=255)
                                                                                                      tAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                                >  1:1085378/1‑151 (MQ=255)
                                                                                                       aaGGAGAAAATGCAGAGTTTTTCGTCACTTTAAATGGCAAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                                <  2:12389/150‑1 (MQ=255)
                                                                                                          gagaaAATGCGGATCTTTTCGTAACTTTAAATGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                                <  1:1031736/147‑1 (MQ=255)
                                                                                                            ggaAGTGCAGACCTTTTCGTCATTTACAATGGTAAGTGCCGCCCCGTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                           <  1:723120/150‑1 (MQ=255)
                                                                                                                aTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAc                                      >  2:576522/1‑151 (MQ=255)
                                                                                                                  gCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                    >  2:115202/1‑151 (MQ=255)
                                                                                                                    agagCTTTTCGTCACGTTAAATGGTAAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                    <  2:121709/149‑1 (MQ=255)
                                                                                                                     gagCTTTTCGTCACGTACAATGGTAAACCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAcc                                     <  1:348415/147‑1 (MQ=255)
                                                                                                                       gCTTTTCGTCACGTTCAAGGGTAAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc                               <  1:1104339/151‑1 (MQ=255)
                                                                                                                       gCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc                               >  2:1087843/1‑151 (MQ=255)
                                                                                                                            tCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                          >  1:510828/1‑151 (MQ=255)
                                                                                                                              gTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgcg                        >  2:1139525/1‑151 (MQ=255)
                                                                                                                                cACGTTCAATGGCCAATCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc                         <  2:243971/148‑1 (MQ=255)
                                                                                                                                 aCGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc                      >  1:37328/1‑150 (MQ=255)
                                                                                                                                  acgtaaaaggtgaaaTGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCgcgc                       <  2:1075226/136‑1 (MQ=255)
                                                                                                                                    ttCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc                   >  1:566286/1‑150 (MQ=255)
                                                                                                                                     tCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcat                 >  1:604885/1‑151 (MQ=255)
                                                                                                                                      cAATGGTCAATGCCACACCCTGTCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcatt                >  1:784529/1‑151 (MQ=255)
                                                                                                                                            gtcatgcccACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattga              <  2:1121113/141‑1 (MQ=255)
                                                                                                                                               aaTGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg         <  2:136360/150‑1 (MQ=255)
                                                                                                                                               aaTGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg         >  1:280054/1‑150 (MQ=255)
                                                                                                                                                     acacCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAc  >  2:184896/1‑151 (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 < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: