Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,965 A→G 7.6% 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,9650AG 7.6% 211.1 / 6.2 82intergenic (+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 (35/38);  new base G (4/3);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

TCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTA  >  NC_000913/1286731‑1287202
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              
tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCt                                                                                                                                                                                                                                           <  1:834573/239‑1 (MQ=255)
tCAGCGCTGCCCGCTTCGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCTCTCGCTGGCACTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCt                                                                                                                                                                                                                                           <  1:93345/239‑1 (MQ=255)
 cAGCGCTGCGCGCTTCGCCCTTCGTGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTc                                                                                                                                                                                                                                          <  1:202853/239‑1 (MQ=255)
  aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCGGACGCTc                                                                                                                                                                                                                                          >  1:53372/1‑238 (MQ=255)
     gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                                                                                                                                      <  1:680407/239‑1 (MQ=255)
       tgggcgcTTCGCCTGTCGGGCTGTTGCTTCCGGCAACGCTCCCGAGCTGTAGCTGAGGCAGATCCGTGGGGGTAGCTTCTCACCGTTCACCGCTTGGGCAGAATATTTGATTGCGGATTAGTTTGAGAATTCCGGTGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCAAGTCg                                                                                                                                                                                                                                    <  2:42228/236‑1 (MQ=11)
       tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                                                                                                                                                                                                                                     <  1:904363/238‑1 (MQ=255)
         cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCGTGGGCAGAATATTTGATTGCGGATTCGTTTGGGAGTTCCGGGACCTGTGTCTGATATACTGCTGGCGGTGGGGTCAAACAGCCCGGCGAACAGCGCCGTTCGGATCTTTGCCGGCGGCACCGCTCTCTCGCTGGAGCCCGAGCCGac                                                                                                                                                                                                                                  >  2:126374/1‑238 (MQ=2)
                   cgttCGGGGTGTTGCCGGCTGTACCGCTCTCTCGCTGGCGCTCAAGTCGAACCTGGGTCGAAGCTTCTCATCCTTCCCCGCTTTGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGa                                                                                                                                                                                                                        <  2:321779/237‑1 (MQ=14)
                   ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGa                                                                                                                                                                                                                        <  1:21805/239‑1 (MQ=255)
                   ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGa                                                                                                                                                                                                                        <  1:334028/239‑1 (MQ=255)
                    cTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCTCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGaa                                                                                                                                                                                                                       <  1:1148742/239‑1 (MQ=255)
                          gTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTtc                                                                                                                                                                                                                  <  1:920788/238‑1 (MQ=255)
                             gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGCGCCGAACCTTCGTCGTAGCTTCTCa                                                                                                                                                                                                               >  1:345560/1‑238 (MQ=255)
                             gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAt                                                                                                                                                                                                              <  1:1292670/239‑1 (MQ=255)
                               tGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGCACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGACTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCCCTCGCTGACGCCCGAGCCGAACCTTAGTCGAAGCTTCTCATc                                                                                                                                                                                                             >  1:639412/1‑238 (MQ=255)
                                    gCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCCTAAGCGCTGAGCGGTTGGCGCGTCGGGTCGTGGCGTCCGCCTACGCGCTCTCGCGGACGACACCAGCGAGCCTTAGCAGATTGTTCTCATCCTTccc                                                                                                                                                                                                       >  2:540775/1‑239 (MQ=14)
                                            gCTTTCTCTCTGGCTCTCAAGCCGACCCGTGGTGGAATCGTCTCACCCTTACCCGCTTGGGCAGAATGTTTGATTGCGGATTCGTTTGAGGATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                                                                                                               <  2:837814/239‑1 (MQ=14)
                                            gCTCTCTCGCTGGGGCTCAGGTCTAACCTGGGGCGAAGGTTTTCATCCTTCCCCGGTTGGGCAGAAGATTTGAGTGCGGGTTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                                                                                                               <  2:1237616/239‑1 (MQ=14)
                                            gCTCTCTCGCTGGCGCGCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                                                                                                               <  1:962397/239‑1 (MQ=255)
                                              tctctcGCTGGCGCTCAAGTCGAACCTTGGTCGAAGTTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                                                                                                                             <  1:70344/239‑1 (MQ=255)
                                              tctcgcGCTGGCGCTCAAGTCGAACCTTGGTCGAAGTTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATGCGTGTGAGAATGCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                                                                                                                             <  1:622273/239‑1 (MQ=255)
                                                    gCTGGCGCTCAAGTCGAGGCGTGGTGGAAGCTTCTCGTCCTTCGCCGTTTGGGCAGAATGGTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAata                                                                                                                                                                                        <  2:1158942/238‑1 (MQ=14)
                                                          gcTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCCCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCGGAATATTTGAtt                                                                                                                                                                                 >  1:587823/1‑239 (MQ=255)
                                                                ttCGACCCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAt                                                                                                                                                                           <  1:858264/238‑1 (MQ=255)
                                                                  cGACCCTGGGTCGAAGTTTTTCCTCCTTCCACGTTTGGGCAGAATATTTGATTGGGGATGCGTTTGAGAATTCGGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAtt                                                                                                                                                                          <  2:1110869/238‑1 (MQ=14)
                                                                        ttGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                                                                                                                   <  1:237469/239‑1 (MQ=255)
                                                                                        aTCGTTCCCCGCAGGGGCATAAGATTTAATTGCGGAGTGGTTTGAGAATTCCGGTGCTTATGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                                                                                                                   <  2:42469/239‑1 (MQ=14)
                                                                                            ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTGGCCCTTCGGGTCGAGTCCTGCGGCAACGCTCTCTGGCTGACGCTCGAGCCGATCCTTGGTAGAAGCTTGTCATCCTTCCCCGCATTGGCCGAATATTTAATTACGGGTTCGTTGGGGAGTGCAGGGGCTGTTGgaa                                                                                                                                               >  2:754684/1‑239 (MQ=11)
                                                                                             tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGGGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGTCATCTCTGCGTCCTTCGCCCTTCGGTTCGTGGCCTGGCGCAATGGCCTCTCGCTGACGCTCGAGCCGAACCTTCAACGTATGTTCTCATGCTTCCCCGTAAGGGCGGGGTAGGTTGTTGCTGATTCGCCTCTGATTTCAGCGACgtttagcag                                                                                                                                              >  2:5125/1‑234 (MQ=9)
                                                                                             tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGCACCTTAGGCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAATCGTTTGAGAATTCAGGGCATTTTGgaat                                                                                                                                              >  2:236383/1‑238 (MQ=14)
                                                                                             tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                                                                                                                   >  2:197588/1‑234 (MQ=255)
                                                                                             tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                                                                                                                   <  1:197588/234‑1 (MQ=255)
                                                                                              ccccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTTGCTCGGGAGTTCAGTGCCTTTTGAAAGt                                                                                                                                             >  1:231435/1‑239 (MQ=255)
                                                                                                   cTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATCTTTGATTGCGACTTCGCTTGGGAGTTCAGGGACTTTTGAAAGTGAtgg                                                                                                                                        >  1:593271/1‑239 (MQ=255)
                                                                                                     tGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGATTGAGAGTTCAGGGACTTTTGAAAGTGAtggtg                                                                                                                                      >  1:1092500/1‑239 (MQ=255)
                                                                                                              gatttgatttgGGCGTCGTTTGAGATTTTCATGATTTGTGTAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                                                               <  2:997307/236‑1 (MQ=14)
                                                                                                                tttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTgag                                                                                                                                                                  <  1:1176690/200‑1 (MQ=255)
                                                                                                                 ttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGACTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                                                                                          <  1:501165/239‑1 (MQ=255)
                                                                                                                 ttGATTGCGGATTCGTTAGAGAATTCAGGGGCTTTTGAAAGTGATGGGGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                                                                                          <  1:742128/239‑1 (MQ=255)
                                                                                                                 tagATTGCGTCGTCGTTCGAGTCTTCGGGGGTTTTGGGAGTGGATGGTGGTGGGGGAAGGATTTGCCTGCGCGGCGCGCTTCTCCCTTCGGGTCGTTGCCTGCGGCGACGCTCTCGCGCTGACGCTCAAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGTTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                                                                                          <  2:610104/237‑1 (MQ=11)
                                                                                                                      tGCGGATTCGTTTGAGAATTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTAGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGCACCTTAGCCGAAGCTTCCCATCCCTCCCCGCATTGGCAGAATATTTGTTTGGGGACTCGGTTCAGCGTTCAGGGCCTTTTGAGAGCGATGGTGGTGGGGGCAGGAGTAc                                                                                                                      >  2:858186/1‑238 (MQ=14)
                                                                                                                      tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGTTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCCCGAGCCGAACCTTAGTCGTAGCTTCTCATCCTTCCCCGCATGGGCAGTATATTTGTGTGCGGGTTGGCTAGGAAGGTCAGGGAGTTTTGAGAGTTATGGCGGTGGGGGAAGGATTaa                                                                                                                      >  2:1056210/1‑237 (MQ=14)
                                                                                                                      tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTTGAAGCTTCTCATCCTTCCCCGCATGGGCAGGATATTTGATTGCGGATTCGCTTGAGCGTTCAGGGGCTTTTGAAAGGGAtgggggtggggggtggtttgg                                                                                                                      >  2:364363/1‑229 (MQ=14)
                                                                                                                        cGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCTTTCGGGTCGTTGCATGTGGCAACGCTCTCTCGCTGACGCTCGAGGCGAACCTTAGTCGAAGCTTCTCATGCTTGCCCGCATGGGCAGGATATTTGATTGCGGATTCGTTTGAGAGTTCAGGGACTTTTGAAAGTGATGGAGGTGGGGGAAGGAttggtcc                                                                                                                   >  2:512704/1‑234 (MQ=14)
                                                                                                                        cGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGCACCTTAGTCGACTCTTCTCATCCTTCCCCGCATGGGCAGAGTATTTACTTTCGGACTCGCTTGCCAGCTCAGAGCTTTTTGATAGTGCTGGTGGTGTGGGAAAGAtgccttt                                                                                                                   >  2:732061/1‑233 (MQ=14)
                                                                                                                        cGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGTTTCGCTTGAGAGTTCAGGGACTTTTGAAAGAGATGGTGGTGGGGGAAGGATTACTCa                                                                                                                   >  1:408097/1‑239 (MQ=255)
                                                                                                                        cGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCa                                                                                                                   >  2:1336240/1‑239 (MQ=255)
                                                                                                                          gATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGAGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCgggg                                                                                                                                                        >  2:1176690/1‑199 (MQ=11)
                                                                                                                                         tgcggttatcactgaggTGATTGTGGTGGGGGAAGAGCTACCTAGCGCTGCGCGCTTCTATCTGCGTGGTGTTGCGTGAGGCAACGGTCTCTCTCTGCCGCTCGAGCCGAAAGTTAAGCATAGCGTGTCATCGTTACCCGCATGGGCATAATATGTGATTGCGGATGCGGATGAGAGTTCAGTGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGcc                                                                                                  <  2:147095/223‑1 (MQ=9)
                                                                                                                                         tCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGc                                                                                                   >  1:133776/1‑238 (MQ=255)
                                                                                                                                                  tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTccc                                                                                                                                                                                                       >  1:588228/1‑129 (MQ=17)
                                                                                                                                                   ttGAAGGGGATGGTGGGGGGGGAAGGATTACTCAGCCCTCGCCGCTCGGCCGTTGCGGCCGTTGCCTGCAGGAGCGCTCTCTCGCTGGCGCTCGAGTCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCg                                                                                        <  2:587823/239‑1 (MQ=11)
                                                                                                                                                   ttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCg                                                                                        <  1:1234501/239‑1 (MQ=255)
                                                                                                                                                        aGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGa                                                                                                                                                                                           >  1:207637/1‑135 (MQ=37)
                                                                                                                                                                           aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGTCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                                                                                                                             >  2:1216862/1‑115 (MQ=17)
                                                                                                                                                                           aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGAATATTTGATTGCGGATTCGCTGGGGAGTTCAGGGACTGTTGATAATGATCGTGGTGGGCTGAGGGTTACTCGGGGGTTCGCGCTTCGCCCTTCGGTCCCTTTCCTGGGCGAACGctctcgc                                                                  >  2:227990/1‑236 (MQ=14)
                                                                                                                                                                           aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGTTCGTTGCCTGCGGCAACGCcctctc                                                                  >  1:781658/1‑238 (MQ=255)
                                                                                                                                                                           aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGACTGCGGCAACGCTCTCTCGCTGACGCTCGCACCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGAATGGTCCCACTATTTACTTGAGGACTCGATTAACCGATCAGGGACTTTGTAAAGTGGTGGTGGTGGGAGCAGGATTACCTAGCGCTGCGACGATCCCCCATCGGCACCTTGCCTGCgggcaccacccccgg                                                                 >  2:1121297/1‑226 (MQ=9)
                                                                                                                                                                                   tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                         <  1:1336240/239‑1 (MQ=255)
                                                                                                                                                                                       cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAg                                                     >  1:786877/1‑239 (MQ=255)
                                                                                                                                                                                       cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAg                                                     >  1:23915/1‑239 (MQ=255)
                                                                                                                                                                                       cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCCCTTGGGCCGACTATTTGATTGGGGGATCGCTTAAGTGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCCTCGCCCTTTGGGTCGTTGCGTGCGGCACCGCTCTCTCGCCgggccccggg                                                     >  2:861456/1‑231 (MQ=14)
                                                                                                                                                                                        gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                                                                                          >  2:958987/1‑169 (MQ=255)
                                                                                                                                                                                        gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                                                                                          <  1:958987/169‑1 (MQ=255)
                                                                                                                                                                                        gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAg                                                     >  1:733456/1‑238 (MQ=255)
                                                                                                                                                                                            cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTAAGAGCTGCGCGCATCGCGCTTCAGGTAGTGGCCGGAGGCCACGGTCTCTCGCGGgcgccccagccgg                                                 >  2:1271486/1‑229 (MQ=14)
                                                                                                                                                                                            cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGa                                                 >  1:394284/1‑238 (MQ=255)
                                                                                                                                                                                                cTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCCTCAGCGCTGCGCGCTTCGCACTTCGGTTCGTTGCCTGCGGCATCGGTCTCTCGCGGGCGCTCgcgcgggacat                                             >  2:1015567/1‑228 (MQ=11)
                                                                                                                                                                                                       cTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGaa                                     <  1:861456/239‑1 (MQ=255)
                                                                                                                                                                                                       cTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGCCGaa                                     <  1:858186/239‑1 (MQ=255)
                                                                                                                                                                                                                     cTGCGGCAACGCTCTCTCGCTGTCGCTCTAGTCGAACCTTAGTCGAAGCTTCTCATCCTTCTCCGTTTGGGCGGAATGTTTGATGGCGGATTCGGTTGAGATTTTAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTc                       <  2:166631/239‑1 (MQ=11)
                                                                                                                                                                                                                       gCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTcc                      >  1:339078/1‑238 (MQ=255)
                                                                                                                                                                                                                               gccttcccGCTGACGCTCGAGCCGAACCTTGGTCGAAGCGTCTCACCCTTCCCCGCAGGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGTGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATggg             <  1:16124/232‑1 (MQ=255)
                                                                                                                                                                                                                               gCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATggg             <  2:23915/239‑1 (MQ=255)
                                                                                                                                                                                                                                 tgtcgcgcgGACGCTCGAGCCGATCTTTAGTCGAAGGTTCTCATCCTTCTGCGCATGGGAAGGATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa           <  2:408097/234‑1 (MQ=14)
                                                                                                                                                                                                                                 tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa           <  2:133776/239‑1 (MQ=255)
                                                                                                                                                                                                                                 tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa           <  2:339078/239‑1 (MQ=255)
                                                                                                                                                                                                                                   tctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGa         <  1:236383/239‑1 (MQ=255)
                                                                                                                                                                                                                                          gACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGCGAGATCAGGGACTTTTGAAAGTGATGGGGGTGGGGGAAGGATTACTCAGCGCTGAGCGGTTCGCCCTTGGGTTAGGTGGCTGCGGCCACGCTCTCTCGTCGGCGCTCGAATGGGACCTTGACAGAGTCATGTCATCCTTcctctgctgggaagggacactc  >  2:146979/1‑219 (MQ=18)
                                                                                                                                                                                                                                          gACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTc                                                                                                                                                             >  2:1043585/1‑84 (MQ=255)
                                                                                                                                                                                                                                          gACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTc                                                                                                                                                             <  1:1043585/84‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              
TCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTA  >  NC_000913/1286731‑1287202

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