Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,312,2130GA38.8% 109.4 / 47.5 81L12F (CTT→TTT) yciAacyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base G (9/40);  new base A (31/0);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 3.82e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTAT  >  NC_000913/1312072‑1312361
                                                                                                                                              |                                                                                                                                                         
cccGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATCGGTAAAAAACGTAAAAAAAAATCGccc                                                                                                                                                    >  2:1315031/1‑149 (MQ=255)
        gTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTCAAAAACGCAAAACAAAAGCGCCCCGCggg                                                                                                                                              >  1:1047077/1‑147 (MQ=255)
              aCGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGGGAAAATACGTAAAAAAAAAAAGCCCTTGGGGGCGGTAtgtg                                                                                                                                    >  2:998312/1‑151 (MQ=255)
              aCGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAATACGGAAAAAAAAAACGCCCTGAGGGGCggttt                                                                                                                                       >  2:248433/1‑144 (MQ=255)
               cgcgACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAAAGTAAAAAAAAAACGGCCCGAgggggggtat                                                                                                                                       >  2:1029873/1‑140 (MQ=255)
                 cgACCGTGGGCAATTTCTTTTGCCAGAATATCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATTGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGtgtt                                                                                                                                  <  1:1105281/150‑1 (MQ=255)
                  gACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGtgttgt                                                                                                                               <  1:959285/151‑1 (MQ=255)
                    ccGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTa                                                                                                                              <  2:61957/150‑1 (MQ=255)
                      gTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATAGCCATTTGTGACATTAACCAACCACCAAATATGTCACCATTGGCATTGGTGTCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGa                                                                                                                           <  1:1169702/150‑1 (MQ=255)
                         ggCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATg                                                                                                                      <  1:117148/151‑1 (MQ=255)
                          gCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATgg                                                                                                                     <  1:230017/151‑1 (MQ=255)
                          gCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATgg                                                                                                                     <  2:853566/151‑1 (MQ=255)
                             aTTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTaa                                                                                                                 <  1:694484/151‑1 (MQ=255)
                               ttCTTTTGTTAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTAGCATTGGTATCGGCGGGCATGGCTAAAGTACGTACAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTaaa                                                                                                                <  2:854031/150‑1 (MQ=255)
                               ttCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAATAAGGAAAAAAAAAACGCCCCGAGGGACGTTATGTGTTTTAGACATGGGaaa                                                                                                                >  2:1006835/1‑150 (MQ=255)
                                   tttGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGACATGGCTAAAAGAAGTAAAAAAAAAAAGCCCCGAGGGGGGGTTTGTGTTGTTGAAAAGGGaaaa                                                                                                               >  2:1302001/1‑147 (MQ=255)
                                     tGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGGTAAAATAAGAAAAAAAAAAATGCCCTGAGGGGCGGTATTTGTTGTTGACATGGGaaaa                                                                                                               >  2:903968/1‑145 (MQ=255)
                                      gCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACtt                                                                                                       <  2:982216/151‑1 (MQ=255)
                                         aGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCAAGGCTAAAGAACGGAAAACAAAACCACCCCGAAGGGCGGTATGTATTGTTGGCAAGGTaaaaac                                                                                                             >  2:647343/1‑141 (MQ=255)
                                               gcgcCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGaaa                                                                                              <  2:493232/151‑1 (MQ=255)
                                                cgccgcCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGGTAAAATACGTAAAACAAAAACGCCCTGAGGGAGGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAAGaaa                                                                                              >  2:1071337/1‑150 (MQ=255)
                                                cgccgcCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCAATGCTAAAATACGTAAAAAAAAATCGCCCTGAGGGGCGGTTTTTGTTGTGGAGAAGGTAAAAACGAATTAAAAGGaaa                                                                                              >  2:493733/1‑150 (MQ=255)
                                                 gccgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGCTAAAATACGTAAAAAAAAAACGCCCCGAGGGGCGTTATGTGTTGTAGACATGGTAAAAACGACTTAAAAAGaaaa                                                                                             >  2:487070/1‑150 (MQ=255)
                                                   cgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATca                                                                                          <  2:890694/151‑1 (MQ=255)
                                                    gccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGGTGGCGAAAAAAAGTAAAACAAAAACGCCCCGAGGGGGGGTTTGTGTTGTAGAGAAGGGAAAAACCACCTAAAAAAAAAAAca                                                                                          >  2:513426/1‑150 (MQ=255)
                                                    gccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTAAGTAAAAAAAAAACGCCCCGAGGGGGGTTTTGTGTTGTAGAAAAGGGAAAAACCACCTAAAAAGaaaa                                                                                             >  2:1067701/1‑147 (MQ=255)
                                                     ccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGAACGTAAAACAAAATCGCCCTGAGGGGCGTTATTTGTTGTTGACATGGTAAAAACGACTTAAAAAGAAAAacac                                                                                         >  2:1065960/1‑150 (MQ=255)
                                                      cAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAg                                                                                       <  2:898625/151‑1 (MQ=255)
                                                       aaTATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAgg                                                                                      <  2:313109/151‑1 (MQ=255)
                                                       aaTATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAACGTAAAAAAAAAACGCCCCGGGGGGCGGTATTTGTTGTTGAAAAGGGAAAAACGACTTAAAAAGAAAAAca                                                                                          >  2:132843/1‑147 (MQ=255)
                                                        atatCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGAAAGCTCAAAAAATGAAAACAAAACTGCCCCGAGGGGCGGTTTGTGTTTTAGGAAAGGGAAAAAAACCTTTAAAAGGAAAAta                                                                                           >  1:749815/1‑145 (MQ=255)
                                                          atCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGgcgc                                                                                   <  2:169052/151‑1 (MQ=255)
                                                            ccATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCAATGCCCAAACAAAGAAAAACAAAAACTCATGAGGGACGGTATGTTTTGTAGAAATGGGAAAAACAACTTTAaaaaaaaaataccggagggc                                                                                 >  1:1168506/1‑139 (MQ=255)
                                                            ccATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCGAAAATACGGAAAAAAAAAACGCCCTGAGGGAGGCTTTGGGTTGTTGAAAAGGGAAAAAAGACTTTAAAAGGaaaaaaca                                                                                         >  2:254216/1‑140 (MQ=255)
                                                             cATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGAt                                                                                <  1:810759/151‑1 (MQ=255)
                                                             cATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAAAGTAAAAAAAAAAAGCCCTGAGGGGCGGTATGTGTTGTAGAAAAGGGAAAAACGACTTAAAAAGaaaaaaaaagggg                                                                                    >  1:24257/1‑138 (MQ=255)
                                                              aTTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAAAAAGTAAAAAAAAATCGCCCCGAGGGGCGCTATTTGTTGTAGAAATGGTAAAAACGACTTAAAAAGAAAAACAAAGGCGCGa                                                                                 >  1:818210/1‑149 (MQ=255)
                                                                ttGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGGTTAAATACGCAAAAAAAAAACGCCCCGAGGGGGGGTTTGTGTTTTTGACATAGGAAAAACGAAATTAAAAAAAAAATACAGGcgcggg                                                                                >  2:15519/1‑146 (MQ=255)
                                                                  gtgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCGGGCCCAAAAACGTAAAACAACAACGCCCTGAGGGGCGATATGTGTTGTGGACAAGGTGGAAACCACCTAAAAAGAAAAAAACAGgcggg                                                                                  >  2:317351/1‑142 (MQ=255)
                                                                    gACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGctact                                                                          <  2:1126062/150‑1 (MQ=255)
                                                                     aCATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTAt                                                                        <  2:202911/151‑1 (MQ=255)
                                                                     aCATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTAt                                                                        <  2:1330146/151‑1 (MQ=255)
                                                                       aTTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGTGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAGCCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGa                                                                      <  1:580070/151‑1 (MQ=255)
                                                                          aaccaaccACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGAttt                                                                   <  2:930319/151‑1 (MQ=255)
                                                                          aaccaaccACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCAAAGCCAAAAAAAGTAAAAAAAAAAACCCCTGAGGGGCGGTTTGTGTTGTTGGAAATGGGAAAACGAATTAAAAAAAAAAAAAAAGGCGCGGAGCTCATATTAttt                                                                   >  1:938908/1‑151 (MQ=255)
                                                                           accaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCTGGCTAAAGAACGTAAAAAAAAAAAGCCCCGAGGGGCGGTATGTGTTGTTGAAATGGTAAAAACGACTTTAAAAGAAAAACACAAGCGCGGTGCTACTATGAtttt                                                                  >  2:981980/1‑151 (MQ=255)
                                                                           accaaacaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGAttt                                                                   <  1:550200/150‑1 (MQ=255)
                                                                            ccaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTc                                                                 <  2:604601/151‑1 (MQ=255)
                                                                                 caccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTgag                                                            <  2:1204216/151‑1 (MQ=255)
                                                                                  accaAAGATGTCACCATTGGCATTGGTATCGGCGGGGCTGGGTAAAAAAAGTAAAAAAAAAACGTCCCGTGGGGGGGTATGTGTTGTTGACAAGGGAAAAAAGAACTAAAAGGGAAAAAACAGGAGCGCTGGGAATATGATTTTTCTgtgg                                                           >  1:138727/1‑151 (MQ=255)
                                                                                   ccaAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGGCAAAGTAAGTAAAACAAAAAAGCCCCGCGGGGGGGTATGTGGTGTGGAAATGGGAAAAAAGACTTAAAAAGAAAAACACAAGCGCGAGGCTACTATGATTTTTCTtag                                                            >  2:699114/1‑149 (MQ=255)
                                                                                    caAAGATGTCACCATTGGCATTGGTATCGGCGGGGACGGGTAAAAAACGTAAAACAAAAACGCCCCGCGGGGCGGTATGTGTTGTAGAAATGGAAAAAAACGATTAAAAGGAAAATCACAAGGGCGACGCTCCTATTATTTTTCTtag                                                            >  1:1045859/1‑148 (MQ=255)
                                                                                      aaGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGaa                                                       <  2:571002/151‑1 (MQ=255)
                                                                                       aGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAAc                                                      <  2:282382/151‑1 (MQ=255)
                                                                                        gATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACa                                                     <  2:1047850/151‑1 (MQ=255)
                                                                                        gATGTCACCATTGGCATTGGTATCGGCGGGCAAGGCTAAAATAAGGAAAACAAAAACGTCCCGAGGGGCGCTATGGGGTGTAGACATGGGAAAAACAACTTTAAAGGAAAAAAAAAGGGGCGGTGCTACTATGATTTTTCTgatggg                                                         >  1:1195065/1‑143 (MQ=255)
                                                                                         aTGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACag                                                    <  2:1209628/151‑1 (MQ=255)
                                                                                           gTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACagag                                                  <  1:908328/151‑1 (MQ=255)
                                                                                           gTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACaga                                                   <  1:777390/150‑1 (MQ=255)
                                                                                              aCCATTGGCATTGGTATCGGCGGGGATGGGGAAAATACCTAAAACAAAATCGCCCTCAGGGGGGGGATGTGTTGTAGACCAGGGAAAAACGACTTTAAAAGAAAAACACAGGGGCGGGGCTACTAAGATTTTTCTTAGGAGaaaa                                                     >  1:1151684/1‑143 (MQ=255)
                                                                                              aCCATTGGCATTGGTATCGGCGGGCATGGCCAAAAAAAGTAAAAAAAAAACGCCCCGAGGGACGGTTTGTGTTTTAGAAAAGGGAAAAACGACCTAAAAAGAAAATCACAAGCGCGATGCTACTATGATTTTTCTGAGGAGAACAgagg                                                 >  2:1052511/1‑149 (MQ=255)
                                                                                                    ggcattggTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAg                                         <  2:285029/151‑1 (MQ=255)
                                                                                                       attggTATCGGCGGGCATGGCTCAAAAACGGAAAAAAAAAACGCCCTGAGGGGCGTTATGTGGTGTGGACAAGGTAAAAACGAATTAAAAAGAAGAACACAAGGGCGAAGATACTTTGATTTTTCTGGGGGGAAAAGAGGAGGTACa                                          >  1:330612/1‑147 (MQ=255)
                                                                                                          ggTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCa                                    <  2:384100/150‑1 (MQ=255)
                                                                                                          ggTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAg                                   <  1:1060512/151‑1 (MQ=255)
                                                                                                               cggcggGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGcccc                               >  2:732978/1‑150 (MQ=255)
                                                                                                                 gcggGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGaa                            <  1:1171001/151‑1 (MQ=255)
                                                                                                                 gcggGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGaa                            <  2:829144/151‑1 (MQ=255)
                                                                                                                  cggGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGa                             >  1:726031/1‑149 (MQ=255)
                                                                                                                    ggCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGa                         <  2:1163721/151‑1 (MQ=255)
                                                                                                                        aTGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTg                      >  2:315491/1‑151 (MQ=255)
                                                                                                                          ggCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGc                    >  2:1338571/1‑151 (MQ=255)
                                                                                                                            cTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGc                    >  1:432456/1‑149 (MQ=255)
                                                                                                                             tAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGt                 >  2:135339/1‑151 (MQ=255)
                                                                                                                              aaaGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTc                >  1:40508/1‑151 (MQ=255)
                                                                                                                                 gTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTgg             <  2:557506/151‑1 (MQ=255)
                                                                                                                                       aaaaCAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAgg        <  1:411259/150‑1 (MQ=255)
                                                                                                                                       aaaaCAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGa       <  2:24261/151‑1 (MQ=255)
                                                                                                                                        aaaCAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGAt      <  1:944728/151‑1 (MQ=255)
                                                                                                                                          aCAAGATCGCCCTGAGGGACGTTATGTGTTGTCGACATGGTAAAACCGACTTAACAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGAttt    >  1:1041409/1‑151 (MQ=255)
                                                                                                                                            aaGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTAt  >  2:783679/1‑151 (MQ=255)
                                                                                                                                              |                                                                                                                                                         
CCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTAT  >  NC_000913/1312072‑1312361

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: