Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,801,056 A→T 8.4% 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,0560AT 8.4% 330.4 / 13.9 126N31I (AAT→ATT) nrdINrdEF cluster assembly flavodoxin
Reads supporting (aligned to +/- strand):  ref base A (62/52);  new base T (3/8);  total (65/61)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 128 total aligned reads displayed.

CTGGTTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTGCT  >  NC_000913/2800907‑2801204
                                                                                                                                                     |                                                                                                                                                    
 tggtTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAAt                                                                                                                                                     >  2:1626492/1‑150 (MQ=255)
  ggtTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTg                                                                                                                                                    <  1:462278/150‑1 (MQ=255)
  ggtTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATg                                                                                                                                                    <  1:1328211/150‑1 (MQ=255)
  ggtTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATg                                                                                                                                                    >  1:1618019/1‑150 (MQ=255)
       ccgtccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCgg                                                                                                                                               <  1:1630699/150‑1 (MQ=255)
       ccgtccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCgg                                                                                                                                               >  1:872746/1‑150 (MQ=255)
       ccgtccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCgg                                                                                                                                               <  1:703179/150‑1 (MQ=255)
        cgtccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCggg                                                                                                                                              <  1:567550/150‑1 (MQ=255)
          tccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGat                                                                                                                                            <  2:1677480/150‑2 (MQ=255)
          tccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGaa                                                                                                                                            <  1:601448/150‑1 (MQ=255)
            cgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGatgc                                                                                                                                          <  2:22780/150‑4 (MQ=255)
            cgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACg                                                                                                                                          >  1:467654/1‑150 (MQ=255)
            cgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACg                                                                                                                                          <  1:1219018/150‑1 (MQ=255)
              gACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCCGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGatgcgc                                                                                                                                        <  1:204949/150‑6 (MQ=255)
              gACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGa                                                                                                                                        >  1:1208407/1‑150 (MQ=255)
              gACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGa                                                                                                                                        >  2:654904/1‑150 (MQ=255)
                  tGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGCTTCa                                                                                                                                    >  2:141847/1‑150 (MQ=255)
                   gATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAg                                                                                                                                   >  1:1491354/1‑150 (MQ=255)
                        aCCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAgc                                                                                                                              >  1:1565840/1‑149 (MQ=255)
                          cGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACg                                                                                                                            >  1:206638/1‑150 (MQ=255)
                           gTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAACCACGCAGCGTGTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCCAGGCGCGGGAACGGATTCAGGTAGACgg                                                                                                                           >  1:69407/1‑149 (MQ=255)
                             cTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGATCAATGAGCGGGAACGGATTCAGGTAGACGAGc                                                                                                                         >  1:1118611/1‑150 (MQ=255)
                              tccATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACACAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGcc                                                                                                                        <  2:1672201/148‑1 (MQ=255)
                              tGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGcc                                                                                                                        >  2:124564/1‑150 (MQ=255)
                               gCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCt                                                                                                                       >  2:575764/1‑150 (MQ=255)
                               gCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGGCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCt                                                                                                                       <  2:1355895/150‑1 (MQ=255)
                                  tCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTAc                                                                                                                    >  1:976496/1‑150 (MQ=255)
                                      gcgcCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATcc                                                                                                                >  2:1517097/1‑150 (MQ=255)
                                        gcCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTg                                                                                                              <  1:963705/150‑1 (MQ=255)
                                         ccACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGa                                                                                                             >  2:1374794/1‑150 (MQ=255)
                                            cacGCTGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCg                                                                                                          <  1:463270/150‑1 (MQ=255)
                                            cacGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCg                                                                                                          >  1:13448/1‑150 (MQ=255)
                                            cacGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCg                                                                                                          >  1:713808/1‑150 (MQ=255)
                                            cacGCGGCCAGGGCATGAGCCAGCTCGGCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGGCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCg                                                                                                          <  2:67831/150‑1 (MQ=255)
                                              cgcgGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTg                                                                                                        >  1:473146/1‑150 (MQ=255)
                                              cgcgGCCAGTACATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTg                                                                                                        >  2:1611636/1‑150 (MQ=255)
                                               gcgGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGc                                                                                                       <  2:615690/150‑1 (MQ=255)
                                                 ggCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGccc                                                                                                     <  2:1535510/150‑1 (MQ=255)
                                                      gTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTa                                                                                                >  1:1469854/1‑150 (MQ=255)
                                                         cATGAGCCAGGTCGTCTACTTCTCCAGCAGATCCGAAAAGACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATACTGATCGTGCCCTCTTAcgg                                                                                             <  2:272691/150‑1 (MQ=255)
                                                         cATGAGCCAGCTCGTATACTTCTACAGCAGCTCCGAAAACACGGAGCGTTTTATCGCACGTTTCGGTCTGCCCCCGGAGCGCATCCCGCTCAATGAGCGGGTACGGATTCACGTAGACGAGCCTTACATCATGATCTTGCCCTCTTAcgg                                                                                             >  1:272691/1‑150 (MQ=255)
                                                             aGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTAcggcggc                                                                                         <  2:1208407/150‑1 (MQ=255)
                                                              gCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTAcggcggcg                                                                                        >  1:383698/1‑150 (MQ=255)
                                                              gCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTAcggcggcg                                                                                        >  1:1417054/1‑150 (MQ=255)
                                                                   cTCGTCTACATCTCCAGCAGCTCCGAAATCACGCATCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGt                                                                                   <  2:72506/150‑1 (MQ=255)
                                                                      gTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACg                                                                                <  1:1380457/150‑1 (MQ=255)
                                                                      gTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTACCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACg                                                                                >  2:1436443/1‑150 (MQ=255)
                                                                       tCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACgg                                                                               >  1:239573/1‑150 (MQ=255)
                                                                        cTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGc                                                                              <  1:1198758/150‑1 (MQ=255)
                                                                                ccAGCAGCTCCGAAAACCCGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGAGCgg                                                                      >  2:911133/1‑150 (MQ=255)
                                                                                ccAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCgg                                                                      <  2:396096/150‑1 (MQ=255)
                                                                                  agcagcTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTg                                                                    <  2:206638/150‑1 (MQ=255)
                                                                                    cagcTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGcc                                                                  >  1:1676342/1‑150 (MQ=255)
                                                                                      gcTCGGATAACTCTCAGGGTTTTATCGAACGTTTAGGTCTTCCCGCTGTGGGATTCCCGCTCAATGAGCGGGAACGGTTTCAGGTAGACGAGCATTAAATCCTGATCGTGCCCTCTTAGGGGTGCGGGGGTACTTATGGCGCGGTgaaac                                                                <  2:1719770/150‑5 (MQ=255)
                                                                                      gcTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCAc                                                                <  2:13448/150‑1 (MQ=255)
                                                                                      gcTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCAc                                                                >  2:855770/1‑150 (MQ=255)
                                                                                       cTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCTCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACg                                                               >  2:511651/1‑150 (MQ=255)
                                                                                       cTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACg                                                               <  1:1018271/150‑1 (MQ=255)
                                                                                        tCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACCGCGCCAGCCGTACCGCCGCCGCCGtaagaggg                                                              <  2:467654/150‑8 (MQ=255)
                                                                                              aaCACGCAGCGTTTTATCGACCGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGCTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTa                                                        <  1:286268/150‑1 (MQ=255)
                                                                                                 acGCAGCTTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGAACGGATTGTGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAAtt                                                     <  1:1597736/150‑1 (MQ=255)
                                                                                                  cGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGCTACGGCTGGCGCGGTGCCACGACAGGTAATTc                                                    >  1:401884/1‑150 (MQ=255)
                                                                                                   gCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCg                                                   <  1:855770/150‑1 (MQ=255)
                                                                                                   gCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCg                                                   >  2:1429241/1‑150 (MQ=255)
                                                                                                        gTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACTGATTCAGGTAGACGAGCATTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCtttt                                              >  1:330794/1‑150 (MQ=255)
                                                                                                        gTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCtttt                                              >  2:1204149/1‑150 (MQ=255)
                                                                                                        gTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCtttt                                              <  2:397049/150‑1 (MQ=255)
                                                                                                         ttttATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCttttt                                             <  2:1491354/150‑1 (MQ=255)
                                                                                                         ttttATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCttttt                                             <  1:172008/150‑1 (MQ=255)
                                                                                                          tttATCGAACGCTTCGGTCAGCCCGCGGCGCGCCTCCCGCGCAATGACCGGGAACGGACTCCGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTa                                            <  2:218225/150‑1 (MQ=255)
                                                                                                            tATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTaaa                                          >  2:508060/1‑150 (MQ=255)
                                                                                                            tATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTaaa                                          >  1:1357677/1‑150 (MQ=255)
                                                                                                            tATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTaaa                                          <  2:1676342/150‑1 (MQ=255)
                                                                                                             aTCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTATGTCGCTTTTTAAac                                         >  1:29486/1‑150 (MQ=255)
                                                                                                             aTCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAac                                         <  1:780092/150‑1 (MQ=255)
                                                                                                             aTCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAac                                         >  1:965793/1‑150 (MQ=255)
                                                                                                             aTCGAACGTTTAGGGCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAac                                         <  1:1380215/150‑1 (MQ=255)
                                                                                                               cGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGAACGGCTGGCGCTGTGCCACGACAGGTAATTCGATTTTTAAacga                                       >  2:1098061/1‑150 (MQ=255)
                                                                                                                    gTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGCGCCACGACAGGTAATTCGCTTTTTAAACGACGAGc                                  <  2:1618019/150‑1 (MQ=255)
                                                                                                                         ggTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAAc                             >  2:24759/1‑150 (MQ=255)
                                                                                                                              gCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCCCGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGc                        >  1:713818/1‑150 (MQ=255)
                                                                                                                              gCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGc                        <  2:971236/150‑1 (MQ=255)
                                                                                                                              gCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGc                        >  2:1703271/1‑150 (MQ=255)
                                                                                                                                    ggTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCt                  <  2:1469854/150‑1 (MQ=255)
                                                                                                                                    ggTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCt                  <  1:654904/150‑1 (MQ=255)
                                                                                                                                     gTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCtt                 <  2:1526303/150‑1 (MQ=255)
                                                                                                                                     gTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCtt                 >  1:1526303/1‑150 (MQ=255)
                                                                                                                                     gTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTAGCTTTTTAAACGACGAGCACAACCGGGCGTTGCtt                 >  2:347070/1‑150 (MQ=255)
                                                                                                                                        cgcATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCGCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTcgc              >  1:79632/1‑150 (MQ=255)
                                                                                                                                         cgttCCCGCTCATTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTcgcg             >  2:425468/4‑150 (MQ=255)
                                                                                                                                          gttCCCGCTCATTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGGGGTACGGCTGGCGCGGTGCCCCGACAGGTAGTTCGATTTTTAAACGACGCGCACAACCGGGCGTTGCTCCgcgg            >  1:1096699/3‑150 (MQ=255)
                                                                                                                                          cATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCgcgg            <  1:380276/150‑1 (MQ=255)
                                                                                                                                          cATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGAGCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCgcgg            <  1:141847/150‑1 (MQ=255)
                                                                                                                                             cccGCTCATTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGt         >  2:1076926/1‑150 (MQ=255)
                                                                                                                                              ccGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGtt        >  1:1476378/1‑150 (MQ=255)
                                                                                                                                              cagatCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCTGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGtt        <  2:1118611/146‑1 (MQ=255)
                                                                                                                                               cGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATGGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTa       >  1:1244662/1‑150 (MQ=255)
                                                                                                                                                gCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTAt      >  1:107091/1‑150 (MQ=255)
                                                                                                                                                  tCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCTTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTg    <  1:24759/150‑1 (MQ=255)
                                                                                                                                                    aaTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTGCt  >  1:967662/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGGTTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTGCT  >  NC_000913/2800907‑2801204

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