Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,232 C→T 7.4% T370T (ACC→ACT icd → isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,196,2320CT 7.4% 248.8 / 5.5 80T370T (ACC→ACTicdisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase
Reads supporting (aligned to +/- strand):  ref base C (35/39);  new base T (4/2);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCG  >  NC_000913/1196084‑1196378
                                                                                                                                                    |                                                                                                                                                  
gtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGAcc                                                                                                                                                    >  5:242983/1‑149 (MQ=255)
   tcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCg                                                                                                                                                   >  8:154282/1‑147 (MQ=255)
   tcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCg                                                                                                                                                   >  6:258210/1‑147 (MQ=255)
     gCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGc                                                                                                                                               >  7:274740/1‑149 (MQ=255)
            gTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGa                                                                                                                                         <  7:154282/148‑1 (MQ=255)
                aaaCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAcctga                                                                                                                                    >  7:324101/1‑145 (MQ=255)
                 aaCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAAt                                                                                                                                   <  6:242983/149‑1 (MQ=255)
                 aaCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGAGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAAt                                                                                                                                   >  3:40580/1‑149 (MQ=255)
                  aCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAAtt                                                                                                                                  <  5:230669/149‑1 (MQ=255)
                  aCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGACGATGCCGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAAtt                                                                                                                                  <  3:289005/149‑1 (MQ=255)
                    aTCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGt                                                                                                                                <  4:172117/149‑1 (MQ=255)
                    aTCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGt                                                                                                                                >  2:381862/1‑149 (MQ=255)
                    aTCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGGCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGt                                                                                                                                <  6:169352/149‑1 (MQ=255)
                         tGACGAATGCGCCCTGTTTGAAGCCACCCACGTTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTTTTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAg                                                                                                                           <  8:324101/149‑1 (MQ=38)
                         tGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAg                                                                                                                           <  1:381862/149‑1 (MQ=255)
                             gAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTAt                                                                                                                       <  5:80835/149‑1 (MQ=255)
                              aaTGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGGAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATg                                                                                                                      <  8:274740/149‑1 (MQ=38)
                                  cgcCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGaa                                                                                                                   <  5:258210/148‑1 (MQ=255)
                                   gcCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAgg                                                                                                                 >  3:255410/1‑149 (MQ=255)
                                   gcCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCt                                                                                                                                           >  3:46604/1‑123 (MQ=255)
                                   gcCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCt                                                                                                                                           <  4:46604/123‑1 (MQ=255)
                                     ccTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTaa                                                                                                                                    <  4:352436/128‑1 (MQ=255)
                                     ccTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTaa                                                                                                                                    >  3:352436/1‑128 (MQ=255)
                                      cggggtGACGCCGCCCACGGTACGGCGCCGAAATATGCCGGTCAGGACAAAGTAAGGGCGGGGGCTATTGTTCTCTCCGCTGAGATGATGCTGCGCGACATGGGGTGGACCGAAGCGGCTGCCGGAATTGGTAAAGGTATGGTAGgcgg                                                                                                              >  2:154469/6‑148 (MQ=255)
                                           tGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTAATCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATc                                                                                                          <  1:313408/148‑1 (MQ=255)
                                              aGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcg                                                                                                      <  4:40580/149‑1 (MQ=255)
                                              aGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCGCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGCCTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcg                                                                                                      >  7:383809/1‑149 (MQ=255)
                                               gCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaatg                                                                                                      >  1:191439/1‑146 (MQ=255)
                                                  acccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTaaa                                                                                                                            <  2:333851/123‑1 (MQ=255)
                                                  acccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTaaa                                                                                                                            >  1:333851/1‑123 (MQ=255)
                                                   cccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaa                                                                                                  >  3:138329/1‑148 (MQ=255)
                                                      acGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAcc                                                                                              >  3:49014/1‑149 (MQ=255)
                                                      acGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAcc                                                                                              >  7:56504/1‑149 (MQ=255)
                                                      acGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTGAAAGGTATGGAAGGCGCAATCAACGCGAAAAcc                                                                                              >  4:233338/1‑149 (MQ=255)
                                                        cgTACTGCGCCGAAATATGCCCGTCATGACAAAGTAAATCCTTGCTCTATTATTCTCTCCGCTTAGATGATGCTGGGCCACATGGGTTGGACCGAAGCGGCTTACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGTAAACCGt                                                                                            <  6:352735/148‑1 (MQ=37)
                                                          tACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGTACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa                                                                                          <  2:171154/149‑1 (MQ=255)
                                                               cgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTAt                                                                                     <  8:56504/149‑1 (MQ=255)
                                                               cgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTAt                                                                                     <  3:233338/149‑1 (MQ=255)
                                                                gcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAACACCGTAACCTAt                                                                                     >  8:357546/1‑148 (MQ=255)
                                                                  cGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGAc                                                                                  >  4:179062/1‑149 (MQ=255)
                                                                      atatGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                                                                                                       >  6:212465/1‑124 (MQ=255)
                                                                      atatGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                                                                                                       <  5:212465/124‑1 (MQ=255)
                                                                      atatGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa                                                                                          >  1:205906/1‑137 (MQ=255)
                                                                      atatGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa                                                                                          <  2:205906/137‑1 (MQ=255)
                                                                                 aGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTATCCTATGACTTCGAGCGTCTgat                                                                    <  4:159046/148‑1 (MQ=255)
                                                                                   gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGACGCGGCTGACTTAATTGTTACCGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatgga                                                                 <  8:34164/149‑1 (MQ=255)
                                                                                   gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatgga                                                                 >  7:21084/1‑149 (MQ=255)
                                                                                         gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGgcgc                                                           >  6:80789/1‑149 (MQ=255)
                                                                                             aTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCTCATGGGGTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTaaa                                                       <  8:383809/149‑1 (MQ=255)
                                                                                                cTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGATAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGCGCGTCTGATGGATGCCGCTAAActc                                                    >  2:331690/1‑148 (MQ=255)
                                                                                                   gCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGa                                                                             >  2:369600/1‑121 (MQ=255)
                                                                                                   gCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGa                                                                             <  1:369600/121‑1 (MQ=255)
                                                                                                       tattatTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCGACGCGAAAACCGTAACCTATGACTTCGAGCGt                                                                         >  3:108616/1‑121 (MQ=255)
                                                                                                       tattatTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCGACGCGAAAACCGTAACCTATGACTTCGAGCGt                                                                         <  4:108616/121‑1 (MQ=255)
                                                                                                        attattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTAt                                                                                     <  8:258427/108‑1 (MQ=255)
                                                                                                        attattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTAt                                                                                     >  7:258427/1‑108 (MQ=255)
                                                                                                         ttattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGt                                           >  2:226861/1‑149 (MQ=255)
                                                                                                           attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgat                                                                    <  3:205387/122‑1 (MQ=255)
                                                                                                           attCTCTCAGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgat                                                                    >  4:205387/1‑122 (MQ=255)
                                                                                                            ttCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCa                                        <  4:138329/149‑1 (MQ=255)
                                                                                                               tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCagag                                     >  5:144949/1‑149 (MQ=255)
                                                                                                                 tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTAt                                                                                     <  5:191603/99‑1 (MQ=255)
                                                                                                                 tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTAt                                                                                     >  6:191603/1‑99 (MQ=255)
                                                                                                                     cTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGt                               <  3:179062/149‑1 (MQ=255)
                                                                                                                      tGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa                                                                                          <  6:326164/89‑1 (MQ=255)
                                                                                                                      tGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa                                                                                          >  5:326164/1‑89 (MQ=255)
                                                                                                                         gatgatGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcg                           >  2:206561/1‑149 (MQ=255)
                                                                                                                          atgatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcgc                          >  7:217157/1‑149 (MQ=255)
                                                                                                                             atgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGat                       <  1:226861/149‑1 (MQ=255)
                                                                                                                                   cgcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatcatc                   <  6:144949/147‑1 (MQ=255)
                                                                                                                                   cgcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatcatc                   <  8:21084/147‑1 (MQ=255)
                                                                                                                                    gcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatggatg                                                               >  3:132359/1‑102 (MQ=255)
                                                                                                                                    gcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatggatg                                                               <  4:132359/102‑1 (MQ=255)
                                                                                                                                      cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcgc                          <  5:16586/137‑1 (MQ=255)
                                                                                                                                      cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcgc                          >  6:16586/1‑137 (MQ=255)
                                                                                                                                         atttgTTGGACCGAAGCGGCTGACTTAATTGTAAAAGGTATGGAAGCCGCAATCAACGACAAAACAGTAACATAAGACTTCGAGAGTCTGATGGATGGCGCTAAACTCCTGAAATCTTCACAGTTTGGTGACGCGATCATCGAAAAAAt           >  2:188818/5‑149 (MQ=255)
                                                                                                                                            ggTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGt         <  1:331690/148‑1 (MQ=255)
                                                                                                                                              ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAAt      >  8:146277/1‑149 (MQ=255)
                                                                                                                                                 gACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTa                                                                                           >  1:33163/1‑61 (MQ=255)
                                                                                                                                                 gACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTa                                                                                           <  2:33163/61‑1 (MQ=255)
                                                                                                                                                  aCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCg  <  1:206561/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCG  >  NC_000913/1196084‑1196378

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