Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,622,348:1 +A 6.8% coding (115/639 nt) purN → phosphoribosylglycinamide formyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,622,3481.A 6.8% 786.4 / 55.1 191coding (115/639 nt)purNphosphoribosylglycinamide formyltransferase 1
Reads supporting (aligned to +/- strand):  ref base . (91/87);  new base A (6/7);  total (97/94)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.81e-01

Only 100 of 195 total aligned reads displayed.

CTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTT  >  NC_000913/2622201‑2622496
                                                                                                                                                    |                                                                                                                                                    
ctctGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑a                                                                                                                                                     >  3:1071397/1‑149 (MQ=255)
ctctGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑a                                                                                                                                                     >  4:277693/1‑149 (MQ=255)
 tctGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑Ac                                                                                                                                                    >  4:1001508/1‑149 (MQ=255)
 tctGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑Ac                                                                                                                                                    >  2:800469/1‑149 (MQ=255)
      ttCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGtt                                                                                                                                               >  1:1202395/1‑149 (MQ=255)
       tCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTc                                                                                                                                              <  4:1482493/149‑1 (MQ=255)
       tCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTc                                                                                                                                              <  2:961083/149‑1 (MQ=255)
         cGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAATCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCgg                                                                                                                                            >  3:1226953/1‑149 (MQ=255)
         cGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCgg                                                                                                                                            <  1:468972/149‑1 (MQ=255)
         cGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCgg                                                                                                                                            <  4:229703/149‑1 (MQ=255)
          gAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGc                                                                                                                                           >  4:1107569/1‑149 (MQ=255)
           aacaacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGcc                                                                                                                                          <  4:651871/149‑1 (MQ=255)
             catcgcgTGGTTATCGAATAATGAATATTGGGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGTCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATATGGCCG‑ACGCGTTCGGCCtt                                                                                                                                        <  2:751586/146‑1 (MQ=255)
             caacGTGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCG‑TCGGCCtt                                                                                                                                        <  3:371796/149‑1 (MQ=255)
             caacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCtt                                                                                                                                        <  3:1107569/149‑1 (MQ=255)
               acGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCG‑TCGGCCttat                                                                                                                                      <  4:582957/149‑3 (MQ=255)
                cgcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGaa                                                                                                                                     <  2:1104708/149‑1 (MQ=255)
                cgcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGaa                                                                                                                                     <  2:1233726/149‑1 (MQ=255)
                 gcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTg                                                                                                                                       >  3:653542/1‑146 (MQ=255)
                     ggTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGACGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACgcgc                                                                                                                                <  3:1453024/149‑1 (MQ=255)
                        tATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGcgcccg                                                                                                                             >  2:1349763/1‑149 (MQ=255)
                         aTCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGcgcccgc                                                                                                                            >  2:1497379/1‑149 (MQ=255)
                              ataatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGc                                                                                                                       >  3:222184/1‑149 (MQ=255)
                                aatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCgg                                                                                                                     <  2:1046181/149‑1 (MQ=255)
                                 atGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCggg                                                                                                                    <  2:153901/149‑1 (MQ=255)
                                   gAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACTCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTa                                                                                                                  <  3:877654/149‑1 (MQ=255)
                                   gAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTa                                                                                                                  >  1:1363177/1‑149 (MQ=255)
                                     atatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTAtt                                                                                                                <  4:1226953/149‑1 (MQ=255)
                                      tatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCG‑TCGGCCttat                                                                                                                                      >  2:309204/1‑124 (MQ=255)
                                      tatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCgg                                                                                                                     <  2:1641138/143‑1 (MQ=255)
                                      tatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTg                                                                                                               >  2:856605/1‑149 (MQ=255)
                                       atTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGc                                                                                                              >  2:1528444/1‑149 (MQ=255)
                                          gtggtgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAAc                                                                                                           <  4:1021555/149‑1 (MQ=255)
                                           tggtgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACg                                                                                                          >  2:579006/1‑149 (MQ=255)
                                              tgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCAt                                                                                                       >  2:945595/1‑149 (MQ=255)
                                              tgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCAt                                                                                                       >  4:704770/1‑149 (MQ=255)
                                                 ttATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCt                                                                                                                                         >  3:454175/1‑112 (MQ=255)
                                                 ttATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCt                                                                                                                                         <  4:454175/112‑1 (MQ=255)
                                                     ttCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAACTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATa                                                                                                      >  4:1255242/1‑143 (MQ=255)
                                                     ttCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAAACAACAAAACTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATa                                                                                                      <  3:1255242/143‑1 (MQ=255)
                                                          gCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGcc                                                                                           >  2:369709/1‑149 (MQ=255)
                                                           cAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCa                                                                                          <  4:381116/149‑1 (MQ=255)
                                                            aaCGGAAGTAATTTACAGGCAATTATTGACGCCTGT‑AAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAgc                                                                                        >  2:996770/1‑149 (MQ=255)
                                                                 aaGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGt                                                                                    <  1:1528444/149‑1 (MQ=255)
                                                                      aTTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGAc                                                                               >  1:784633/1‑149 (MQ=255)
                                                                       tttACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACa                                                                              <  2:1363177/149‑1 (MQ=255)
                                                                       tttACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACa                                                                              <  1:1622087/149‑1 (MQ=255)
                                                                              tcAATTATTGCCGACTTTAACACCAACAAAATTAAAGGCACTGTACGAGTAGTTTTCAGCAATAAGTCCG‑ACCCTTTCGGCCTTGAACGCGCCCGCCAGGCGGGACTTGCAACACTTACGCTCATCGCCAGCGCGTGTGCCCGTCgagt                                                                       >  4:108943/2‑146 (MQ=255)
                                                                              gCAATTATTGACGCCTGTAAACCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCG‑TCGGCCttatt                                                                                                                                     <  1:58568/87‑4 (MQ=38)
                                                                              gCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGa                                                                                <  4:958947/140‑1 (MQ=255)
                                                                              gCAATTATTGAAGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGa                                                                                >  3:958947/1‑140 (MQ=255)
                                                                               cAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGaa                                                                      >  2:1613086/1‑149 (MQ=255)
                                                                                        aCGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGAc                                                             <  3:1572894/149‑1 (MQ=255)
                                                                                          gCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCg                                                           >  3:1646351/1‑149 (MQ=255)
                                                                                           ccTGTAAAACCAACAAGATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgg                                                          <  4:552746/149‑1 (MQ=255)
                                                                                            cTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCtt                                                                                                                                        <  2:494494/70‑1 (MQ=255)
                                                                                             tGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGa                                                        <  3:421587/149‑1 (MQ=255)
                                                                                             tGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGa                                                        <  2:833344/149‑1 (MQ=255)
                                                                                               tAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTc                                                                           >  1:96564/1‑128 (MQ=255)
                                                                                               tAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGt                                                      >  3:546865/1‑149 (MQ=255)
                                                                                                  aaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAGCGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGa                                                   <  3:270751/149‑1 (MQ=255)
                                                                                                  aaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGa                                                   <  3:1002829/149‑1 (MQ=255)
                                                                                                   aCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGGTTGACAGTCGTGAAGCCTATGACAGGGAGTTGAt                                                  >  1:341714/1‑149 (MQ=255)
                                                                                                    ccAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGAtt                                                 <  3:303522/149‑1 (MQ=255)
                                                                                                      aacaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGc                                                                    <  1:813101/128‑1 (MQ=255)
                                                                                                      aacaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCa                                               >  2:6082/1‑149 (MQ=255)
                                                                                                      aacaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCa                                               >  1:855492/1‑149 (MQ=255)
                                                                                                      aacaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCa                                               >  1:1031823/1‑149 (MQ=255)
                                                                                                       acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATAc                                                                                                     <  1:1301407/94‑1 (MQ=255)
                                                                                                       acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGt                                                      <  1:869250/141‑1 (MQ=255)
                                                                                                       acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCAt                                              <  1:6082/149‑1 (MQ=255)
                                                                                                       acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCAt                                              >  4:1325388/1‑149 (MQ=255)
                                                                                                         aaaaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGa                                            >  1:621488/1‑149 (MQ=255)
                                                                                                         aaaaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGa                                            >  3:1631333/1‑149 (MQ=255)
                                                                                                          aaaTTAAAGGCACCGTACGGGCAGTTTTCATCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                         <  4:1160664/119‑1 (MQ=255)
                                                                                                          aaaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                         >  3:1160664/1‑119 (MQ=255)
                                                                                                             ttAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATc                                        >  4:948366/1‑149 (MQ=255)
                                                                                                             ttAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATc                                        >  3:719542/1‑149 (MQ=255)
                                                                                                                aaGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGAc                                     <  1:297827/149‑1 (MQ=255)
                                                                                                                  ggCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAg                                                       <  2:433599/129‑1 (MQ=255)
                                                                                                                  ggCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAg                                                       >  1:433599/1‑129 (MQ=255)
                                                                                                                  ggCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGtt                                                     <  3:1432830/131‑1 (MQ=255)
                                                                                                                  ggCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACAt                                   >  2:664237/1‑149 (MQ=255)
                                                                                                                     aCCGTACGGGCAGTTTTCAGCAATAAGGTCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTa                                >  4:1471180/1‑149 (MQ=255)
                                                                                                                     aCCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGa                                      <  1:191199/143‑1 (MQ=255)
                                                                                                                         tACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCa                            >  4:1585324/1‑149 (MQ=255)
                                                                                                                              gCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGa                       >  4:1629710/1‑149 (MQ=255)
                                                                                                                                       aGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGAtt                                                 <  4:870714/114‑1 (MQ=255)
                                                                                                                                          aaTAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgct           >  1:1011076/1‑149 (MQ=255)
                                                                                                                                          aaTAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgct           >  4:409931/1‑149 (MQ=255)
                                                                                                                                          aaTAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgct           >  2:247558/1‑149 (MQ=255)
                                                                                                                                           aTAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctg          >  2:762798/1‑149 (MQ=255)
                                                                                                                                           aTAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctg          >  3:808517/1‑149 (MQ=255)
                                                                                                                                            tAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctgg         <  1:996770/149‑1 (MQ=255)
                                                                                                                                            caaGGCCGAACGCG‑TCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctgg         >  4:432169/2‑149 (MQ=255)
                                                                                                                                             aaGGCCGAACGCG‑TCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctggc        >  3:858715/1‑149 (MQ=255)
                                                                                                                                              aGGCCGAACGCG‑TCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctggct       >  4:31090/1‑149 (MQ=255)
                                                                                                                                                 ccG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGt    >  4:1367685/1‑149 (MQ=255)
                                                                                                                                                  cG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGtt   >  1:901946/1‑149 (MQ=255)
                                                                                                                                                   g‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGttt  >  1:1051998/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTT  >  NC_000913/2622201‑2622496

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑