Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,965 A→T 6.7% intergenic (+439/+101) narI → / ← rttR nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,9650AT 6.7% 687.2 / 4.5 255intergenic (+439/+101)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (93/126);  new base T (6/12);  total (106/149)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 259 total aligned reads displayed.

CTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGT  >  NC_000913/1286816‑1287114
                                                                                                                                                     |                                                                                                                                                     
    tCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                                                                                                                   >  2:8158975/1‑150 (MQ=255)
     ccTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                                                                                                                   >  1:4482773/1‑149 (MQ=255)
     ccTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                                                                                                                   <  2:4482773/149‑1 (MQ=255)
      cTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCg                                                                                                                                                 >  1:438646/1‑150 (MQ=255)
       gtCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                                                <  2:4779861/149‑1 (MQ=255)
         ccccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                                              <  1:3321668/150‑1 (MQ=255)
              ggtGGGCAGAATATTTTATTGCGGATTAGTTTGAGCATTCCGGGGGTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAc                                                                                                                                         <  1:5157075/148‑1 (MQ=14)
              cTTGGTCAGAATACTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAACGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAc                                                                                                                                         <  2:3165548/150‑1 (MQ=255)
              cTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAc                                                                                                                                         <  1:4660528/150‑1 (MQ=255)
               ttGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAcc                                                                                                                                        <  1:1684259/150‑1 (MQ=255)
                  ggCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCACACGCCCTGACCCACCAGCCCAACCATa                                                                                                                                     >  1:3100219/1‑150 (MQ=14)
                    cAGAATATTTGATTGCGGATTCGTTTTAGATTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGt                                                                                                                                   <  1:6009775/150‑1 (MQ=255)
                     aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCACTCGCTGACGCTAGAGCCGAACCGTAGTc                                                                                                                                  >  2:5746485/1‑150 (MQ=255)
                      gAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAGGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCg                                                                                                                                 <  2:736114/150‑1 (MQ=255)
                        atatTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGaa                                                                                                                               <  1:2955/150‑1 (MQ=255)
                         tatTTGAGGGCGGATGGGTTTGAGCATTCGGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAg                                                                                                                              <  2:6009602/150‑1 (MQ=14)
                             tGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTtc                                                                                                                          <  1:5912876/150‑1 (MQ=255)
                               aTTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTtctc                                                                                                                        <  1:3853072/150‑1 (MQ=255)
                                ttGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGGTGCGGGAACGCTCTGTCGCTGTCGGTCGTGTGGAAGCTTTGTCGATGGTTGTCa                                                                                                                       >  2:3116915/1‑150 (MQ=11)
                                 tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCATCGGGTCGTTGCCAGCGGCAAACCACACACGCAGAAGCACGATCCCAACCATATTCGAAGcgtaacaa                                                                                                                      >  1:6201254/1‑143 (MQ=14)
                                   cGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCATCGGGTCGTTGCCTGCGGCAACCCTCAAACGCAGACCCTCGAGCCCAACCAAAGACCAAGCTTCACAAcc                                                                                                                    >  1:5359965/1‑150 (MQ=14)
                                      agtCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTc                                                                                                                 <  1:7504714/148‑1 (MQ=255)
                                        tCGTTTGAGAATTCCGGGGCTTTTGACAGTGAGGGTGGTGGGGGTTGGTTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTccc                                                                                                               <  1:2060397/150‑1 (MQ=14)
                                        tCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTccc                                                                                                               <  2:4832163/150‑1 (MQ=255)
                                          gTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCCGCAACCCGCTCACCCTGACCCTCGAGCCGCACCATAGTCGAAGCTTCTCAACCGAccccc                                                                                                             >  1:8227768/1‑149 (MQ=14)
                                           tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGc                                                                                                            >  2:7938848/1‑150 (MQ=255)
                                            ttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCGCTCTCGCTGTCGCTCGAGGCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                                           >  2:783376/1‑150 (MQ=11)
                                              gagaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGACGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATg                                                                                                         >  2:7511391/1‑150 (MQ=255)
                                                 aaTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGc                                                                                                      <  2:7778760/150‑1 (MQ=255)
                                                          gCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                             <  2:5102885/150‑1 (MQ=38)
                                                            ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGa                                                                                           <  1:5302064/150‑1 (MQ=38)
                                                              ttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTACCCGCATGGGCAGGATATTTGAtt                                                                                         >  2:1788493/1‑150 (MQ=255)
                                                               aggaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTg                                                                                        <  1:758454/147‑1 (MQ=255)
                                                                 aacgTGATGGGGGTGGGGGACGGATTACTCCGCGCGGCGCGCTTCGCCCTTCGGTTCGTGGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCg                                                                                      <  2:6762546/147‑1 (MQ=14)
                                                                 aaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCGCGAGGCCTACCTGAGTCGAAGCTTCTCAACCGTACCCCCCAGGGCCGAAAATTTGATTGCg                                                                                      >  1:6664609/1‑150 (MQ=11)
                                                                     tGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAtt                                                                                  >  2:6813542/1‑150 (MQ=39)
                                                                      gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTc                                                                                 <  1:6166562/150‑1 (MQ=39)
                                                                      gATGGTGGGGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTc                                                                                 <  2:153965/150‑1 (MQ=255)
                                                                        tggtggtgtGTTAAGGATTGTTCAGCTCTGATCTCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGc                                                                               <  1:8041643/150‑1 (MQ=11)
                                                                         ggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCt                                                                              <  2:3975014/150‑1 (MQ=40)
                                                                          gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCtt                                                                             <  1:2825580/150‑1 (MQ=255)
                                                                             gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgag                                                                          <  2:4144099/150‑1 (MQ=255)
                                                                             gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgag                                                                          <  1:1065605/150‑1 (MQ=255)
                                                                             gtgctgtaaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCAAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgag                                                                          <  1:2243612/143‑1 (MQ=255)
                                                                              tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgaga                                                                         <  1:7169295/150‑1 (MQ=255)
                                                                                ggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGt                                                                       <  1:7747624/150‑1 (MQ=255)
                                                                                  ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCCTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTc                                                                     <  1:6927716/150‑1 (MQ=255)
                                                                                   gAAGGATTACTCAGCGCGGCGCGCTTCGCCCTGCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                                                    <  2:8227768/150‑1 (MQ=255)
                                                                                    aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                   <  1:422607/150‑1 (MQ=255)
                                                                                    aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                   <  2:905364/150‑1 (MQ=255)
                                                                                    aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                   <  1:6463846/150‑1 (MQ=255)
                                                                                     aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                                  <  1:1864772/150‑1 (MQ=255)
                                                                                     aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                                  <  1:5632688/150‑1 (MQ=255)
                                                                                        aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCTCTCGCGCGGACGCGCGAGCCGAACCTGAGTCGAAGCGTCTCATACGTCCCCGCCAGGGCAGGAGATATGATAGCGGATTCGCTTGGGAGttgcgggcc                                                               >  2:7618462/1‑143 (MQ=14)
                                                                                         ttACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACATTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGTGAGTTCAGGGACt                                                              <  1:3397843/150‑1 (MQ=255)
                                                                                            cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCGCGCTGACGCGCGAGCCGAACCGTAGTCGAAGCGTCTCATCCTTCCCCGCAAGGGCCGAAGATTTGATTGCGGATTCGCGGGGGAGTTCAGGGACtttt                                                           >  2:959818/1‑150 (MQ=14)
                                                                                               aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaa                                                        <  1:2362842/150‑1 (MQ=255)
                                                                                               aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAACTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGATTGGTAGTTCAGGGACTTTTGaa                                                        >  2:6182891/1‑150 (MQ=11)
                                                                                                ggggtGCGCTCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                       <  1:963880/146‑1 (MQ=255)
                                                                                                ctgcTGCGCTCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                       <  1:4558487/148‑1 (MQ=255)
                                                                                                 cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACttt                                                            <  1:477998/144‑1 (MQ=255)
                                                                                                   cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTg                                                    <  1:1286998/150‑1 (MQ=255)
                                                                                                     gcgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCCAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAT                                                  <  2:77582‑M1/150‑1 (MQ=255)
                                                                                                      cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGGCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTGGTGAGTTCAGGGACTTTTGAAAGTGATG                                                 >  2:4786247‑M1/1‑150 (MQ=255)
                                                                                                       gcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGGCGAACCGTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAATATTTGATTGCGGATTCGGTTGGGAGTTCAGGGACTTTTGTAAGTGATGG                                                >  2:7143948‑M1/1‑150 (MQ=255)
                                                                                                           ttCGCCCTTCGGGTCGTTGCCGGCCGCAACGCTCACACGCAGACGCACGAGCCGAACCATAGTCGAATCATCTCAACAATAACCCCAATGGCCGAAAATTTGATTTAGGATTCGCTTGAGAGTTAAGGGCATTTTGTCAGAGAtgggggg                                            >  2:2976174/1‑149 (MQ=12)
                                                                                                           gtcGCCCGTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGT                                            <  1:2558215‑M1/149‑1 (MQ=255)
                                                                                                            tCGCCCTTCGGGTAGTTGACTGCCGAACCGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTG                                           <  1:7563640‑M1/150‑1 (MQ=255)
                                                                                                             cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGG                                          <  2:6692731‑M1/150‑1 (MQ=255)
                                                                                                               cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGG                                        <  2:4389379‑M1/150‑1 (MQ=255)
                                                                                                                ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCATTGACTTTTGAAAGTGATGGTGGTGGGGG                                       <  1:3613279‑M1/150‑1 (MQ=255)
                                                                                                                  ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGGGATGGTGGTGGGGGAA                                     <  2:6046243‑M1/150‑1 (MQ=255)
                                                                                                                  ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAA                                     <  2:2948271‑M1/150‑1 (MQ=255)
                                                                                                                  ttCGGGTCGTTGCCTGCGGCAACGCTCTCACGCTGACGCTCGAGGCGAACCTTTGTAGAAGCTTATCATACATACCCGCATGGGCAGAATATTTTATTTCGGTTTAGGTTGGGAGTTAAGGGACTTTTGAAAGTGATTGTGGTGGGGGaa                                     >  2:6980518/1‑150 (MQ=11)
                                                                                                                    cGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGG                                   <  1:7755253‑M1/150‑1 (MQ=255)
                                                                                                                     gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGA                                  <  2:2464161‑M1/150‑1 (MQ=255)
                                                                                                                     gggTCGTTGCCTGCGGCAACGCTCTCTCGCGGACGCTCGAGCCCAACCTTAGTCGAAGCTTATCAACCTTCCCCGCATGGGCAGAATATTTGATTTCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGA                                  >  2:6351763‑M1/1‑150 (MQ=255)
                                                                                                                      ggTCGTTGCCTGCGGCAACGCTCTCACGCTGACGCACGAGACGAACCTTAGTCGAAGCTTCTCAAACTTCCCCCCATGGGCAGAATATTTAAATGCGGATTCGTTTGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAT                                 >  2:3441022‑M1/1‑150 (MQ=255)
                                                                                                                       agcgtTGCCGCAGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTACCCCGCATGGGCAGAATATTTGATTGCGGATTCGCGTGAGAGTTCAGGGACTTTTGAAAGGGATGGTGGTGGGGGAAGGATT                                >  2:518475‑M1/3‑150 (MQ=255)
                                                                                                                          gTTGCCCGCGGCAACGCTCTCTCGCTGACGCTCGCGCCGAACCGGCGTCGAAGCTTCTCCTCCGTCCCCGCATGGGCAGACTACTTGATTGCGGATTCGCTTGAGTGTTCAGGGACTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACt                             <  2:1712093/150‑1 (MQ=17)
                                                                                                                           ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTattt                            >  2:3487994/1‑147 (MQ=255)
                                                                                                                            tGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGTAGCTTCTCATCCTTCCCCGCATGGGCAGTATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCa                           <  1:5058417/150‑1 (MQ=255)
                                                                                                                              ccTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCaca                         >  2:6437332/1‑148 (MQ=255)
                                                                                                                               cTGCGGCAACGCTCTCTCGCTGGCGCTCAAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcg                        >  2:8115146/1‑150 (MQ=255)
                                                                                                                                tGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                       >  2:5354074/1‑150 (MQ=255)
                                                                                                                                    gCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcg                   >  2:4210565/1‑150 (MQ=255)
                                                                                                                                         gCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCGGCGCGCtt              >  2:2634156/1‑150 (MQ=255)
                                                                                                                                         gCTCTCTCGCCGACGCTCGCGCCGAACCCTAGTCGCAGCTTCTCACCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGAATTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCtt              <  2:7263916/150‑1 (MQ=17)
                                                                                                                                          ctctctCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc             <  1:947210/150‑1 (MQ=255)
                                                                                                                                          ctctctCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc             >  2:6072695/1‑150 (MQ=255)
                                                                                                                                          ctctctCGCTGACCCTCGAGCCGAACCATAGTCGAAGCGTCTCATCCATCACCCCAGGGGCAGAATATTTGATTTCGGATTCGATTGAGAGTTCAGGGACTTTTTAAATTGATGGTGGTGGGGGAAGGATTACTAACCGCGGCGCGCtta             >  1:3655191/1‑149 (MQ=17)
                                                                                                                                           tctctcGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg            <  1:5834981/150‑1 (MQ=255)
                                                                                                                                           tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg            <  2:1466317/150‑1 (MQ=255)
                                                                                                                                           tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg            <  1:8089250/150‑1 (MQ=255)
                                                                                                                                             ccgcGCTGCCGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGGTTTGGCAGTTTATTTGATTGCGGATTCGCTTGAGAGTTCAGGGTATTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGcc          <  1:2170824/147‑1 (MQ=14)
                                                                                                                                                 gCTGTCGCTCTAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc      <  2:541489/150‑1 (MQ=255)
                                                                                                                                                 gCTGTCGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc      <  2:1946180/150‑1 (MQ=255)
                                                                                                                                                 gCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc      >  2:5460080/1‑150 (MQ=255)
                                                                                                                                                 gCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGCGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc      <  1:7330340/150‑1 (MQ=255)
                                                                                                                                                 gCTGACGCTCGAGCCGAACCATAGTCGAATCTTCTCATCCTTCCCCGCACGCGCTGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc      <  2:2141524/150‑1 (MQ=17)
                                                                                                                                                     |                                                                                                                                                     
CTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGT  >  NC_000913/1286816‑1287114

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