Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 398,532 C→T P221S (CCC→TCC)  yaiW → surface‑exposed outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913398,5320CT97.6% 130.2 / NA 42P221S (CCC→TCC) yaiWsurface‑exposed outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (20/20);  minor base G (0/1);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTCGGTGCGTACAGAAAAGCAATTGAGCGCGATTTTTGACGATCTGATCAACATGGTGCCAATGGGGCAGACGCTGTTTGGTTCGCTCAATCCGGTGCGCACCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATCCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGCTGGCACTGGACGGTGATTTAATTCGTTACGACAGTAAAGAACCCGGGAAAACCGAACTGGCAACGCGCAAACTGGCTGCAAAACTGGGGATGAGCGACAGCGAAATTCGCCGTCAGTTAGAGAAAGG  >  NC_000913/398286‑398770
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
gATTCGGTGCGTACAGAAAAGCAATTGAGCGCGATTTTTGACGATCTGATCAACATGGTGCCAATGGGGCAGACGCTGTTTGGTTCGCTCAATCCGGTGCGCACCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGc                                                                                                                                                                                                                                            >  1:348848/1‑251 (MQ=255)
      gTGCGTACAGAAAAGCAATTGAGCGCGATTTTTGACGATCTGATCAACATGGTGCCAATGGGGCAGACGCTGTTTGGTTCGCTCAATCCGGTGCGCACCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTta                                                                                                                                                                                                                                      <  1:389597/251‑1 (MQ=255)
       tGCGTACAGAAAAGCAATTGAGCGCGATTTTTGACGATCTGATCAACATGGTGCCAATGGGGCAGACGCTGTTTGGTTCGCTCAATCCGGTGCGCACCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCCGTtat                                                                                                                                                                                                                                     >  2:203694/1‑251 (MQ=255)
          gTACAGAAAAGCAATTGAGCGCGATTTTTGACGATCTGATCAACATGGTGCCAATGGGGCAGACGCTGTTTGGTTCGCTCAATCCGGTGCGCACCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAgc                                                                                                                                                                                                                                  <  1:213913/251‑1 (MQ=255)
                       aTTGAGCGCGATTTTTGACGATCTGATCAACATGGTGCCAATGGGGCAGACGCTGTTTGGTTCGCTCAATCCGGTGCGCACCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATAcc                                                                                                                                                                                                                     <  1:203694/251‑1 (MQ=255)
                       aTTGAGCGCGATTTTTGACGATCTGATCAACATGGTGCCAATGGGGCAGACGCTGTTTGGTTCGCTCAATCCGGTGCGCACCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATAcc                                                                                                                                                                                                                     <  1:324487/251‑1 (MQ=255)
                          gAGCGCGATTTTTGACGATCTGATCAACATGGTGCCAATGGGGCAGACTCTGTTTGGTTCGCTCAAGCCGGTGCGCACCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTTGAAAATGGACGGTACGGTCCGTCAGGAAGTCTTCAGTCGCCGCTGCGGGTTGTGGTTTGGTACTTACCATTTACTGAGTTATTCCGCCAGTTATAGCGCACCGATATACCGt                                                                                                                                                                                                                   <  1:32626/250‑1 (MQ=255)
                          gAGCGCGATTTTTGACGATCTGATCAACATGGTGCCAATGGGGCAGACGCTGTTTGGTTCGCTCAATCCGGTGCGCACCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGt                                                                                                                                                                                                                   <  2:219653/250‑1 (MQ=255)
                                 aTTTTTGACGATCTGCTCAACATGGTGCCAATGGGGCAGACGCTGTTTGGTTCGCTCAATCCGGTGCGCACCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGCCGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGAGGGGTGGGGGTTGGGACTTACCCTTTTATGAATTATTCCGCCAGGTATAGCGCACCGACATACCGTTTTGCTg                                                                                                                                                                                                            >  2:325574/1‑250 (MQ=255)
                                      tGACGATCTGATCAACATGGTGCCAATGGGGCAGACGCTGTTTGGTTCGCTCAATCCGGTGCGCACCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGCAGTCTTCAGTCGCCGCGGCGGGGTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTTTAGCGCACCGat                                                                                                                                                                                                                          >  2:221307/1‑231 (MQ=255)
                                      tGACGATCTGATCAACATGGTGCCAAGGGGGCAGACGCTGTTTGGTTCGCTCAATCCGGTGCGCACCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGat                                                                                                                                                                                                                          <  1:221307/231‑1 (MQ=255)
                                                           gCCAATGGGGCAGACGCTGTTTGGTTCGCTCAATCCGGTGCGCCCCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTCCAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCAACCATATACCGTTTTGCTGATTTTAACGCTGGCTGGTAcgcccgcc                                                                                                                                                                                 >  1:370612/1‑251 (MQ=255)
                                                           gCCAATGGGGCAGACGCTGTTTGGTTCGCTCAATCCGGTGCGCACCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGcc                                                                                                                                                                                 <  1:279239/251‑1 (MQ=255)
                                                            ccAATGGGGCAGACGCTGTTTGGTTCGCTCAATCCGGTGCGCACCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCCGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCg                                                                                                                                                                                >  1:281938/1‑251 (MQ=255)
                                                                                         cAATCCGGTGCGCACCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGt                                                                                                                                                    >  1:267201/1‑250 (MQ=255)
                                                                                           aTCCGGTGCGCACCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTa                                                                                                                                                   >  1:165801/1‑249 (MQ=255)
                                                                                              cGGTGCGCACCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGcc                                                                                                                                              <  1:70021/251‑1 (MQ=255)
                                                                                                   ggcACCGGTGGGCCGATGCAGGTCAGCATGGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGTTGTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAgcgg                                                                                                                                         <  2:351633/250‑1 (MQ=255)
                                                                                                           ggggCCGATGAAGGTCAGCATGGCTTTGGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGGACAGTCCGTCAGGAAGTTTTCATTCGCCGCGGCGGGTTGTGGTGTGGTTCTTACCATTTACTGAATTATTCCGCCAGTGATAGCGCACCGATATACCTGTTTGCTGATTTTAACGCTGGCTGGGACGCCAGCCGTAATGCCGCGTTGCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGc                                                                                                                                 <  2:368726/250‑1 (MQ=255)
                                                                                                                    gCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCCTTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCCGTAAGGCCAgcggcg                                                                                                                                       >  1:243739/1‑236 (MQ=255)
                                                                                                                    gCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAgcggcg                                                                                                                                       <  2:243739/236‑1 (MQ=255)
                                                                                                                      aGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCCGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGCTGGcccttgac                                                                                                                      >  2:108500/1‑244 (MQ=255)
                                                                                                                          cAGCATTGCTTTTGCCGAACAGCAGACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCGGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGCGGGCACTGGACGGTg                                                                                                                  <  2:267201/251‑1 (MQ=255)
                                                                                                                                      tGCCGAACAGCATACCAACGGGGATCCGTGGAAAATGGACGGTACAGTCCGTCAGGACGTCTTCAGTCGCCGCGGCCGGTTGTGGGTTGGGACTTTCCCTTTACTGAATTATTCCGCCAGGTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTCCGCCAGCCCGAATACCGCGTTTCCGAACGCGGGCAGTAAGGCCAGCGGCGTTAAGCTGGCACTGGTCGGAGATTTAATTCGtt                                                                                                      >  2:123579/1‑251 (MQ=255)
                                                                                                                                         ggagcAGCATACCAAGGGGTAGCCGTGGAAAATGTACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGGTGTGGTTTGGTACTGACCATTTACTGAATTATTCCTCCAGTTATAGCGCACCGATATTCCGTTTTGCTGATTTTAACGCTGGGTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCTGTCAGTAAGGCCAGCGGCGTGAAGCTGGCACTGGACGGTGATTTAATTCGTTACg                                                                                                   <  2:259126/247‑1 (MQ=255)
                                                                                                                                           aaCAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGCTGGCACTGGACGGGGATTTAATTCGTTACGAc                                                                                                 >  1:10562/1‑251 (MQ=255)
                                                                                                                                                  tACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGCTGGCACTGGACGGTGATTTAATTCGTTACGACAGTaa                                                                                            >  1:127732/1‑249 (MQ=255)
                                                                                                                                                    ccAAAGGGTATCCGTGGAAAATGGACGGTTCAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGCTGGCACTGGACGGTGATTTAATTCGTTACGACAGTAAAGaa                                                                                        <  2:146864/251‑1 (MQ=255)
                                                                                                                                                       aaGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGCTGGCACTGGACGGTGATTTAATTCGTTACGACAGTAAAGAAccc                                                                                     <  2:10562/251‑1 (MQ=255)
                                                                                                                                                                          ggACGGAACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTACGCCAGTTATAGCGCACCGATATAACGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCCGAACGCGGTCCGTTAGGCCAGCGGCGTGAAGCTGGCACTGGACGGGGATTTAATTTGTTACCACAGTAAAAAACCCGGGAAACACGAACTGGCa                                                                   >  1:393099/1‑250 (MQ=255)
                                                                                                                                                                              ggTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCCGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCCGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCCGCGGCGTGAAGCTGGCACTGGACGGTGGTTTAATTCGTTACGACAGGAAAGAACCCGGGAAAACCGAAATGGCAAcgcg                                                              >  2:158535/1‑251 (MQ=255)
                                                                                                                                                                                    ggccGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGCTGGCACTGGACGGTGATTTAATTCGTTACGACAGTAAAGAACCCGGGAAAACCGAACTGGCAACGCGCAAACt                                                        <  2:127732/249‑1 (MQ=255)
                                                                                                                                                                                          aaggaatgtgtCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCGTTTACTGAATTATGCCGCCAGTTATAGCCCACCGAAATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGCTGGCACTGGACGGTGATTTAATTCGTTACTACAGTAAAGAACCCGGGAGAACCGAGCTGGCAACGCGCAAACTAGCTGc                                                  <  2:264464/241‑1 (MQ=255)
                                                                                                                                                                                               aGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGCTGGCACTGGACGGTGATTTAATTCGTTACGACAGTAAAGAACCCGGGAAAAACGAACTGGCAACGCGCAAACTGGCTGCaaaa                                              >  2:41489/1‑250 (MQ=255)
                                                                                                                                                                                                  cTTCAGTCGCCGCGGCGGGTTGTGGGTTGGTAATTACCATTTACTGGATTATTCCGCCAGTTATAGCGCACCGATATACCGTTTTTCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAAAGCGGTCAGTAAGGCCAGCGGCGTGAAGCTGGCACTGGACCGTGATTTAATTCGTTACGACAGTAAAAAACCCGGGAAAACCGAACTGGCAACCCGCAAAATGGGTGCcaaaattg                                          >  2:211293/1‑247 (MQ=255)
                                                                                                                                                                                                      aGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGCTGGCACTGGACGGTGATTTAATTCGTTACGACAGTAAAGAACCCGGGAAAACCGAACTGGCAACGCGCAAACTGGCTGCAAAACTGGGGAt                                      <  1:158535/251‑1 (MQ=255)
                                                                                                                                                                                                                      tggggTTTGGGACTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGCTGGCACTGGACGGTGATTTCATTCGTTACGACAGTAAAGAACCCGGGAAAAACGAACTGGCAACGCGCAAACTGGCTGCAAAACTGGGGATGAGCGACAGCGAAAtt                      >  2:18133/4‑251 (MQ=255)
                                                                                                                                                                                                                              ggttcTTACCATTTACTGAATTATTCCGCCAGTTATAGCGCACGGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGCTGGCACTGGACGGTGATTTAATTCGTTACGACAGTAAAGAACCCGGGAAAACCGAACTGGCAACGCGCAAACTGGCTGCAAAACTGGGGATGAGCGACAGCGAAATTCGCCGTCa              <  2:113807/247‑1 (MQ=255)
                                                                                                                                                                                                                               gTACTTACCATTTACTGAATTATCCCCCCGTTTATCGCGCACGGATATACCGTTTTGCTGATTTTAACGTGGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGCTGGCACTGGACGGTGATTTAATTCGTTACGACAGTAAAGAACCCGGGAAAACCGAACTGGCAACGCGCAAACTGGCTGCAAAACTGGGGATGAGCGACAGCGAAATTCGCCGTCAg             <  2:221610/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ttACCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGCTGGCACTGGACGGTGATTTAATTCGTTACGACAGTAAAGAACCCGGGAAAACCGAACTGGCAACGCGCAAACTGGCTGCAAAACTGGGGATGAGCGACAGCGAAATTCGCCGTCAGtt           >  1:156164/1‑249 (MQ=255)
                                                                                                                                                                                                                                     aCCATTTACTGAATTATTCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGCTGGCACTGGACGGTGATTTAATTCGTTACGACAGTAAAGAACCCGGGAAAACCGAACTGGCAACGCGCAAACTGGCTGCAAAACTGGGGATGAGCGACAGCGAAATTCGCCGTCAGTTag         >  2:357486/1‑249 (MQ=255)
                                                                                                                                                                                                                                           tACTGAATTATTCCGCCAGTTATAGCGCACCGATATTCCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGCTGGCACTGGACGGTGATTTAATTCGTTACGACAGTAAAGAACCCGGGAAAACCGAACTGGCAACGCGCAAACTGGCTGCAAAACTGGGGATGAGCGACAGCGAAATTCGCCGTCAGTTAGAGAAAgg  <  2:169850/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
GATTCGGTGCGTACAGAAAAGCAATTGAGCGCGATTTTTGACGATCTGATCAACATGGTGCCAATGGGGCAGACGCTGTTTGGTTCGCTCAATCCGGTGCGCACCGGTGGGCCGATGCAGGTCAGCATTGCTTTTGCCGAACAGCATACCAAAGGGTATCCGTGGAAAATGGACGGTACAGTCCGTCAGGAAGTCTTCAGTCGCCGCGGCGGGTTGTGGTTTGGTACTTACCATTTACTGAATTATCCCGCCAGTTATAGCGCACCGATATACCGTTTTGCTGATTTTAACGCTGGCTGGTACGCCAGCCGTAATGCCGCGTTTCAGAACGCGGTCAGTAAGGCCAGCGGCGTGAAGCTGGCACTGGACGGTGATTTAATTCGTTACGACAGTAAAGAACCCGGGAAAACCGAACTGGCAACGCGCAAACTGGCTGCAAAACTGGGGATGAGCGACAGCGAAATTCGCCGTCAGTTAGAGAAAGG  >  NC_000913/398286‑398770

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