Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 436,419 T→C G262G (GGT→GGC)  thiL → thiamine monophosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913436,4190TC93.8% 208.3 / 2.4 68G262G (GGT→GGC) thiLthiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base C (37/23);  total (40/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

CTGACGGTTTGATTTCCGATCTCGGGCATATCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACC  >  NC_000913/436274‑436565
                                                                                                                                                 |                                                                                                                                                   
ctGACGGTTTGATTTCCGATCTCGGGCATATCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGaa                                                                                                                                                  <  1:1000694/149‑1 (MQ=255)
ctGACGGTTTGATTTCCGATCTCGGGCATATCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGa                                                                                                                                                   >  2:76349/1‑148 (MQ=255)
  gACGGTTTGATTTCCGATCTCGGGCATATCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAg                                                                                                                                                 <  1:1155969/148‑1 (MQ=255)
  gACGGTTTGATTTCCGATCTCGGGCATATCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGAt                                                                                                                                               <  1:8857/150‑1 (MQ=255)
    cGGTTTGATTTCCGATCTCGGGCATATCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTAc                                                                                                                                            >  1:909622/1‑151 (MQ=255)
     ggTTTGATTTCCGATCTCGGGCATATCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGAtca                                                                                                                                             >  1:1064953/1‑147 (MQ=255)
     ggTTTGATTTCCGATCTCGGGCATATCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACg                                                                                                                                           >  1:1002133/1‑151 (MQ=255)
      gTTTGATTTCCGATCTCGGGCATATCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGa                                                                                                                                          <  2:420484/151‑1 (MQ=255)
       tttGATTTCCGATCTCGGGCATATCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACg                                                                                                                                           >  1:445308/1‑149 (MQ=255)
        ttGATTTCCGATCTCGGGCATATCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGt                                                                                                                                        >  1:611928/1‑151 (MQ=255)
             ttCCGATCTCGGGCATATCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTgtgt                                                                                                                                   >  1:323869/1‑151 (MQ=255)
               ccGATCTCGGGCATATCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGttt                                                                                                                                 >  2:574792/1‑151 (MQ=255)
                      cGGGCATATCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACt                                                                                                                             <  1:1101727/148‑1 (MQ=255)
                          cATATCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCgg                                                                                                                      >  2:9626/1‑151 (MQ=255)
                           atatCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGa                                                                                                                     <  1:490069/151‑1 (MQ=255)
                            tatCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGa                                                                                                                     >  1:412552/1‑150 (MQ=255)
                             atCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAAc                                                                                                                   >  2:1018296/1‑151 (MQ=255)
                               cGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACt                                                                                                                  <  1:532075/150‑1 (MQ=255)
                                     aGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑Gt                                                                                                          <  2:700004/151‑1 (MQ=255)
                                        cAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGgc                                                                                                       >  2:548671/1‑151 (MQ=255)
                                         aGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGgcg                                                                                                      <  2:52996/151‑1 (MQ=255)
                                           cGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGgcgcg                                                                                                    >  1:681699/1‑151 (MQ=255)
                                           cGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGgcgcg                                                                                                    <  1:179961/151‑1 (MQ=255)
                                            gACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGgcgcgc                                                                                                   >  1:1127523/1‑151 (MQ=255)
                                             aCTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCt                                                                                                  >  1:815020/1‑151 (MQ=255)
                                                 cggcgCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGAt                                                                                              <  2:1247096/151‑1 (MQ=255)
                                                     gcACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTgg                                                                                          <  2:483438/151‑1 (MQ=255)
                                                               aCCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGaca                                                                                <  1:257237/151‑1 (MQ=255)
                                                               aCCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGac                                                                                 >  1:302850/1‑150 (MQ=255)
                                                                 cTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACAcc                                                                              >  2:1159114/1‑151 (MQ=255)
                                                                  tGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACAcc                                                                              >  1:1083647/1‑150 (MQ=255)
                                                                  tGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCt                                                                             >  1:95916/1‑151 (MQ=255)
                                                                  aggCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACAcc                                                                              >  1:463316/2‑150 (MQ=255)
                                                                   ggCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACAcc                                                                              <  2:1069709/149‑1 (MQ=255)
                                                                     cATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTggg                                                                          <  1:1175438/151‑1 (MQ=255)
                                                                      aTTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTggg                                                                          >  2:342905/1‑150 (MQ=255)
                                                                        tgctgcCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGt                                                                       >  2:699143/1‑151 (MQ=255)
                                                                         gctgcCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTa                                                                      <  2:399911/151‑1 (MQ=255)
                                                                           tgcCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTGcc                                                                    >  2:240226/1‑151 (MQ=255)
                                                                           tgcCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTAcc                                                                    >  2:112877/1‑151 (MQ=255)
                                                                            gcCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTAc                                                                     <  2:801328/149‑1 (MQ=255)
                                                                            gcCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCg                                                                   <  2:396555/151‑1 (MQ=255)
                                                                                ttttCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGttt                                                                >  1:418854/1‑150 (MQ=255)
                                                                                    cTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTg                                                           >  2:894822/1‑151 (MQ=255)
                                                                                         gcgcTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATc                                                       <  2:763372/150‑1 (MQ=255)
                                                                                          cgcTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCgg                                                     >  2:435289/1‑151 (MQ=255)
                                                                                              ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc                                                   <  2:1216230/149‑1 (MQ=255)
                                                                                              ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc                                                   <  2:1137177/149‑1 (MQ=255)
                                                                                                ctcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCaaa                                                >  1:977453/1‑150 (MQ=255)
                                                                                                ctcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAt                                               >  2:1058951/1‑151 (MQ=255)
                                                                                                  cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAt                                               <  2:1238369/149‑1 (MQ=255)
                                                                                                     cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccg                                          >  2:45326/1‑151 (MQ=255)
                                                                                                     cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccg                                          >  2:529816/1‑151 (MQ=255)
                                                                                                       tGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccgcc                                        >  1:722452/1‑151 (MQ=255)
                                                                                                                   aCAGGCGCTGCGCTGGGCGCTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAggg                            >  2:825322/1‑151 (MQ=255)
                                                                                                                      ggCGCTGCGCTGGGCGCTCTCGGGCGGGGAAGATTGCGGGTTGTGTTTCACTGTGCCGGAACTGACCC‑GGGGCGCGCGGGTTGGGGCTCGCGGCCCCCTGGGCGTACCTTTTCCTTGTAGGGGGAAAATGACCGCCGATATCGAAGGGCtt                         <  2:52406/151‑1 (MQ=255)
                                                                                                                                tgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTAtt                >  1:884843/1‑150 (MQ=255)
                                                                                                                                 gggCGCTCTGGGGGGGGGAAGATGAAGAGTTGTTTTTCACTGTCCCCGAA‑TGAACCGGGGGGGCGCTGGATGGGGCTCTGGAACACGTGGGCGCCCCGTTTCCCTTTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTa                  <  2:55040/147‑1 (MQ=255)
                                                                                                                                 gggCGCTCTCTGGCGGCAAAGATTAGGAGTTGTTTTTCCCTGTGCCGAAACTGACCC‑GGGGCGCGCTGGATGTGGCTCCCGGACCCCGGGGAGTCCCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGtt                     <  2:740571/144‑1 (MQ=255)
                                                                                                                                 gggCGCTCTCGGGGGGGGAATATTACGTTTTGTTTTTCACTGTGCCGGACCTGACCC‑GGGGCGCGCTGGTTGTGGCTCTCGGACCCCGGGGCGTCCTTTTTCCCTTTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGtttt                   <  1:1092183/146‑1 (MQ=255)
                                                                                                                                 gggCGCTCTCGGGGGGCAAAGTTTACGGTTTGTGTTTCTGTGTCCCGGAACTGACCC‑GGGGCGCGCTGGATGTGCCTCTGGGACCCCGGGGCGTCCTTTTTCCTTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTAtt                <  2:1246633/149‑1 (MQ=255)
                                                                                                                                 gggCCCTTTCGGGCGGGAAATATTACGAGTTGTTTTCCGTGGTGCCGGAACTGACCC‑GGGGCGCGTTGGTTGGGGTTCGGGAACCCTGGGGCCTCCCTTTTCCCTTTAGCGGAAAATTGACCGCCGATATCGAAGGGCTTTGTTTTATTc               <  2:61375/150‑1 (MQ=255)
                                                                                                                                   gcgcTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg            >  1:642858/1‑151 (MQ=255)
                                                                                                                                     gcTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGAc          >  2:1260756/1‑151 (MQ=255)
                                                                                                                                     gcTCTCTGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGAc          >  2:122897/1‑151 (MQ=255)
                                                                                                                                      ctctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg         >  2:48439/1‑151 (MQ=255)
                                                                                                                                      ctctctGGCGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg         >  1:1055656/1‑151 (MQ=255)
                                                                                                                                            ggcggcGAAAATTACGTTTTGTGTTCCGTTGTCCGGGAACTCCGCC‑GTGGCGCGCGGGATGGGGCTTTCGGACCCCTGGGCTTTCGTTTTCCCGTTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg      <  2:802227/148‑1 (MQ=255)
                                                                                                                                              cGGCGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAcc  >  1:953385/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTGACGGTTTGATTTCCGATCTCGGGCATATCGTGAAAGCCAGCGACTGCGGCGCACGTATTGACCTGGCATTGCTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACC  >  NC_000913/436274‑436565

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