Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,726,162 T→C 7.2% Y165C (TAC→TGC) ‡ PP_2387 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,726,1620TC 7.2% 195.6 / 4.3 69Y165C (TAC→TGC) ‡PP_2387hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/37);  new base C (3/2);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CGGGGGCGACAGTGGCCAGCTGACCGGGGGCGCGGCCACAGGAGGGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGG  >  NC_002947/2726023‑2726308
                                                                                                                                           |                                                                                                                                                  
cGGGGGCGACAGTGGCCAGCTGACCGGGGGCGCGGCCACAGGAGGGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGAAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTg                                                                                                                                            <  2:184654/148‑1 (MQ=255)
 gggggCGACAGTGGCCAGCTGACCGGGGGCGCGGCCACAGGAGGGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGc                                                                                                                                          >  1:2741/1‑149 (MQ=255)
     gCGACAGTGGCCAGCTGACCGGGGGCGCGGCCACAGGAGGGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGtgc                                                                                                                                      <  2:2741/149‑1 (MQ=255)
        aCAGTGGCCAGCTGACCGGGGGCGCGGCCACAGGAGGGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGtgctgc                                                                                                                                   >  6:82629/1‑149 (MQ=255)
          agcggCCAGCTGACCGGGGGCGCGGCCACAGGAGGGGTCAGCGCTGGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCgg                                                                                                                                 <  4:222162/146‑1 (MQ=255)
              gCCAGCTGACCGGGGGCGCGGCCACAGGAGGGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCg                                                                                                                             >  2:279301/1‑149 (MQ=255)
               ccAGCTGCCCGGGCGCGCGGCCACAGTCGGGGTCCGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGt                                                                                                                            <  8:252800/149‑1 (MQ=255)
               ccAGCTGACCGGGGGCGCGGCCACAGGAGGGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGt                                                                                                                            <  3:275169/149‑1 (MQ=255)
                cAGCTGACCGGGGGCGCGGCCACAGGAGGGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGATGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGtt                                                                                                                           <  2:286534/149‑1 (MQ=255)
                  gCTGACCGGGGGCGCGGCCACAGGAGGGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTgg                                                                                                                         <  5:161419/149‑1 (MQ=255)
                         gggggCGCGGCCACAGGAGGGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGc                                                                                                                  <  1:279301/149‑1 (MQ=255)
                         gggggCGCGGCCACAGGAGGGGTCAGCGCTTGTCCGAGGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGc                                                                                                                  <  4:272455/149‑1 (MQ=255)
                           gggCGCGGCCACAGGAGGGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGa                                                                                                                <  2:401034/149‑1 (MQ=255)
                           gggCGCGGCCACAGGAGGGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGa                                                                                                                <  5:158063/149‑1 (MQ=255)
                           gggCGCGGCCACAGGAGGGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGa                                                                                                                <  3:13889/149‑1 (MQ=255)
                               gcgGCCACAGGAGGGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGAc                                                                                                            <  3:32779/149‑1 (MQ=255)
                               gcgGCCACAGGAGGGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCATCGCTAAAACAGTTAAGAAAGGCCTGCTAGGCCTTGGCGTGCTGAGGGTCGTTGGTACCTGCGAGGAc                                                                                                            >  1:331715/1‑149 (MQ=255)
                                   ccACAGGAGGGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTca                                                                                                        <  3:223575/149‑1 (MQ=255)
                                       aggaggGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTcacaca                                                                                                    <  3:330702/149‑1 (MQ=255)
                                       aggaggGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTcaca                                                                                                      >  2:358485/1‑147 (MQ=255)
                                       aggaggGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTcaca                                                                                                      >  8:2458/1‑147 (MQ=255)
                                            gggTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCt                                                                                               >  2:353333/1‑149 (MQ=255)
                                              gTCAGCTCTTGTCCGATGCGCCGCCACAGAGACTGCCGACGAGGTGTTGCGCAGGGCTGCATGCATTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACATGGCTTc                                                                                             <  2:331715/149‑1 (MQ=255)
                                              gTCAGCGCTTGTCCTCTGCGCCGCCACAGCGACTGCCGACGCGGTGGTTCGCAGGGCTGCATGCCTTCTCTAAAACTGTTGAGAAAGGCCTGGTAGTCCTTGGCGTGCTGCGGGTCTTTGGTACCTGCGAGGACTTCACACAGGGCTTc                                                                                             <  7:277520/149‑1 (MQ=255)
                                              gTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCGTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTc                                                                                             >  6:265597/1‑149 (MQ=255)
                                              gTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTc                                                                                             >  6:54453/1‑149 (MQ=255)
                                              gTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTc                                                                                             >  7:150616/1‑149 (MQ=255)
                                                     cTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGcc                                                                                      <  5:301643/149‑1 (MQ=255)
                                                            gATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAACGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCCCACAGGGCTTCGTGGGCCAGGCTGt                                                                               >  7:276560/1‑149 (MQ=255)
                                                            gATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAgg                                                                                   >  8:281155/1‑145 (MQ=255)
                                                            gATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAgg                                                                                   <  7:281155/145‑1 (MQ=255)
                                                             aTGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTTCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCa                                                                                     <  1:326262/142‑1 (MQ=255)
                                                             aTGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTTCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCa                                                                                     >  2:326262/1‑142 (MQ=255)
                                                                cgccgcCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGa                                                                            <  7:2458/148‑1 (MQ=255)
                                                                cgccgcCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGa                                                                            <  3:117862/148‑1 (MQ=255)
                                                                cgccgcCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTGGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGa                                                                            <  1:90244/148‑1 (MQ=255)
                                                                        caGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGa                                                                   <  3:153892/149‑1 (MQ=255)
                                                                            gACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGcac                                                               >  4:175162/1‑149 (MQ=255)
                                                                              cTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGcacca                                                             >  4:149475/1‑149 (MQ=255)
                                                                                     cgtggtggtGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATg                                                      <  1:149181/146‑1 (MQ=255)
                                                                                          tggtgCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGa                                                 <  3:175162/149‑1 (MQ=255)
                                                                                            gtgCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCGCACAGGGCTTc                                                                                             >  8:205817/1‑103 (MQ=255)
                                                                                            gtgCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCGCACAGGGCTTc                                                                                             <  7:205817/103‑1 (MQ=255)
                                                                                               cgcAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTgg                                            <  5:265597/149‑1 (MQ=255)
                                                                                               cgcAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTgg                                            >  6:129456/1‑149 (MQ=255)
                                                                                                gcAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGc                                           <  1:48114/149‑1 (MQ=255)
                                                                                                      cTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAgg                                                                                   <  3:27212/103‑1 (MQ=255)
                                                                                                      cTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAgg                                                                                   >  4:27212/1‑103 (MQ=255)
                                                                                                         cATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGt                                  >  8:86005/1‑149 (MQ=255)
                                                                                                         cATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGt                                  >  1:102838/1‑149 (MQ=255)
                                                                                                            gCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCt                               >  3:92512/1‑149 (MQ=255)
                                                                                                            gCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCt                               >  7:341851/1‑149 (MQ=255)
                                                                                                             ccTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCtg                              <  5:54453/149‑1 (MQ=255)
                                                                                                                 cGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCtggtg                           <  5:156688/148‑1 (MQ=255)
                                                                                                                      aaaCTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCAc                                                                      >  6:202918/1‑100 (MQ=255)
                                                                                                                      aaaCTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCAc                                                                      <  5:202918/100‑1 (MQ=255)
                                                                                                                      aaaCTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCa                      <  1:358485/148‑1 (MQ=255)
                                                                                                                          tGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCtggt                            >  6:44950/1‑138 (MQ=255)
                                                                                                                          tGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCtggt                            <  5:44950/138‑1 (MQ=255)
                                                                                                                          tGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGaa                  >  1:324089/1‑148 (MQ=255)
                                                                                                                            ttGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAAca                >  6:95919/1‑148 (MQ=255)
                                                                                                                            ttGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAAc                 <  3:149475/147‑1 (MQ=255)
                                                                                                                            ttGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAAc                 <  4:92512/147‑1 (MQ=255)
                                                                                                                               gccaaGGCCTACCAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAAcatcat            >  6:252364/4‑149 (MQ=255)
                                                                                                                                ccaaGGCCTACCAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAAcatca             >  2:297684/3‑147 (MQ=255)
                                                                                                                                ccaaGGCCTACCAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAAcatca             <  1:297684/145‑1 (MQ=255)
                                                                                                                                 caaGGCCTACCAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCtggt                            <  5:407922/130‑1 (MQ=255)
                                                                                                                                 caaGGCCTACCAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCtggt                            >  6:407922/2‑131 (MQ=255)
                                                                                                                                         ggTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTgg  >  7:148219/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
CGGGGGCGACAGTGGCCAGCTGACCGGGGGCGCGGCCACAGGAGGGGTCAGCGCTTGTCCGATGCGCCGCCACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGG  >  NC_002947/2726023‑2726308

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: ‑