Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,611,827 T→C V124A (GTG→GCG)  hyfJ → putative processing element hydrogenase 4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,611,8270TC98.6% 256.6 / ‑4.2 72V124A (GTG→GCG) hyfJputative processing element hydrogenase 4
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (27/44);  total (28/44)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

GCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGTGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGTTTGCCCCGCCACAAGGAATAACAATTGAAGCGGTAAACGGAATGCTCG  >  NC_000913/2611681‑2611972
                                                                                                                                                  |                                                                                                                                                 
gCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTgcgc                                                                                                                                                 >  2:792123/1‑149 (MQ=255)
gCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCa                                                                                                                                                <  2:714543/150‑1 (MQ=255)
 cTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAgg                                                                                                                                              >  1:1415541/1‑151 (MQ=255)
  tGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAgg                                                                                                                                              >  1:710565/1‑150 (MQ=255)
   gggggTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGaa                                                                                                                                            <  1:165217/151‑1 (MQ=255)
    ggggTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGaa                                                                                                                                            >  1:1061837/1‑150 (MQ=255)
           tgaCTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCat                                                                                                                                    >  1:952484/1‑151 (MQ=255)
              cTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCatct                                                                                                                                  >  2:425200/1‑150 (MQ=255)
                tttGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCatcta                                                                                                                                 <  2:1371816/149‑1 (MQ=255)
                 ttGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCatctatc                                                                                                                               <  1:937473/150‑1 (MQ=255)
                      gagaTTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATg                                                                                                                          <  1:666689/150‑1 (MQ=255)
                      gagaTTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATg                                                                                                                          <  2:407994/150‑1 (MQ=255)
                       agaTTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATg                                                                                                                          >  2:18434/1‑149 (MQ=255)
                       agaTTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGc                                                                                                                        >  1:1392216/1‑151 (MQ=255)
                       agaTTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGc                                                                                                                        <  1:1255614/151‑1 (MQ=255)
                         aTTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCa                                                                                                                      <  1:838605/151‑1 (MQ=255)
                            gTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGaa                                                                                                                   <  2:1186929/151‑1 (MQ=255)
                             ttATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGaa                                                                                                                   <  2:1032393/150‑1 (MQ=255)
                              tATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGaaa                                                                                                                  >  2:287437/1‑150 (MQ=255)
                              tATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGaaa                                                                                                                  <  1:1043933/150‑1 (MQ=255)
                              tATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGaaa                                                                                                                  <  1:936476/150‑1 (MQ=255)
                               aTTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGaaaaa                                                                                                                <  1:1102036/151‑1 (MQ=255)
                               aTTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGaaaa                                                                                                                 <  1:32685/150‑1 (MQ=255)
                               aTTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGaaaa                                                                                                                 >  1:659677/1‑150 (MQ=255)
                                 tGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAAtt                                                                                                              >  1:463072/1‑151 (MQ=255)
                                 tGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAAt                                                                                                               <  2:1415180/150‑1 (MQ=255)
                                  gACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAAtt                                                                                                              <  1:325379/150‑1 (MQ=255)
                                  gACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAAtt                                                                                                              <  2:667122/150‑1 (MQ=255)
                                   aCGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGc                                                                                                            >  1:424111/1‑151 (MQ=255)
                                   aCGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGc                                                                                                            <  1:813734/151‑1 (MQ=255)
                                     gccagccaCACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGt                                                                                                          <  2:41598/151‑1 (MQ=255)
                                      ccagccaCACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTg                                                                                                         <  1:444252/151‑1 (MQ=255)
                                          ccaCACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGaa                                                                                                     <  2:586548/151‑1 (MQ=255)
                                           cacaCCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAg                                                                                                    <  2:639633/151‑1 (MQ=255)
                                            acacCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAgg                                                                                                   <  1:1374604/151‑1 (MQ=255)
                                              acCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGAt                                                                                                 >  1:261731/1‑151 (MQ=255)
                                               ccGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAgg                                                                                                   <  1:810203/148‑1 (MQ=255)
                                                 gCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGAttt                                                                                               >  2:1117364/1‑150 (MQ=255)
                                                 gCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTtc                                                                                              <  2:501599/151‑1 (MQ=255)
                                                    cTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTtctc                                                                                            >  2:1123098/1‑150 (MQ=255)
                                                            ccGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATtaata                                                                                     >  1:109292/1‑149 (MQ=255)
                                                             cGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAgg                                                                                  <  1:1218915/151‑1 (MQ=255)
                                                               ggCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGAc                                                                                <  2:291970/151‑1 (MQ=255)
                                                                     cgcgCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGtt                                                                            <  2:710486/149‑1 (MQ=255)
                                                                       cgcCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATg                                                                        <  2:1339963/151‑1 (MQ=255)
                                                                         ccACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGc                                                                      <  1:792201/151‑1 (MQ=255)
                                                                          cACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCt                                                                     <  1:282602/151‑1 (MQ=255)
                                                                              ggCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATg                                                                  <  2:246433/150‑1 (MQ=255)
                                                                                cGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTc                                                                >  2:1402032/1‑150 (MQ=255)
                                                                                   aTGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGAc                                                            <  2:1235141/151‑1 (MQ=255)
                                                                                   aTGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGAc                                                            <  2:728540/151‑1 (MQ=255)
                                                                                          gACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCg                                                      <  2:886940/150‑1 (MQ=255)
                                                                                             gTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATg                                                   <  2:953001/150‑1 (MQ=255)
                                                                                                 gtgtgGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGttt                                                <  1:1219108/149‑1 (MQ=255)
                                                                                                 gtgtgGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGtt                                                 <  2:215381/148‑1 (MQ=255)
                                                                                                 gtgtgGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGTTTGc                                              <  1:798819/151‑1 (MQ=255)
                                                                                                    tgGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGTGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGTTTGcccc                                           >  1:1250266/1‑151 (MQ=255)
                                                                                                    tgGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGTTTGcccc                                           >  2:1015762/1‑151 (MQ=255)
                                                                                                     ggCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGTTTGcccc                                           <  2:339966/150‑1 (MQ=255)
                                                                                                      gCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGTTTGCCCCGc                                         <  1:291111/151‑1 (MQ=255)
                                                                                                         ggCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGTTTGCCCCGCca                                       <  2:1307639/150‑1 (MQ=255)
                                                                                                                       aattcaTTCATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGTTTGCCCCGCCACAAGGAATAACAAtt                        >  1:1184785/1‑151 (MQ=255)
                                                                                                                        attcattcATCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGTTTGCCCCGCCACAAGGAATAACAATTg                       >  1:1368265/1‑151 (MQ=255)
                                                                                                                         ttcattcaTCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGTTTGCCCCGCCACAAGGAATAACAATTg                       >  1:589195/1‑150 (MQ=255)
                                                                                                                            attcatCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGTTTGCCCCGCCACAAGGAATAACAATTGAAGc                   >  1:434267/1‑151 (MQ=255)
                                                                                                                             ttcatCTGTTGGAAGAAATTGCGCAGGAACCGGCCTTCTATCTGATGGAAAAAAAAATGGCGTGAGAAGTATTTCTCATAAAAAGGGATTGTTAATGGTTATGGAAGACGGGGCGTTTTTGCCCCCGCCACAAGGAATAACAATTGAAGc                   <  2:815721/150‑1 (MQ=255)
                                                                                                                              tcatCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGTTTGCCCCGCCACAAGGAATAACAATTGAAGCgg                 >  1:1141596/1‑151 (MQ=255)
                                                                                                                              tcatCTGTTGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGTTTGCCCCGCCACAAGGAATAACAATTGAAGCgg                 >  1:800004/1‑151 (MQ=255)
                                                                                                                                 tCTGTTGGATGAAATTGCGAAGGACCCGCCCATCTATCTGAGGGCAAAAAAAATGGCGGGAGAAGTTTTTCTAATAAAAAGGGACTGTTGATGGTTATGTAAGAGGGGGGGTTTTTTCCCCCGCCACAAGGAATAACAATTGAAGCg                  <  2:627977/147‑1 (MQ=255)
                                                                                                                                     ttGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGTTTGCCCCGCCACAAGGAATAACAATTGAAGCGGTAAACgg          >  2:1148141/1‑151 (MQ=255)
                                                                                                                                      tGGATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTAATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGTTTGCCCCGCCACAAGGAATAACAATTGAAGCGGTAAACg           >  1:187499/1‑149 (MQ=255)
                                                                                                                                        gATGAAATTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGTTTGCCCCGCCACAAGGAATAACAATTGAAGCGGTAAACGGAAt       >  1:245397/1‑151 (MQ=255)
                                                                                                                                             aaTTGCGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGTTTGCCCCGCCACAAGGAATAACAATTGAAGCGGTAAACGGAATGCTcg  >  2:1417362/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCTGGGGGTGATGACTTTTGGTGAGATTGTTATTGACGCCAGCCACACCGCCCTGTTGACCCGGGCATTCGCGCCACTGGCGGATGACGCGACGTCTGTGTGGCAGGCGCGTAGCATTCAATTCATTCATCTGTTGGATGAAATTGTGCAGGAACCGGCCATCTATCTGATGGCCAGAAAAATTGCGTGAGAAGGATTTCTCATTAATAAGGACTGTTGATGGCTATGTCAGACGAGGCGATGTTTGCCCCGCCACAAGGAATAACAATTGAAGCGGTAAACGGAATGCTCG  >  NC_000913/2611681‑2611972

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