Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,4410AC38.8% 115.6 / 22.7 67G384G (GGT→GGGyqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (38/3);  new base C (0/26);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 1.35e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGA  >  NC_000913/2984292‑2984588
                                                                                                                                                     |                                                                                                                                                    
gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCacac                                                                                                                                                     >  1:355290/1‑151 (MQ=255)
gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCacac                                                                                                                                                     <  1:1332184/151‑1 (MQ=255)
 aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAGCAGATTAATGACGGCAGATATTTATCGGTAAGGAGTAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAcc                                                                                                                                                    <  2:352530/151‑1 (MQ=255)
 aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAcc                                                                                                                                                    >  1:1141464/1‑151 (MQ=255)
    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCt                                                                                                                                                  <  1:163119/150‑1 (MQ=255)
      gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTgg                                                                                                                                                <  1:878704/150‑1 (MQ=255)
      gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGc                                                                                                                                              >  1:822435/1‑151 (MQ=255)
          tGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGccccc                                                                                                                                          >  2:473247/1‑151 (MQ=255)
                tGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGAt                                                                                                                                    >  2:649207/1‑151 (MQ=255)
                       aTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGc                                                                                                                             >  1:1093590/1‑151 (MQ=255)
                           gggAGCCAAAAAAATTTATGTCGCTTTGCCACCCATATTTATTAAAAAAATAAATGACGAAAGATTTTTACCGGAAGGGAGAAAATGCAGAGCTTTTGGTCACTTTCATTGGTAAAGCCCCCCCCCGGGCCCCCGCCGATACAAAGCGtt                                                                                                                          <  2:1232250/150‑1 (MQ=255)
                                   aaaaaaTTTATGCCGTCTTCCCAGCCATCTTTATTAAAAAGATTAATGACGGCAGATATTTACGGGTAAGGAGAAAATGCAGAGTTTTTCTTCACGTTAAATGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTcc                                                                                                                  <  1:1244985/148‑1 (MQ=255)
                                   aaaaaaTGTATGTCGACTTGCCAGCCATCATTATTAAACAGATTAAGGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCCCTTTAAATGTTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                                 <  1:572257/149‑1 (MQ=255)
                                   aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                                 >  1:653869/1‑151 (MQ=255)
                                    gAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCtt                                                                                                                >  1:1004700/1‑151 (MQ=255)
                                     aaaaTGTATGTCGTCTTGCCACCCATAATTATTAAAAAGATTAATGCCGGAAGATTTTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCCCGTTAAAGGATAAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                               <  1:50530/151‑1 (MQ=255)
                                       attgtttgtCGTCTTGCCAGCCTTTATTATAAAAGAGATAAAGGACGGAAGATTTTTAGCGGAAAGGAGAAAAGGAAGAGCTTTTCGTCACGTTAATTGATAAAGCCCCCCCCGTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                                                             <  1:190210/149‑1 (MQ=255)
                                         tgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg                                                                                                           >  2:219535/1‑151 (MQ=255)
                                             tgtctttttGCCAGCCCTAATTATACAAGATTTTAAGGGCGCTATTTTTTTTCCGGTAAGGAGAAAACGCAGATCTATTCTCCGCTTTAAAGGGTCATTACCACCCCCCGGCCCCCGCCGATTCAAAGCGTTGCCAGTCCTTTGCGGGc                                                                                                         <  1:298003/142‑1 (MQ=255)
                                             tgccgccTTGCCAGCAATATTTATACAAAAAATAAATGCCGGAAAATTTTTAGGGGAAAGGAGAAAATGCAGCCTTTTCCGCCACTTTAATTGTCAATCCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                                                             <  2:252264/139‑1 (MQ=255)
                                              gtcgtcTTGCCAGCAATAATTATTCAAAAAATTAATGACGGCAGATTTTTACGGGAAAGGAGAAAAGGCAGAGTTATTCGTCACGTTAAAGGGCAATTCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCa                                                                                                        <  1:309193/149‑1 (MQ=255)
                                                cGACTTGCCAGCAATAATTATTAAACAGATAAATGACGGAAGTTTTTTATGGGAAAGGAGAAAATGCAGAGCTTTTCGTCCCTTTCAATGGTAAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc                                                                                                    <  2:33465/151‑1 (MQ=255)
                                                 gtcTTGCCAGCCATTTTTATTCAACAGATAAAGGACGGAAGATTTTTATGGGAAGGGAGAAAATGCAGATTTTTTCGCCACGTTAAAGGGTCAACCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc                                                                                                    <  1:211531/148‑1 (MQ=255)
                                                  aCTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                                   >  1:193879/1‑150 (MQ=255)
                                                    ttGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt                                                                                                >  2:727434/1‑151 (MQ=255)
                                                    ttGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt                                                                                                >  2:1068449/1‑151 (MQ=255)
                                                     tGCCAGCCATCTTTATTCAACAGATTAATGACGGCAGATATTTATGGGTAAGGAGAAAATGCAGAGTTATTCGTCACTTTCAATGGTAAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGt                                                                                                  <  1:472573/148‑1 (MQ=255)
                                                        cagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                            >  2:305618/1‑151 (MQ=255)
                                                            catcatTATCCAAAAGTTAAAGGTCGCCAATTTTTTACCGGTAGGAAGGAAATGAAAGCCTTTTCGTCAGGTTCAATGGTCAATGCCCCCCCCTGCCCTCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt                                                                                         <  1:860018/150‑1 (MQ=255)
                                                              taattatcCAACAGATTAAGGACGGAAGATTGTTAGCGGTAGGGAGAAAAGGCAGGCCTTTTCGCCACGTTAAATGTTCAGCACCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt                                                                                        <  2:1272635/147‑1 (MQ=255)
                                                               catTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                     >  1:1092893/1‑151 (MQ=255)
                                                                       aaCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                                >  2:784224/1‑148 (MQ=255)
                                                                       aaCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                                >  2:97756/1‑148 (MQ=255)
                                                                       aaCAGATTAATGACGGAAGATTTTTAGCGGAAAGGAGAAAATGCAGAGTTTTTCGTCCCTTTAAAGGGTCAACCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                              <  1:42076/150‑1 (MQ=255)
                                                                               aaTGACGGCAGATTTTTATGGAGAAGGAGGAAACGCAAAGCTATTCGCCACGTTCAGTGGTAAGCCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa                                                                                   <  2:859243/137‑1 (MQ=255)
                                                                                 tGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                    >  1:967334/1‑150 (MQ=255)
                                                                                    cGGAAGATTTTTACCGGAAAGGAGAAAATGCAGACCTTTTCGTCGGTTTCAAGGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc                                                                <  1:1258065/151‑1 (MQ=255)
                                                                                           tatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                         >  2:709332/1‑151 (MQ=255)
                                                                                           tatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                         >  1:195555/1‑151 (MQ=255)
                                                                                            atTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCg                                                        >  1:106622/1‑151 (MQ=255)
                                                                                             tttATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                       >  1:687997/1‑151 (MQ=255)
                                                                                               tATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                      >  2:263834/1‑150 (MQ=255)
                                                                                                atgggTAAGGAGAAAATGCAGACCTTTTCGTCACGTTCAATGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg                                                    <  2:171616/148‑1 (MQ=255)
                                                                                                 tCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc                                                   >  1:1209945/1‑151 (MQ=255)
                                                                                                    gTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGcc                                                >  2:526129/1‑151 (MQ=255)
                                                                                                      aaGGAGAAAATGCAGACCTTTTCGTCGCGTTAAAGGGTCAACGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                               <  1:571027/150‑1 (MQ=255)
                                                                                                        ggaggaAATGCAGAGCTATTCGCCCCTTTAAATGGTCAACGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTAt                                             <  2:1175548/150‑1 (MQ=255)
                                                                                                        ggagaaAATGCAGACCTTTTCGTCCCTTTAAATGTTAAAGCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                                  <  1:856171/145‑1 (MQ=255)
                                                                                                         gaggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAGCCAGGATTCGGCAACCGGAAGCGCCTATCg                                           >  1:732385/1‑151 (MQ=255)
                                                                                                          aggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                          >  1:1144525/1‑151 (MQ=255)
                                                                                                           ggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg                                         >  2:447757/1‑151 (MQ=255)
                                                                                                            gaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg                                         >  2:1309796/1‑150 (MQ=255)
                                                                                                             aaaTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                       >  2:309866/1‑151 (MQ=255)
                                                                                                              aaTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                       >  2:1095587/1‑150 (MQ=255)
                                                                                                              aaTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                       >  2:355898/1‑150 (MQ=255)
                                                                                                                 gCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAcc                                    >  1:1238333/1‑150 (MQ=255)
                                                                                                                   agagCTATTCGTCGCGTTAAATGTTCAACCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGa                                  <  2:956990/150‑1 (MQ=255)
                                                                                                                   agagCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                   >  2:792818/1‑149 (MQ=255)
                                                                                                                           tCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                         >  2:112079/1‑151 (MQ=255)
                                                                                                                            cGTCCCTTTCAGTGGTCAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                          <  1:978059/149‑1 (MQ=255)
                                                                                                                                  tttaagcgatgaaTGCCCCCCCCCGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGAt                           <  2:550898/131‑1 (MQ=255)
                                                                                                                                     cAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcatt               >  2:276398/1‑151 (MQ=255)
                                                                                                                                     cAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcatt               >  2:231495/1‑151 (MQ=255)
                                                                                                                                      agtggaaaaTGCCCCCCCCCGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcatt               <  1:219541/143‑1 (MQ=255)
                                                                                                                                      aaTGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg              >  2:381707/1‑151 (MQ=255)
                                                                                                                                      aaTGGTCAAGCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacca                 <  2:481538/148‑1 (MQ=255)
                                                                                                                                         ggTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgaca           <  1:864818/151‑1 (MQ=255)
                                                                                                                                             aaTGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg        >  2:539978/1‑150 (MQ=255)
                                                                                                                                             aaTGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg        >  2:881649/1‑150 (MQ=255)
                                                                                                                                                  caccCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGa  <  1:1095108/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGA  >  NC_000913/2984292‑2984588

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