Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,4410AC40.5% 128.1 / 18.2 75G384G (GGT→GGGyqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (37/7);  new base C (0/30);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 4.39e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984291‑2984591
                                                                                                                                                      |                                                                                                                                                       
cGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCcaca                                                                                                                                                         <  1:888030/151‑1 (MQ=255)
  aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAcc                                                                                                                                                       >  2:591789/1‑151 (MQ=255)
   atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAccc                                                                                                                                                      >  2:537196/1‑151 (MQ=255)
   atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAccc                                                                                                                                                      >  2:638401/1‑151 (MQ=255)
         cATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                                                                                <  1:391314/150‑1 (MQ=255)
         cATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                                                                                <  2:1007111/150‑1 (MQ=255)
            gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGccccc                                                                                                                                             <  1:823317/150‑1 (MQ=255)
             agaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgc                                                                                                                                           >  1:285820/1‑151 (MQ=255)
               aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgccg                                                                                                                                         >  1:548983/1‑151 (MQ=255)
                aTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGa                                                                                                                                        <  2:283903/151‑1 (MQ=255)
                      aaaTGAGGGAGCCAAAAAAATGTATGCCGACTTGCCAGCAATCATTATACAACAGATTAATGACGGCAGATATTTAGCGGAAGGGAGGAAATGCAGAGCTATTCGCCACTTAAAATGGTCAATCCCCCCCCCGGCCCCCCGCCGATACaaa                                                                                                                                  <  2:1070258/151‑1 (MQ=255)
                           tgggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGtt                                                                                                                             >  2:710171/2‑151 (MQ=255)
                           aGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGtt                                                                                                                             >  2:519965/1‑151 (MQ=255)
                            ggggcccaaaaaaaTTTATGTCGTCTTGCCAGCCATTATTATTAAACATATTAATGACGACAAATTTTTACCGGTAAGGAGAAAATGAAGAGTTTTTCGCCACTTTAAATGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGtt                                                                                                                             <  2:807667/140‑1 (MQ=255)
                                 ccaaaaaaaTTTATGTCGCCTTCCCACCCATCATTATTCAACAGATTAATGCCGGCAGATTTTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGCCACTTTCAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCa                                                                                                                         <  1:415621/144‑1 (MQ=255)
                                 ccaaaaaaaTGTATGTCGATTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTTTTCGTCACTTTCAATGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCa                                                                                                                         <  1:378734/144‑1 (MQ=255)
                                 caaaGAAATTGTATGTCGCCTTCCCGGCCATTATTATACAAAAAATAAAGGACGGAAGATTTTTATGGGAAAGGAGGAAATGCAGACTTTTTCGTCAGTTTAAAGGGTAAACCCCAACCCCGGCCCCCCGCCGATACAAAGCGTTGCCa                                                                                                                         <  1:941536/146‑1 (MQ=255)
                                    aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                                    >  2:714346/1‑151 (MQ=255)
                                    aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATAAATGACGGCAGATTTTTATCGGTAAGGAGAAAATGCAGAGTTTTTCGTCACTTTAAATGGTCAACCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                                    <  1:403093/151‑1 (MQ=255)
                                      aaaaTTTTTGTCGTTTTCCCAGCCATAATTATTAAACAAATTAATGACGGCATATTTTTAGCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGTCATTGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                                  <  2:434505/151‑1 (MQ=255)
                                      aaaaTGTTTGTCGTTTAGCCAACCATTATTATTAAACAGATTAAGGCCGGCAGATTTTAATCGGAAAGGAGAAAAGGCAGTGCTTTTCGTCACTTTCAGTGGTAAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                                  <  2:540495/151‑1 (MQ=255)
                                      aaaaTGTTTGTCGACTGGCCAGCCATAATTATCAAACAGATTAATGACGGCAGATTTTTATCGGTAAGGAGGAAAGGCAGAGCTTTTCGTCACGTTAAGTGGTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCtt                                                                                                                   <  2:537773/150‑1 (MQ=255)
                                      aaaaTGTATTTCGACTTGCCAGCCATAATTATTAAAAAATTAAATGGCGGCAGATTTTTATGGGAAAGAAGGAAGTGAGGAGTTATTCGAAACTTTAATTGGTCAATGCACCACCGTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                                  <  2:1014867/151‑1 (MQ=255)
                                      aaaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGAAGATATTTATCGGAAAGGAGAAAATGCAGAGCTTTTCGTCACGTTCAAGGGTCAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                                  <  1:161047/151‑1 (MQ=255)
                                       aaagttatgtCCCCTTGCCAGCAATCATTATAAAACAGTTTAACGGAGGTACATATTTCGCGGAAGGGAGAAATTGAGGAGCTTTTCTCCACGTAAAATGGTAAGCCCCACCCCCCGTCACCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                                                                 <  1:1023364/146‑1 (MQ=255)
                                          tgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg                                                                                                              >  1:677579/1‑151 (MQ=255)
                                          tgtatgtCGACTTGCCAGCCATAATTATTAAAAAGATAAATGACGGCAGATTTTTATGGGTAAGGAGAAAATGCAGAGTTTTTCGTCACTTTAAAGGGTCAATCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                               <  1:80743/150‑1 (MQ=255)
                                              tgccgccTTGCCAGCCATCTTTATTAAAAAGATAAATGCCGGCAGATATTTATGGGAAAGGAGAAAAGGCAGAGTTATTCGTCACTTTCAATGGTCAAGCCCCCACCTTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                                                                <  1:337756/139‑1 (MQ=255)
                                               gtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                                                                         >  2:100679/1‑151 (MQ=255)
                                                     ttGCCGGCCATCATTATTTAAGAGAGAAATGACGGCGGATATTTATCGGTAAGGAGGATATGCAGAGCTATTAATAAAGTTCAATGGTCAATGCCACACCCTGGCGCCCGCAGATACAAAGCGTTGCCAGTACTTTGCGGGCATTACGttt                                                                                                   >  2:235628/1‑151 (MQ=255)
                                                      tGCCAGCCATATTTATTCAAAAGATTAATGACGGCAGATATTTATCGGAAGGGAGAAAATGCAGAGCTTTTCGTCACGTTCAATGTTCAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGt                                                                                                     <  1:137555/148‑1 (MQ=255)
                                                       gacagccaTCATTATCAAAAAGATTAAAGCCGGCAGATATTCAGCGGTAAGAAGGAAAGGCGGGCATATTCGTCGCGTTAAATTGCCAATCCCACTCCGTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTa                                                                                                        <  2:425144/142‑1 (MQ=255)
                                                        ccagccaTCATTATTCAACAGATTAATGACGGAAGATATTTATCGGTAAGGAGGAAATGCAGAGTTTTTCGTCACGTTCAATGTTAAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAc                                                                                                <  1:80340/151‑1 (MQ=255)
                                                        caaacaatatttattCAACAAATTAAGGACGAAAGATATTTATGGGTAAGGAGGAAATGCAGAGCTATTCCTCACGTTCAATGGTAACCCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                                  <  1:14316/139‑1 (MQ=255)
                                                          agcaatattttttCAACATATTAATGACGGAATATTTTTATCGGTAAGGAGAAAAGGCAGTCTTTTTCCCTCCGTTAATTGGTCACTGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTTCTTTGCGGGCATTACGTTTCa                                                                                                 <  2:181846/137‑1 (MQ=255)
                                                             catcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGCTACCAGTCCTTTGCGGGCATTACGTTTCACCAt                                                                                             >  2:1296437/1‑149 (MQ=255)
                                                               tcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                          >  1:1150278/1‑150 (MQ=255)
                                                                    attCAAAAGATAAATGACGGAAGATATTTATCGGAAAGGAGGAAATGCAGAGCTTTTCGTCACGTTAAATGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa                                                                                      <  1:881756/149‑1 (MQ=255)
                                                                     ttCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAct                                                                                                      >  2:448921/1‑131 (MQ=255)
                                                                        aaaaGATTAATGACGGAAGATTTTTATCGGAAGGGAGAAAATGCAGAGCTTTTCGTCACGTTCAATGGTAAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAct                                                                                                      <  1:448722/126‑2 (MQ=255)
                                                                         aCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                                 >  2:995356/1‑149 (MQ=255)
                                                                                 aTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGTCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                       >  2:737654/1‑151 (MQ=255)
                                                                                    aCGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAtt                                                                    >  2:175928/1‑151 (MQ=255)
                                                                                      ggCAGTTATTTATGGGAAAGGAGAAAATAAAGACCTATTCGTCCCTTCCAAGGGTAATTGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                                  <  1:35633/151‑1 (MQ=255)
                                                                                      gaaatatatTTATCGGAAAGGAGGAAATGCAGAGCTTTTCTTCACGTTAAAGGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa                                                                      <  2:1182908/142‑1 (MQ=255)
                                                                                              tttATCGGAAAGGAGAAAATGCAGAGCTTTTCGTCACTTTCAATGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCa                                                               <  1:882182/146‑1 (MQ=255)
                                                                                               gtATCGGAAAGGAGGAAATGCAGAGAAAAAAGGCACGTCCACTTGTCAATGCCAGAGGCTGGCACCCACCGATACAAAACCTTGACAGTCCTTTGCGGGCATTAAGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                            >  2:818051/2‑148 (MQ=255)
                                                                                                tagcGGAAAGGAGAAAATGCAGAGTTTTTCGTCACGTTCATGGGTAAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                         <  2:443634/147‑1 (MQ=255)
                                                                                                      tAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                                  >  2:1266102/1‑151 (MQ=255)
                                                                                                        aggaggAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTAt                                                >  1:148827/1‑151 (MQ=255)
                                                                                                         ggaggaAATGCAGAGCTTTTCGTCGCGTTAAAGGCTAAAGCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTAt                                                <  1:1085734/150‑1 (MQ=255)
                                                                                                           aggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                             >  2:480095/1‑151 (MQ=255)
                                                                                                            gaaaaTGCAGAGCTATTCGTCACTTTAAATGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                             <  1:1329787/148‑1 (MQ=255)
                                                                                                             gaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGt                                           >  2:24090/1‑151 (MQ=255)
                                                                                                               aaTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                          >  1:879118/1‑150 (MQ=255)
                                                                                                                    agagCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAcc                                       >  1:1311258/1‑148 (MQ=255)
                                                                                                                    agagCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                      >  1:678834/1‑149 (MQ=255)
                                                                                                                    agagCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                      >  1:1267247/1‑149 (MQ=255)
                                                                                                                           ttCGTCGCGTTCAATGGTCAAGCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGAt                              <  2:753733/150‑1 (MQ=255)
                                                                                                                           ttCGCCCCGTTCAAGGGTAAACCCCCCACCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt                                   <  1:698975/145‑1 (MQ=255)
                                                                                                                            tcgtcgCGTTCAGTGGTCAATGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGa                               <  2:77075/148‑1 (MQ=255)
                                                                                                                            tCGTCATTTACAATGTCCAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                             <  2:583950/150‑1 (MQ=255)
                                                                                                                            tCGTCAGGTTAAATGTTAATGCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                             <  2:1244837/150‑1 (MQ=255)
                                                                                                                            tCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                            >  1:1166328/1‑151 (MQ=255)
                                                                                                                            tCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                            >  1:157130/1‑151 (MQ=255)
                                                                                                                            tCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                             >  2:887582/1‑150 (MQ=255)
                                                                                                                            tCGTCACGTTCAATGGTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGa                               <  1:118281/148‑1 (MQ=255)
                                                                                                                                cACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc                        >  1:506685/1‑151 (MQ=255)
                                                                                                                                 ccGTTCAATGGTCAAGCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc                           <  2:353403/146‑1 (MQ=255)
                                                                                                                                    ttCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc                     >  1:406939/1‑150 (MQ=255)
                                                                                                                                      cAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcatt                  >  1:601051/1‑151 (MQ=255)
                                                                                                                                       aaTGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg                 >  1:565108/1‑151 (MQ=255)
                                                                                                                                        aTGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg                 >  2:1324443/1‑150 (MQ=255)
                                                                                                                                        aTGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgac               >  1:831095/1‑151 (MQ=255)
                                                                                                                                              aaTGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg           >  2:1080454/1‑150 (MQ=255)
                                                                                                                                              aaTGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacat             >  2:972803/1‑148 (MQ=255)
                                                                                                                                                  ccccacCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCg      <  1:1105026/148‑1 (MQ=255)
                                                                                                                                                     cacCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGa     <  1:791414/149‑1 (MQ=255)
                                                                                                                                                      acCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTTCGATCGCGCCACCATTGACATTGACCCGACGc  >  1:478061/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
CGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984291‑2984591

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