Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,318,079 A→T 7.8% intergenic (‑73/+134) rimP ← / ← metY ribosome maturation factor for 30S subunits/tRNA‑Met

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,318,0790AT 7.8% 308.1 / 14.0 90intergenic (‑73/+134)rimP/metYribosome maturation factor for 30S subunits/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (42/41);  new base T (4/3);  total (46/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

AACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGCC  >  NC_000913/3317935‑3318226
                                                                                                                                                |                                                                                                                                                   
aaCCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATc                                                                                                                                                 >  4:937709/1‑149 (MQ=255)
 aCCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCg                                                                                                                                                <  3:12210/149‑1 (MQ=255)
      ttCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCt                                                                                                                                           >  3:1656479/1‑149 (MQ=255)
           aaCCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTtata                                                                                                                                      <  4:1656479/149‑1 (MQ=255)
                        aaCTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATaa                                                                                                                                   <  1:1140227/139‑1 (MQ=255)
                        aaCTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATaa                                                                                                                                   >  2:1140227/1‑139 (MQ=255)
                            ggCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAAc                                                                                                                     <  3:937709/149‑1 (MQ=255)
                                         tctcTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCggg                                                                                                                                              >  2:531952/1‑111 (MQ=255)
                                         tctcTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCggg                                                                                                                                              <  1:531952/111‑1 (MQ=255)
                                         tctcTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCt                                                                                                                          <  3:1569431/131‑1 (MQ=255)
                                         tctcTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCt                                                                                                                          >  4:1569431/1‑131 (MQ=255)
                                           tctGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAAc                                                                                                                     >  4:560346/1‑134 (MQ=255)
                                           tctGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAAc                                                                                                                     <  3:560346/134‑1 (MQ=255)
                                           tctGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTg                                                                                                      >  1:434821/1‑149 (MQ=255)
                                            ctGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTgg                                                                                                     <  2:434821/149‑1 (MQ=255)
                                               ttAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCt                                                                                                                          >  2:1494103/1‑125 (MQ=255)
                                               ttAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCt                                                                                                                          <  1:1494103/125‑1 (MQ=255)
                                                          tctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAAccc                                                                                                                           <  4:1041048/113‑1 (MQ=255)
                                                          tctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAAccc                                                                                                                           >  3:1041048/1‑113 (MQ=255)
                                                          tctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCgg                                                                                                          <  3:763040/130‑1 (MQ=255)
                                                          tctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCgg                                                                                                          >  4:1739746/1‑130 (MQ=255)
                                                          tctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCgg                                                                                                          <  3:1739746/130‑1 (MQ=255)
                                                          tctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCgg                                                                                                          >  4:763040/1‑130 (MQ=255)
                                                          tctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGaa                                                                                       >  4:417383/1‑149 (MQ=255)
                                                           ctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCggg                                                                                                                                              >  4:1556631/1‑93 (MQ=255)
                                                           ctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCggg                                                                                                                                              <  3:1556631/93‑1 (MQ=255)
                                                               tgtgGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCgg                                                                                                                                               >  1:733723/1‑88 (MQ=255)
                                                               tgtgGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCgg                                                                                                                                               <  2:733723/88‑1 (MQ=255)
                                                               tgtgGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTc                                                                                            <  1:1242637/139‑1 (MQ=255)
                                                               tgtgGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTc                                                                                            >  2:1242637/1‑139 (MQ=255)
                                                                 tgGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCg                                                                                                           <  3:46288/122‑1 (MQ=255)
                                                                 tgGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCg                                                                                                           >  4:46288/1‑122 (MQ=255)
                                                                  ggACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAGTTGCGGttata                                                                                                     <  4:1270625/127‑5 (MQ=255)
                                                                  ggACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAGTTGCGGttata                                                                                                     >  3:1270625/1‑123 (MQ=255)
                                                                        gCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCggggtt                                                                                                                                           >  2:223763/1‑81 (MQ=255)
                                                                        gCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCggggtt                                                                                                                                           <  1:223763/83‑3 (MQ=255)
                                                                          ccACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCaaa                                                                          >  3:610681/1‑146 (MQ=255)
                                                                          ccACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCaaa                                                                          <  4:610681/146‑1 (MQ=255)
                                                                            aCCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCa                                                                     <  3:417383/149‑1 (MQ=255)
                                                                               cccAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCttt                                                                                                                                         >  3:733866/1‑78 (MQ=255)
                                                                               cccAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCttt                                                                                                                                         <  4:733866/78‑1 (MQ=255)
                                                                                ccAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTta                                                                                                                                        >  4:59446/1‑78 (MQ=255)
                                                                                ccAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTta                                                                                                                                        <  3:59446/78‑1 (MQ=255)
                                                                                                    cTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAttt                                                                                 >  1:1783956/1‑113 (MQ=255)
                                                                                                    cTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAttt                                                                                 <  2:1783956/113‑1 (MQ=255)
                                                                                                    cTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCAc                                                                    >  3:1303718/1‑126 (MQ=255)
                                                                                                    cTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCAc                                                                    <  4:1303718/126‑1 (MQ=255)
                                                                                                      aaaGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTg                                                                                                      <  1:1787373/90‑1 (MQ=255)
                                                                                                      aaaGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTg                                                                                                      >  2:1787373/1‑90 (MQ=255)
                                                                                                           ccAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCa                                                                                          <  1:1502756/97‑1 (MQ=255)
                                                                                                           ccAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCa                                                                                          >  2:1502756/1‑97 (MQ=255)
                                                                                                           ccAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATAttt                                      >  3:1894935/1‑149 (MQ=255)
                                                                                                              gTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGaa                                   >  2:1304216/1‑149 (MQ=255)
                                                                                                               ttATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGaa                                                       >  4:771953/1‑128 (MQ=255)
                                                                                                               ttATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGaa                                                       <  3:771953/128‑1 (MQ=255)
                                                                                                                 aTTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCa                                                                     <  4:1749606/112‑1 (MQ=255)
                                                                                                                 aTTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCa                                                                     >  3:1749606/1‑112 (MQ=255)
                                                                                                                 aTTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGaaaaa                                <  4:1894935/149‑1 (MQ=255)
                                                                                                                          tCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATAttt                                      >  3:1422723/1‑134 (MQ=255)
                                                                                                                          tCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATAttt                                      <  4:1422723/134‑1 (MQ=255)
                                                                                                                             gATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCa                                                                     >  4:598586/1‑100 (MQ=255)
                                                                                                                             gATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCa                                                                     <  3:598586/100‑1 (MQ=255)
                                                                                                                             gATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTCCGAAGGCCGAAGTCTTCACAGTATATTTGaaaaa                                >  2:828586/1‑137 (MQ=255)
                                                                                                                             gATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGaaaaa                                <  1:828586/137‑1 (MQ=255)
                                                                                                                               tAACAAAAAACCCCGATAAATCGGGGCTTTACATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGa                                    >  3:135056/1‑131 (MQ=255)
                                                                                                                               tAACAAAAAACCCCGATAAATCGGGGCTTTACATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGa                                    <  4:135056/131‑1 (MQ=255)
                                                                                                                                   aaaaaaCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTa                                                                  >  1:1812929/1‑97 (MQ=255)
                                                                                                                                   aaaaaaCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTa                                                                  <  2:1812929/97‑1 (MQ=255)
                                                                                                                                   aaaaaaCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGaa                                                       <  1:1250855/108‑1 (MQ=255)
                                                                                                                                   aaaaaaCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGaa                                                       >  2:1250855/1‑108 (MQ=255)
                                                                                                                                    aaaaaCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTgg             >  4:262027/1‑149 (MQ=255)
                                                                                                                                     aaaaCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGt            >  4:1434074/1‑149 (MQ=255)
                                                                                                                                      aagccccGATTTATCGGGGTTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGaga                 >  1:560790/4‑141 (MQ=255)
                                                                                                                                      aagccccGATTTATCGGGGTTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGaga                 <  2:560790/140‑3 (MQ=255)
                                                                                                                                      aaaCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGa                             <  3:1412927/131‑1 (MQ=255)
                                                                                                                                      aaaCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGa                             <  3:208431/131‑1 (MQ=255)
                                                                                                                                      aaaCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGa                             >  4:208431/1‑131 (MQ=255)
                                                                                                                                      aaaCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGa                             >  4:1412927/1‑131 (MQ=255)
                                                                                                                                       aaCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGaa                                                       >  4:1827436/1‑104 (MQ=255)
                                                                                                                                       aaCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGaa                                                       <  3:1827436/104‑1 (MQ=255)
                                                                                                                                        gccccGATTTATCGGGGTTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAg                                             >  2:1696753/2‑113 (MQ=255)
                                                                                                                                        gccccGATTTATCGGGGTTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAg                                             <  1:1696753/112‑1 (MQ=255)
                                                                                                                                        aCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGaa                                                                                       >  2:1770594/1‑71 (MQ=255)
                                                                                                                                        aCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGaa                                                                                       <  1:1770594/71‑1 (MQ=255)
                                                                                                                                          cccGATTTATCGGGGTTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTa                                                                  >  2:1199657/1‑90 (MQ=255)
                                                                                                                                          cccGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACttt                                                                                             >  4:711310/1‑63 (MQ=255)
                                                                                                                                          cccGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACttt                                                                                             <  3:711310/63‑1 (MQ=255)
                                                                                                                                           ccGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCggg      >  3:1216505/1‑149 (MQ=255)
                                                                                                                                              aTAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGc   <  3:1434074/149‑1 (MQ=255)
                                                                                                                                               tAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGcc  <  3:262027/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGCC  >  NC_000913/3317935‑3318226

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