Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 342,144 T→C 8.7% I340I (ATT→ATC yahJ → putative deaminase with metallo‑dependent hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913342,1440TC 8.7% 195.5 / 8.2 69I340I (ATT→ATCyahJputative deaminase with metallo‑dependent hydrolase domain
Reads supporting (aligned to +/- strand):  ref base T (29/34);  new base C (2/4);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

AAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGCCAATCTGTACGCGCAGCTCTATATTCG  >  NC_000913/341999‑342290
                                                                                                                                                 |                                                                                                                                                  
aaaCGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCa                                                                                                                                                >  2:360847/1‑150 (MQ=255)
  aCGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCAATCGGCAcc                                                                                                                                              <  1:1003256/150‑2 (MQ=255)
     gTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTg                                                                                                                                           <  1:1359236/150‑1 (MQ=255)
        gagaAACCGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCAATCGGCAccgtagag                                                                                                                                        <  2:754915/150‑8 (MQ=255)
        gagaAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCat                                                                                                                                        <  1:666117/150‑1 (MQ=255)
        gagaAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCAATCGGCAccgtcgag                                                                                                                                        <  1:1302573/150‑8 (MQ=255)
           aaaaCGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATg                                                                                                                                     <  2:645242/150‑1 (MQ=255)
                gCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCt                                                                                                                                >  1:676171/1‑150 (MQ=255)
                  cacaGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCa                                                                                                                              <  2:320505/150‑1 (MQ=255)
                        tGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACTGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGt                                                                                                                        <  1:1047157/150‑1 (MQ=255)
                        tGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGt                                                                                                                        >  1:1063004/1‑150 (MQ=255)
                         gAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGtt                                                                                                                       >  1:1247763/1‑150 (MQ=255)
                           aGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGc                                                                                                                     <  1:886424/150‑1 (MQ=255)
                           aGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGc                                                                                                                     <  1:360847/150‑1 (MQ=255)
                            gggCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCa                                                                                                                    >  2:276556/1‑150 (MQ=255)
                             ggCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCAc                                                                                                                   <  1:49611/150‑1 (MQ=255)
                             ggCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCAc                                                                                                                   <  1:1060697/150‑1 (MQ=255)
                              gCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACg                                                                                                                  <  2:1149769/150‑1 (MQ=255)
                              gCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACg                                                                                                                  <  1:936996/150‑1 (MQ=255)
                                  gCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACaa                                                                                                              >  1:91986/1‑150 (MQ=255)
                                     gACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAgg                                                                                                           <  1:759875/150‑1 (MQ=255)
                                       ccATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCg                                                                                                         >  1:884718/1‑150 (MQ=255)
                                        cATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGt                                                                                                        >  1:659294/1‑150 (MQ=255)
                                              tCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGt                                                                                                  <  2:793834/150‑1 (MQ=255)
                                              tCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGt                                                                                                  <  2:788826/150‑1 (MQ=255)
                                               cACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGtg                                                                                                 >  1:273259/1‑150 (MQ=255)
                                                  gCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGtgatg                                                                                              <  2:95196/150‑1 (MQ=255)
                                                     tttGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACt                                                                                           <  2:1063004/150‑1 (MQ=255)
                                                     tttGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACt                                                                                           >  2:1133736/1‑150 (MQ=255)
                                                                  cGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCg                                                                              <  1:1380788/150‑1 (MQ=255)
                                                                  cGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCg                                                                              >  2:1401815/1‑150 (MQ=255)
                                                                  cGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCg                                                                              >  2:1327367/1‑150 (MQ=255)
                                                                        aCGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCGCAAACCGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCg                                                                        <  1:726893/150‑1 (MQ=255)
                                                                           aGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGc                                                                                                                                    <  1:319933/87‑1 (MQ=255)
                                                                           aGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGc                                                                                                                                    >  2:319933/1‑87 (MQ=255)
                                                                             cAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCAc                                                                   <  1:276556/150‑1 (MQ=255)
                                                                              aaCAGGTAGATGAACTGGCGAACCGGATGGAGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACt                                                                  >  2:697120/1‑150 (MQ=255)
                                                                                 aGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGt                                                               <  2:659294/150‑1 (MQ=255)
                                                                                  ggTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTc                                                              >  2:1420126/1‑150 (MQ=255)
                                                                                            cTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGt                                                    >  2:993053/1‑150 (MQ=255)
                                                                                            cTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGt                                                    >  2:507997/1‑150 (MQ=255)
                                                                                                 gAACCGGATGTTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTggg                                               <  2:426187/150‑1 (MQ=255)
                                                                                                 gAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTggg                                               >  2:182119/1‑150 (MQ=255)
                                                                                                     cGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGGCTGGGCGAc                                           <  2:749075/150‑1 (MQ=255)
                                                                                                      ggATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTa                                                                           >  2:152949/1‑117 (MQ=255)
                                                                                                      ggATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTa                                                                           <  1:152949/117‑1 (MQ=255)
                                                                                                      ggATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACa                                          <  1:1133736/150‑1 (MQ=255)
                                                                                                             ggtgCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGACGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGa                                   >  2:513682/1‑150 (MQ=255)
                                                                                                               tgCAACAAATTTCTATCTCCTCGACGGTGCCGATTGGCACGCTGCATATTCCGCTCAAACATTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGaaa                                 <  2:977633/150‑1 (MQ=255)
                                                                                                                gCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGaaaa                                >  1:1354790/1‑150 (MQ=255)
                                                                                                                  aacaaATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAg                              >  2:344629/1‑150 (MQ=255)
                                                                                                                   acaaATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGc                             <  1:697120/150‑1 (MQ=255)
                                                                                                                   acaaATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGc                             <  2:1247763/150‑1 (MQ=255)
                                                                                                                   acaaATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGc                             <  1:325411/150‑1 (MQ=255)
                                                                                                                   acaaATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGc                             >  1:763340/1‑150 (MQ=255)
                                                                                                                    caaATTTCTATCGCCTCGACGGTTCCGATTGGCACGCTGCATATTCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGcc                            <  2:1048320/150‑1 (MQ=255)
                                                                                                                      aaTTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGCCaa                          >  1:609083/1‑150 (MQ=255)
                                                                                                                         ttCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGCCAATCt                       >  1:928864/1‑150 (MQ=255)
                                                                                                                           cTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGCCAATCTGt                     <  1:890681/150‑1 (MQ=255)
                                                                                                                                 ccTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGCCAATCTGTACgcgc               <  2:884718/150‑1 (MQ=255)
                                                                                                                                 ccTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGCCAATCTGTACgcgc               >  2:958629/1‑150 (MQ=255)
                                                                                                                                   tgcagcgTGCCAATCGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGCCAATCTGTACGCGCAg             >  1:125467/7‑150 (MQ=255)
                                                                                                                                      agcgTGCCAATCGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGCCAATCTGTACGCGCAGctc          <  2:125467/147‑1 (MQ=255)
                                                                                                                                       cGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGCCAATCTGTACGCGCAGctct         <  2:91986/150‑1 (MQ=255)
                                                                                                                                         gTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGCCAATCTGTACGCGCAGCTCtat       >  2:1305978/1‑150 (MQ=255)
                                                                                                                                          tGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGCCAATCTGTACGCGCAGCTCtata      >  2:946990/1‑150 (MQ=255)
                                                                                                                                           gCCAATCGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGCCAATCTGTACGCGCAGCTCTatat     >  2:703673/1‑150 (MQ=255)
                                                                                                                                             cGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGCCAATCTGTACGCGCAGCTCTATATTc   >  1:387743/1‑150 (MQ=255)
                                                                                                                                              gATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGCCAATCTGTACGCGCAGCTCTATATTCg  <  2:370895/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGCATATGCCGCTCAAACAGTTGCACGACAAAGGCGTAAAAGTGATGACTGGCACTGACAGCGTTATCGACCACTGGTCGCCTTATGGTCTGGGCGACATGCTGGAAAAAGCCAATCTGTACGCGCAGCTCTATATTCG  >  NC_000913/341999‑342290

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