Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 346,120 T→C 5.7% N576D (AAC→GAC) ‡ PP_0287 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947346,1200TC 5.7% 262.7 / 4.7 87N576D (AAC→GAC) ‡PP_0287hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (37/45);  new base C (3/2);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.16e-01

ACGCCATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTT  >  NC_002947/345972‑346260
                                                                                                                                                    |                                                                                                                                              
aCGCCATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGtt                                                                                                                                                <  2:110733/149‑1 (MQ=255)
   ccATTGCGTACCACCAGGCTGCCTCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGAt                                                                                                                                             <  2:319488/149‑1 (MQ=255)
    cATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATg                                                                                                                                            >  7:33139/1‑149 (MQ=255)
    cATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGGGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATg                                                                                                                                            >  8:347494/1‑149 (MQ=255)
     aTTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATgg                                                                                                                                           >  5:71336/1‑149 (MQ=255)
     aTTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATgg                                                                                                                                           >  7:5781/1‑149 (MQ=255)
     aTTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATgg                                                                                                                                           >  6:284549/1‑149 (MQ=255)
      ttGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGc                                                                                                                                          <  5:261355/149‑1 (MQ=255)
         cGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTTGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGaa                                                                                                                                       <  6:276913/149‑1 (MQ=255)
         cGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGaa                                                                                                                                       <  8:33139/149‑1 (MQ=255)
         cGTACCAACAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGAGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCg                                                                                                                                         >  4:136834/1‑147 (MQ=255)
          gTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAg                                                                                                                                      >  3:119630/1‑149 (MQ=255)
          gTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCAACGATGGCgtgc                                                                                                                                      <  5:325926/149‑4 (MQ=255)
           tACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGt                                                                                                                                     >  5:169166/1‑149 (MQ=255)
              caccaGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTg                                                                                                                                  >  7:218846/1‑149 (MQ=255)
               accaGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCTCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTgc                                                                                                                                 <  2:241371/149‑1 (MQ=255)
                           cgcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAg                                                                                                                     <  5:52326/149‑1 (MQ=255)
                           cgcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAg                                                                                                                     <  8:5781/149‑1 (MQ=255)
                            gcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTTACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                    <  8:68313/149‑1 (MQ=255)
                            gcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGCCCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                    <  1:76671/149‑1 (MQ=255)
                            gcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                    <  3:136834/149‑1 (MQ=255)
                            gcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                    <  1:359439/149‑1 (MQ=255)
                            gcAGCTCCTGGAAGGGGGTGTCCATGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                    >  7:233281/1‑149 (MQ=255)
                              aGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGgt                                                                                                                  <  8:218846/149‑1 (MQ=255)
                               gCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGgtg                                                                                                                 <  7:138838/149‑1 (MQ=255)
                                   cTGGAAGGGGGTGTGCTTTCCGCGTGGTTCGCGGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTc                                                                                                               <  4:285421/147‑1 (MQ=255)
                                   cTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGa                                                                                                             >  6:27817/1‑149 (MQ=255)
                                     ggAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAg                                                                                                                     >  7:1621/1‑139 (MQ=255)
                                     ggAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAg                                                                                                                     <  8:1621/139‑1 (MQ=255)
                                      gAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGTCGTTCAGGGTGTCGACCa                                                                                                          <  7:84854/149‑1 (MQ=255)
                                      gAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGAc                                                                                                            <  5:284549/147‑1 (MQ=255)
                                      gAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCa                                                                                                          <  4:119630/149‑1 (MQ=255)
                                        aGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGt                                                                                                        >  1:273703/1‑149 (MQ=255)
                                           gggTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCt                                                                                                     >  5:201395/1‑149 (MQ=255)
                                                tCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTg                                                                                              <  1:124864/149‑1 (MQ=255)
                                                 ccTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                    <  8:59991/128‑1 (MQ=255)
                                                  cTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTgg                                                                                             <  6:169166/148‑1 (MQ=255)
                                                  cTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCAACGATGGCgtgctg                                                                                                                                    >  7:27672/1‑106 (MQ=255)
                                                  cTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCAACGATGGCgtgctg                                                                                                                                    <  8:27672/111‑6 (MQ=255)
                                                     gCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGtt                                                                                                                        <  5:12101/120‑1 (MQ=255)
                                                     gCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGtt                                                                                                                        >  6:12101/1‑120 (MQ=255)
                                                     gCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGt                                                                                         >  3:88991/1‑149 (MQ=255)
                                                      ccGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGt                                                                                         >  2:258070/1‑148 (MQ=255)
                                                          gTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCg                                                                                    <  2:273703/149‑1 (MQ=255)
                                                           tGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCgg                                                                                   <  1:258070/149‑1 (MQ=255)
                                                             gTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTg                                                                                 <  8:233281/149‑1 (MQ=255)
                                                              ttCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTgg                                                                                >  4:135321/1‑149 (MQ=255)
                                                                  ctgctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACATCTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGt                                                                            <  7:347494/149‑1 (MQ=255)
                                                                  ctgctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGt                                                                            <  8:272137/149‑1 (MQ=255)
                                                                   tgctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGt                                                                            >  8:286454/1‑148 (MQ=255)
                                                                      tCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGt                                                                            >  8:132351/1‑145 (MQ=255)
                                                                      tCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGt                                                                            >  2:31364/1‑145 (MQ=255)
                                                                          cGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCAc                                                                    <  6:243376/149‑1 (MQ=255)
                                                                           gTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACt                                                                   <  7:286454/149‑1 (MQ=255)
                                                                           gTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACt                                                                   >  4:340733/1‑149 (MQ=255)
                                                                              ttGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTggt                                                                >  1:98398/1‑149 (MQ=255)
                                                                                  gATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTggtcagg                                                            <  5:27817/149‑1 (MQ=255)
                                                                                  gATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTggtcagg                                                            <  2:349871/149‑1 (MQ=255)
                                                                                   aTTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGt                                                                                  <  8:16316/126‑1 (MQ=255)
                                                                                   aTTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGt                                                                                  >  7:16316/1‑126 (MQ=255)
                                                                                   aTTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTggtcaggg                                                           >  7:291347/1‑149 (MQ=255)
                                                                                        gCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTggtc                                                               >  8:152508/1‑140 (MQ=255)
                                                                                        gCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTggtc                                                               <  7:152508/140‑1 (MQ=255)
                                                                                                     ttGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGt                                                                                  <  5:12291/108‑1 (MQ=255)
                                                                                                     ttGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGt                                                                                  >  6:12291/1‑108 (MQ=255)
                                                                                                     ttGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACcg                                         >  6:167612/1‑149 (MQ=255)
                                                                                                         cacaACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACcggcgg                                     <  6:71336/149‑1 (MQ=255)
                                                                                                          acaACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCggcggc                                    >  3:277029/1‑149 (MQ=255)
                                                                                                          acaACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCggcggc                                    >  3:361130/1‑149 (MQ=255)
                                                                                                            aaCTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTggtcagg                                                            >  6:172423/1‑123 (MQ=255)
                                                                                                            aaCTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTggtcagg                                                            <  5:172423/123‑1 (MQ=255)
                                                                                                                 ctgTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGtt                             >  8:284412/1‑149 (MQ=255)
                                                                                                                      ccAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCtt                         <  7:132351/148‑1 (MQ=255)
                                                                                                                      ccAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCtt                         <  4:88991/148‑1 (MQ=255)
                                                                                                                       cAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTg                        <  5:167612/148‑1 (MQ=255)
                                                                                                                             ggCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCg                  <  7:284412/148‑1 (MQ=255)
                                                                                                                               cATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTTAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTc               >  5:96964/1‑149 (MQ=255)
                                                                                                                                   gACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTtg            <  3:135321/148‑1 (MQ=255)
                                                                                                                                   gACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCATGGCCTTCGGTCTtga           >  8:77679/1‑149 (MQ=255)
                                                                                                                                    aCCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGc                                                  <  4:197264/109‑1 (MQ=255)
                                                                                                                                    aCCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGc                                                  >  3:197264/1‑109 (MQ=255)
                                                                                                                                    aCCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTtgat          <  5:74372/149‑1 (MQ=255)
                                                                                                                                     aaacttcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTGCTTTCCTTCGGTCTtgat          >  1:265751/7‑148 (MQ=255)
                                                                                                                                        ggaCGCCTTCGTTGTTGGGGAAGTTTGCGCTGCCGTTGAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATGAGTGGTCAGGGTGAGCAGGCCCTTGACCTGCTGCGTCTGTTCCTTGGCTTCGGTCTTGATGAAg      <  4:361130/147‑1 (MQ=255)
                                                                                                                                        cttcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCACGGGGTGTTACGGGTGGGTGTCAGATAAGTGGTCAGGGTCAGCAGGCCCTTGACCGGGGGCGTCGGTTCCTTGCGTGCGGTATttatga        >  1:106088/4‑149 (MQ=255)
                                                                                                                                           cGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTTAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAgtgt   <  6:96964/149‑1 (MQ=255)
                                                                                                                                            gCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGtt  >  1:360537/1‑149 (MQ=255)
                                                                                                                                            gCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCGTCGGTCTTGATGAAGTGtt  >  6:375543/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
ACGCCATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTT‑‑CTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTT  >  NC_002947/345972‑346260

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