Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,974,082 T→A 8.2% K523* (AAG→TAG) ‡ flgK ← flagellar hook‑associated protein FlgK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,974,0820TA 8.2% 226.3 / 7.9 73K523* (AAG→TAG) ‡flgKflagellar hook‑associated protein FlgK
Reads supporting (aligned to +/- strand):  ref base T (32/35);  new base A (3/3);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTGCG  >  NC_002947/4973936‑4974227
                                                                                                                                                  |                                                                                                                                                 
gTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACttgtt                                                                                                                                                >  2:894802/1‑150 (MQ=255)
   ggCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTc                                                                                                                                             <  2:1408596/150‑1 (MQ=255)
   ggCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTc                                                                                                                                             <  1:644559/150‑1 (MQ=255)
   ggCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTc                                                                                                                                             <  1:1321611/150‑1 (MQ=255)
   ggCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACAAGTTGTc                                                                                                                                             <  2:1030632/150‑1 (MQ=255)
     cATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACAAGTTGTcca                                                                                                                                           <  1:1353901/150‑2 (MQ=255)
           aCGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAg                                                                                                                                     <  1:911777/150‑1 (MQ=255)
              aCAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCt                                                                                                                                           >  2:271162/1‑141 (MQ=255)
              aCAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCt                                                                                                                                           <  1:271162/141‑1 (MQ=255)
                    tAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACttgtt                                                                                                                                                <  1:1032254/130‑1 (MQ=255)
                    tAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACttgtt                                                                                                                                                >  2:1032254/1‑130 (MQ=255)
                    tAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGc                                                                                                                            >  2:640397/1‑150 (MQ=255)
                    tAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGc                                                                                                                            <  1:455540/150‑1 (MQ=255)
                            tCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTa                                                                                                                    >  1:203983/1‑150 (MQ=255)
                            tCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTa                                                                                                                    <  1:640397/150‑1 (MQ=255)
                             cATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTAc                                                                                                                   <  1:554737/150‑1 (MQ=255)
                             cATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTAc                                                                                                                   <  2:1127753/150‑1 (MQ=255)
                              aTTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCCCCACCCTGCCGGTAcc                                                                                                                  >  2:483938/1‑150 (MQ=255)
                              aTTGGCCGCAGGGATGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTAcc                                                                                                                  >  2:1059442/1‑150 (MQ=255)
                                     gCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTgg                                                                                                           <  1:894802/150‑1 (MQ=255)
                                        gggCTGCCACCGATTCTGGTCTCGACGGCAAAGGTACGCGGCTTCCTTGTATTATCCAGTATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGCACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTa                                                                                                        <  2:37302/150‑1 (MQ=255)
                                             gCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCtt                                                                                                   <  2:203983/150‑1 (MQ=255)
                                              ccACCGATTGTGGTCTCGGCGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTg                                                                                                  >  1:321118/1‑150 (MQ=255)
                                              ccACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTg                                                                                                  <  2:1104610/150‑1 (MQ=255)
                                                       gtgGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAAc                                                                                         <  1:778122/150‑1 (MQ=255)
                                                           tctcGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCt                                                                                                             <  2:1412235/126‑1 (MQ=255)
                                                           tctcGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCt                                                                                                             >  1:1412235/1‑126 (MQ=255)
                                                            ctcGACGGAAAAGGTTGGTGTTTTTCCTCTACCATCCAGAAGCGGGTCGCCGTCCGCATCTAGCATCGGCACTTCGATGGACAACTTGCTTTCCTTCTCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCtt                                                                                    <  2:275961/150‑1 (MQ=255)
                                                             tcGACGTCAAAGGGACGCGGGTTACGGCTTGCATCCAGAGTCGGGTCTCCGGCCGCCTCCAGGACCGGCACGACTCTGGACAACTAGTTGTCCTGCCCAGGGACAACGCTGACGGTACCAACCTGGGTACCCTTGGCGTCGAACAGCttt                                                                                   >  2:239994/1‑149 (MQ=255)
                                                             tcGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTg                                                                                   >  2:985372/1‑150 (MQ=255)
                                                               gACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTa                                                                                 >  1:1008149/1‑150 (MQ=255)
                                                                aCGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAg                                                                                <  2:1202639/150‑1 (MQ=255)
                                                                   gCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCCACCAGCATGGGCACGGCCATGGGCAACATGTTGTCCTGCCAAGGGACCGCGCTGCCAGGACCAATCTGTGTCCCCTTGACGTCGATCAGCTTGTAGccc                                                                             >  2:499618/1‑150 (MQ=255)
                                                                     aaaGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTg                                                                           >  1:948519/1‑150 (MQ=255)
                                                                        ggTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGa                                                                        >  2:1360130/1‑150 (MQ=255)
                                                                         gTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTa                                                                                 >  2:1321902/1‑140 (MQ=255)
                                                                         gTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTa                                                                                 <  1:1321902/140‑1 (MQ=255)
                                                                               gggTTACCGCTACCATCCAGAATCTGGTCGCCGGCAGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGGCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGGGCCACCg                                                                 >  1:472869/1‑150 (MQ=255)
                                                                                  ttACCGCTACCATCCAGCATCGGGCCGCCGGCAGCAACCCGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTg                                                              <  2:1008149/150‑1 (MQ=255)
                                                                                         tACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGt                                                       >  2:1360335/1‑150 (MQ=255)
                                                                                             aTCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTa                                                                                                                    >  1:578877/1‑85 (MQ=255)
                                                                                             aTCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTa                                                                                                                    <  2:578877/85‑1 (MQ=255)
                                                                                              tCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGaa                                                                                          <  2:408197/110‑1 (MQ=255)
                                                                                              tCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGaa                                                                                          >  1:408197/1‑110 (MQ=255)
                                                                                              tCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTCTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAac                                                  <  1:1360130/150‑1 (MQ=255)
                                                                                                 aGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATATGGGTCCCCGTGTCGTGGAACAGCTTGTAGCTCTGGGAGCCACCGCTGGCCGCGTCGAAcacc                                               >  2:336992/1‑150 (MQ=255)
                                                                                                   aaTCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAAcaccac                                             <  1:1059442/150‑1 (MQ=255)
                                                                                                      tccatcgccgtgccgatgctggatgcggccggcgaTGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGATCCACCGCTGGCAGCGTCGAACACCAcgcg                                          >  2:254664‑M2/36‑150 (MQ=255)
                                                                                                       ccatcgccgtgccgatgctggatgcggccggcgaTGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTg                                                              <  1:1081614‑M2/95‑1 (MQ=255)
                                                                                                       ccatcgccgtgccgatgctggatgcggccggcgaTGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTg                                                              >  2:1081614‑M2/35‑129 (MQ=255)
                                                                                                         atcgccgtgccgatgctggatgcggccggcgaTGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTg                                                              >  2:1291429‑M2/33‑127 (MQ=255)
                                                                                                         atcgccgtgccgatgctggatgcggccggcgaTGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTg                                                              <  1:1291429‑M2/95‑1 (MQ=255)
                                                                                                           CGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCg                                     >  2:948519‑M1/1‑150 (MQ=255)
                                                                                                               GGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGACACCTCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGGCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCAt                                 >  1:1155382‑M1/1‑150 (MQ=255)
                                                                                                                gccgCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATc                                >  1:517912/1‑150 (MQ=255)
                                                                                                                gccgCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTAGAACACCACGCGTACCGGCATc                                >  2:617654/1‑150 (MQ=255)
                                                                                                                  cgCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTgg                                                                                                           <  1:796256/73‑1 (MQ=255)
                                                                                                                  cgCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTgg                                                                                                           >  2:796256/1‑73 (MQ=255)
                                                                                                                   gCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaa                             <  2:321118/150‑1 (MQ=255)
                                                                                                                    cATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGaa                                                                                          >  1:516377/1‑88 (MQ=255)
                                                                                                                    cATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGaa                                                                                          <  2:516377/88‑1 (MQ=255)
                                                                                                                        cAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaatcgcg                        <  1:1103369/150‑1 (MQ=255)
                                                                                                                         aGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaatcgcga                       >  1:1167469/1‑150 (MQ=255)
                                                                                                                          gCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaatcgcgaa                      >  1:1109276/1‑150 (MQ=255)
                                                                                                                           cATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaatcgcgaat                     <  1:617654/150‑1 (MQ=255)
                                                                                                                            aTCGGAACGGAGATGGACAACTTGTTGTCCTGCCCAGGCACCACGATGCCGGTAGCAATCTGGGTCCCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaatcgcgaatc                    <  2:293188/150‑1 (MQ=255)
                                                                                                                                gCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGa                                                                        <  2:717481/94‑1 (MQ=255)
                                                                                                                                gCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGa                                                                        >  1:717481/1‑94 (MQ=255)
                                                                                                                                         aggACAACAAGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCa                                                           >  1:25/2‑98 (MQ=255)
                                                                                                                                         aggACAACAAGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCa                                                           <  2:25/97‑1 (MQ=255)
                                                                                                                                         aggACAACAAGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCa       >  2:832261/2‑150 (MQ=255)
                                                                                                                                          ggACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGActgctg            <  2:531680/144‑1 (MQ=255)
                                                                                                                                          ggACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGActgctg            >  1:531680/1‑144 (MQ=255)
                                                                                                                                            aCAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAAcaccac                                             <  2:411071/109‑1 (MQ=255)
                                                                                                                                            aCAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAAcaccac                                             >  1:411071/1‑109 (MQ=255)
                                                                                                                                              aaCTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTGCg  <  2:1109276/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTGCG  >  NC_002947/4973936‑4974227

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