Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,801,056 A→T 11.5% N31I (AAT→ATT)  nrdI → NrdEF cluster assembly flavodoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,801,0560AT11.5% 135.6 / 8.5 60N31I (AAT→ATT) nrdINrdEF cluster assembly flavodoxin
Reads supporting (aligned to +/- strand):  ref base A (34/18);  new base T (3/4);  total (37/23)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

TGGTTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTGC  >  NC_000913/2800908‑2801203
                                                                                                                                                    |                                                                                                                                                   
tggtTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGGCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCCGCGTTTTATCGAGCGGTTAGGTCTGCCAGGGGAGCGCATCCCGCTCaa                                                                                                                                                     >  2:463385/1‑149 (MQ=255)
  gtTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCTTGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGGTATATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATg                                                                                                                                                   >  2:695032/1‑149 (MQ=255)
     tccgtccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCGCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGc                                                                                                                                                <  2:943766/149‑1 (MQ=255)
         tccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGa                                                                                                                                            <  2:556357/149‑1 (MQ=255)
         tccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGGTTAGGTCTGCCCGCGGTGCGCATCCCGCTAAATGAGCGGGa                                                                                                                                            >  1:585574/1‑149 (MQ=255)
          ccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGat                                                                                                                                           <  1:15245/149‑2 (MQ=255)
           cgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGGCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAAc                                                                                                                                          <  2:381800/149‑1 (MQ=255)
            ggACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGGTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGatgc                                                                                                                                         <  2:425942/149‑4 (MQ=255)
               cATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTCCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGGCTGCCCGCGGTGCGCAGCCCGCTCAATGAGCGGGAACGGAt                                                                                                                                      <  1:1012168/149‑1 (MQ=255)
                aTGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGAtt                                                                                                                                     <  1:407667/149‑1 (MQ=255)
                 tGATTAACTGTCTGCATCCAGCGCCACACGCGGCCAGTGCTTGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTc                                                                                                                                    <  2:569043/149‑1 (MQ=255)
                  gATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCa                                                                                                                                   >  1:1019618/1‑149 (MQ=255)
                   aTTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAACTCGTCTACTTCTCCAGCAGCTCCGAACACACGCAGCGTTTTATCGAACGTTTAGGGCTGCCCGCGGGGCGCATCCCGCTCAATGAGCGGGAACGGATTCAg                                                                                                                                  <  2:772198/149‑1 (MQ=255)
                        ccGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGGCTGCCCGCGGTGCGCATCCCGCTCAATGCGCGGGAACGGAGTCAGGTAgc                                                                                                                             >  1:570619/1‑148 (MQ=255)
                        ccGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATAGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGa                                                                                                                             >  2:648912/1‑149 (MQ=255)
                                aTCCAGCGCCACACGCGGCCAGTGCATGCGCCAGCTCGTCTACTTCTCCCGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGAACGGATTCAGGTAGACGAGCCtt                                                                                                                     >  1:339988/1‑149 (MQ=255)
                                 tCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCGCAATGAGCGGGAACGGATTCAGGTAGACGAGCCtga                                                                                                                    >  1:328621/1‑147 (MQ=255)
                                  ccAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTAc                                                                                                                   <  1:556006/149‑1 (MQ=255)
                                    aGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGCACGTTTAGGTCTGCCCGCGGTGCGCCTCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACAt                                                                                                                 >  1:931260/1‑149 (MQ=255)
                                    aGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGt                                                                                                                                <  1:784747/134‑1 (MQ=255)
                                    aGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGt                                                                                                                                >  2:784747/1‑134 (MQ=255)
                                     gcgcCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATc                                                                                                                <  2:892160/149‑1 (MQ=255)
                                     gcgcCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACCCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGCCGCGCCTTACATc                                                                                                                >  2:337703/1‑149 (MQ=255)
                                          acacGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGAt                                                                                                           >  1:834022/1‑149 (MQ=255)
                                           cacGCGGCCAGTGCTGGAGCGAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATc                                                                                                          <  1:591359/149‑1 (MQ=255)
                                           cacGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATc                                                                                                          >  2:39285/1‑149 (MQ=255)
                                              gcgGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGat                                                                                                                                           <  1:366219/113‑2 (MQ=255)
                                              gcgGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGat                                                                                                                                           >  2:366219/1‑112 (MQ=255)
                                                 gCCTGTGCATGAGCCAGATCGTCAAGATCTCCAGGAGATCCGAAAAGACGGAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCTATGATCGGGAACTGATTCAGGTAGACTAGCCTTACATCATGATCGTGccc                                                                                                    <  2:163555/149‑1 (MQ=255)
                                                     gTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTCCCCTCtt                                                                                                <  2:80646/149‑1 (MQ=255)
                                                        cATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTAcg                                                                                             >  2:5485/1‑149 (MQ=255)
                                                                   tCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGt                                                                                  <  1:379106/149‑1 (MQ=255)
                                                                   tCGTCTACTTCTCAAGCAGCTCCGATAAGACGCAGCGTTTTATCGAACGTTTAGGTCTGCCAGGGGTGCGGATACCGCTGAATGAGGGGGTACGCATTCAGGTAGACGAGCCGTACCTCCTGCGCGTGCGCGCGGAcggcggcggcggg                                                                                  >  2:562511/1‑148 (MQ=255)
                                                                     gTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACCTTTAGGTCTGCCCGAGGTGCCCATCCCGCTCAATGAGCGGGAATGGATTCAGGTAGCCGAGCCTTACATCCTGATCGAGACCACTTACCGCGGCTGCGGTAc                                                                                >  1:462841/1‑149 (MQ=255)
                                                                           ttCTCCAGCAGCTCCGAAACCACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTgg                                                                          >  1:179916/1‑149 (MQ=255)
                                                                           ttCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCGGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTgg                                                                          >  2:949353/1‑149 (MQ=255)
                                                                                 agcagcTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGt                                                                    >  2:431058/1‑149 (MQ=255)
                                                                                  gcagcTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTg                                                                   >  1:405435/1‑149 (MQ=255)
                                                                                     gcTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGccc                                                                >  2:475620/1‑148 (MQ=255)
                                                                                      cTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTTCCCGCGGTGCGCATCCCGATCAATGAGCGGGAACGTATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCAc                                                               >  1:363800/1‑149 (MQ=255)
                                                                                        ccGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGCGGGGGAACGGATTCAGGTAGACGAGCCTTACTTCCGGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCGACgc                                                             >  2:701023/1‑148 (MQ=255)
                                                                                             aaCACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGt                                                        >  2:834734/1‑149 (MQ=255)
                                                                                                     gCGTTTTATCGATCGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGAGGCGGTACGGCTGGCGCGGTGCAACGACAGGTAATTCGCt                                                <  2:1008296/149‑1 (MQ=255)
                                                                                                     gCGGTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCGCACTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCt                                                <  2:163410/149‑1 (MQ=255)
                                                                                                        ttttATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCGGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCtttt                                             >  1:838914/1‑149 (MQ=255)
                                                                                                           tATCGAACGTTTAGGGCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTaa                                          <  1:362630/149‑1 (MQ=255)
                                                                                                             tCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCCGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAac                                        >  1:659418/1‑149 (MQ=255)
                                                                                                                    tttAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGc                                 >  1:791832/1‑149 (MQ=255)
                                                                                                                    tttAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGc                                 <  1:695032/149‑1 (MQ=255)
                                                                                                                         gTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCACATACGGCGGCGGCGGTTCGGCTGGCGCGGTGCCCAGACAGGTAATTCGCTTTTTAAACGACGATCACAAc                            >  2:670969/1‑149 (MQ=255)
                                                                                                                         gTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGGAGACGCGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAACCGACGAGCACAAc                            >  1:207759/1‑149 (MQ=255)
                                                                                                                            tGCCCGCGGGGCGCAGCCCGCGCAATGAGCGGGAACGGAGTCAGGTAGACGAGCCGTACATCCTGAGCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCgg                         <  2:570619/149‑1 (MQ=255)
                                                                                                                             gCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCggg                        >  1:121337/1‑149 (MQ=255)
                                                                                                                               ccGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGCGCCTAACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGGGCCACGAAAGGTAATTCGCTTTTGAAACGTCGAGCACAACCGGGCg                      >  2:667047/1‑149 (MQ=255)
                                                                                                                                cgcgGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCTTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGCCAGGTAATTCGCTTATTAAACGACGAGCCCAACCGGGCGt                     >  2:646562/1‑149 (MQ=255)
                                                                                                                                         gttCCCGCTCATTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCTGCGGCGCTACGCCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTcgag            >  1:484616/3‑147 (MQ=255)
                                                                                                                                          aTCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCgcgg           >  2:1004835/1‑149 (MQ=255)
                                                                                                                                            cccGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTTCATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTTAACGACGAGCACAACCGGGCGTTGCTTCgcggcg         >  1:959631/1‑149 (MQ=255)
                                                                                                                                            cccGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCgcggcg         >  1:776524/1‑149 (MQ=255)
                                                                                                                                             ccGCTCAATGAGCGGGAACGGATTCAGGTAGCCGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGt        >  2:333017/1‑149 (MQ=255)
                                                                                                                                             ccGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCGTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGt        <  1:39285/149‑1 (MQ=255)
                                                                                                                                                   aaTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTGc  >  1:277293/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGGTTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTGC  >  NC_000913/2800908‑2801203

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