Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,314,779:1 +T 10.9% coding (937/1209 nt) aspC → aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,314,7791.T10.9% 325.2 / 37.0 92coding (937/1209 nt)aspCaminotransferase
Reads supporting (aligned to +/- strand):  ref base . (35/47);  new base T (6/4);  total (41/51)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

CAAGCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTG  >  NC_002947/4314631‑4314928
                                                                                                                                                     |                                                                                                                                                     
cAAGCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATgg                                                                                                                                                      >  1:553066/1‑150 (MQ=255)
   gCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCt                                                                                                                                                   <  1:1273539/150‑1 (MQ=255)
    cTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTg                                                                                                                                                  <  1:1636926/150‑1 (MQ=255)
     tCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTgg                                                                                                                                                 >  2:148583/1‑150 (MQ=255)
      cAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGa                                                                                                                                                <  2:982725/150‑1 (MQ=255)
      cAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGa                                                                                                                                                <  2:1308679/150‑1 (MQ=255)
      cAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGa                                                                                                                                                <  2:293215/150‑1 (MQ=255)
        gCCTGCTGCATGAAAGTGGGCAGCCGCGTATTCTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAgg                                                                                                                                              <  1:744472/150‑1 (MQ=255)
         ccTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGgt                                                                                                                                             <  1:936264/150‑1 (MQ=255)
          ctgctgCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGgtg                                                                                                                                            <  1:1487013/150‑1 (MQ=255)
                 aTGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACg                                                                                                                                     >  2:2341522/1‑150 (MQ=255)
                 aTGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACg                                                                                                                                     >  2:1288305/1‑150 (MQ=255)
                  tGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGc                                                                                                                                    <  2:1698961/150‑1 (MQ=255)
                  tGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGc                                                                                                                                    >  2:1374750/1‑150 (MQ=255)
                        ggggCGGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATGCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACgg                                                                                                                              <  1:182919/149‑1 (MQ=255)
                               ccGCTTATTTTACTCGATGCCGTGACCTGCATCCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTc                                                                                                                       <  1:1298636/150‑1 (MQ=255)
                               ccGCGTATTTTACTCGATGCCGTGACCTGCATCCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTc                                                                                                                       <  1:711086/150‑1 (MQ=255)
                                 gcgTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCtt                                                                                                                     >  1:1476667/1‑150 (MQ=255)
                                  cgTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCttt                                                                                                                    <  1:148583/150‑1 (MQ=255)
                                  cgTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCctt                                                                                                                    <  1:696445/150‑1 (MQ=255)
                                      ttttACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGttt                                                                                                                >  2:754741/1‑150 (MQ=255)
                                       tttACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCCATGTGCCTGGa                                                                                                                                                <  1:2156324/117‑1 (MQ=255)
                                       tttACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCCATGTGCCTGGa                                                                                                                                                >  2:2156324/1‑117 (MQ=255)
                                          aCTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGt                                                                                                            <  2:553066/150‑1 (MQ=255)
                                            tCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGc                                                                                                          >  2:2232960/1‑150 (MQ=255)
                                             ggtttccGTGTGCTGCATTCTCGACGACAGAGGCTTAGACGCACGCCTTTGTCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCg                                                                                                         <  2:396820/145‑1 (MQ=255)
                                               aTGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGAt                                                                                                       >  1:813164/1‑150 (MQ=255)
                                                        ccTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCgg                                                                                              <  1:1339803/150‑1 (MQ=255)
                                                            cATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCtgt                                                                                          <  2:2149826/150‑1 (MQ=255)
                                                             aTTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGt                                                                                                            >  1:1614299/1‑131 (MQ=255)
                                                             aTTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGt                                                                                                            <  2:1614299/131‑1 (MQ=255)
                                                              ttCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGa                                                                                        <  2:1272754/150‑1 (MQ=255)
                                                              ttCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAAAGGTGGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATACTTTGTTTACGTGCGATGTCCTGCGGcgggga                                                                                        >  1:1780278/1‑145 (MQ=255)
                                                                 tcGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGgt                                                                                                                                             <  1:1659704/94‑1 (MQ=255)
                                                                 tcGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGgt                                                                                                                                             >  2:1659704/1‑94 (MQ=255)
                                                                 tcGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCa                                                                                     >  2:1366005/1‑150 (MQ=255)
                                                                      gacCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGc                                                                                >  2:1282169/1‑150 (MQ=255)
                                                                      gacCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGATCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGt                                                                                                            >  2:356559/1‑122 (MQ=255)
                                                                      gacCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGATCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGt                                                                                                            <  1:356559/122‑1 (MQ=255)
                                                                       acCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCa                                                                               <  1:754741/150‑1 (MQ=255)
                                                                        ccGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAg                                                                              <  1:2383186/150‑1 (MQ=255)
                                                                         cGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACgg                                                                                                                              <  1:1484709/101‑1 (MQ=255)
                                                                         cGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACgg                                                                                                                              >  2:1484709/1‑101 (MQ=255)
                                                                         cGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAgg                                                                             >  2:107126/1‑150 (MQ=255)
                                                                             gCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACAt                                                                         <  2:844879/150‑1 (MQ=255)
                                                                              cTTCGCCGAACGCCTTTGCCAGTGCATATCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACAtt                                                                        >  1:2132582/1‑150 (MQ=255)
                                                                              cTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGa                                                                            >  1:445353/1‑146 (MQ=255)
                                                                              cTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGa                                                                            <  2:445353/146‑1 (MQ=255)
                                                                               ttCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGGCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTc                                                                       <  1:1366005/150‑1 (MQ=255)
                                                                                 cGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGTAGTATTGTAACGTTGTCTCCACGCACATGGCCTGGAGGGGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGc                                                                     <  1:2341522/150‑1 (MQ=255)
                                                                                   ccGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTCCACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCa                                                                   <  1:1374750/150‑1 (MQ=255)
                                                                                       aCGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGc                                                               <  1:1288305/150‑1 (MQ=255)
                                                                                            tttGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTaaa                                                                                                                                       <  1:1282189/73‑1 (MQ=255)
                                                                                            tttGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTaaa                                                                                                                                       >  2:1282189/1‑73 (MQ=255)
                                                                                             ttGCCAGTGCATAGCTCAGCGACGTGCGGTACTGGCCGAAATTCCTCAGGATCATGCCCTGGAGCTTTTAAACGCTCACGCCACATCCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGa                                                         <  2:1476667/150‑1 (MQ=255)
                                                                                             ttGCCAGGGCATAACTCAGCGACGTGCTGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGa                                                         <  1:799471/150‑1 (MQ=255)
                                                                                                ccAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGAACTTGCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACtt                                                      >  1:1956957/1‑150 (MQ=255)
                                                                                                 cAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACttt                                                     >  1:1567300/1‑150 (MQ=255)
                                                                                                         aGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTaa                                                                                                                                        <  2:1980673/59‑1 (MQ=255)
                                                                                                         aGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTaa                                                                                                                                        >  1:1980673/1‑59 (MQ=255)
                                                                                                         aGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATa                                             <  2:47357/150‑1 (MQ=255)
                                                                                                            tCAGCGACGTGCGGTAATAGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCt                                          >  2:1141535/1‑150 (MQ=255)
                                                                                                                 gACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCa                                     <  2:813164/150‑1 (MQ=255)
                                                                                                                      gCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAg                                >  1:185375/1‑150 (MQ=255)
                                                                                                                       cGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGc                                                 <  2:2323392/132‑1 (MQ=255)
                                                                                                                       cGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGc                                                 >  1:2323392/1‑132 (MQ=255)
                                                                                                                         gTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTg                             <  1:107126/150‑1 (MQ=255)
                                                                                                                          tAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGGCATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTgg                            <  1:1282169/150‑1 (MQ=255)
                                                                                                                               gTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTg                             <  1:1732063/144‑1 (MQ=255)
                                                                                                                               gTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATc                       >  1:227092/1‑150 (MQ=255)
                                                                                                                               gTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTACGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTg                             >  2:1732063/1‑144 (MQ=255)
                                                                                                                                tCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATcc                      >  1:313121/1‑150 (MQ=255)
                                                                                                                                 cGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCa                                                             <  2:1918341/110‑1 (MQ=255)
                                                                                                                                 cGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCa                                                             >  1:1918341/1‑110 (MQ=255)
                                                                                                                                  gATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTg                    >  1:597525/1‑150 (MQ=255)
                                                                                                                                   atatTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGa                   <  2:1956957/150‑1 (MQ=255)
                                                                                                                                     cacctcCAGGCCATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTc                                                                                                    <  1:868510/64‑1 (MQ=37)
                                                                                                                                     cacctcCAGGCCATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTc                                                                                                    >  2:868510/4‑67 (MQ=37)
                                                                                                                                      ttCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGa                <  1:2232960/150‑1 (MQ=255)
                                                                                                                                      acctcCAGGCCATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGa                                                                            <  1:1245704/88‑1 (MQ=255)
                                                                                                                                      acctcCAGGCCATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGa                                                                            >  2:1245704/3‑90 (MQ=255)
                                                                                                                                        ctcCAGGCCATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGAtt              >  2:1737660/1‑150 (MQ=255)
                                                                                                                                        ctcCAGGCCATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGAtt              >  2:811338/1‑150 (MQ=255)
                                                                                                                                         tcCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCtgt                                                                                          >  1:454607/1‑73 (MQ=255)
                                                                                                                                         tcCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCtgt                                                                                          <  2:454607/73‑1 (MQ=255)
                                                                                                                                         tcCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCa                                                                               <  1:180904/84‑1 (MQ=255)
                                                                                                                                         tcCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCa                                                                               >  2:180904/1‑84 (MQ=255)
                                                                                                                                         tcCAGGCCATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGa                                                                                        <  1:708965/75‑1 (MQ=255)
                                                                                                                                         tcCAGGCCATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGa                                                                                        >  2:708965/1‑75 (MQ=255)
                                                                                                                                           cAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTccc           >  2:120230/1‑150 (MQ=255)
                                                                                                                                              gCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTAt        <  2:227092/150‑1 (MQ=255)
                                                                                                                                                acaTGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCa      >  2:1692402/1‑150 (MQ=255)
                                                                                                                                                  aTGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAgg    >  1:2181828/1‑150 (MQ=255)
                                                                                                                                                    ggCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTg  >  2:2362517/1‑150 (MQ=255)
                                                                                                                                                    ggCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTg  >  1:1680281/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAAGCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATGGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTG  >  NC_002947/4314631‑4314928

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