Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,713,655 T→G 5.7% F254V (TTT→GTT)  srmB → ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,713,6550TG 5.7% 187.4 / 3.3 70F254V (TTT→GTT) srmBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base G (2/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

GGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGACCGAAGGTCGCGTAAACGTACTGGT  >  NC_000913/2713507‑2713803
                                                                                                                                                    |                                                                                                                                                    
ggAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTgtgt                                                                                                                                                      <  2:196104/149‑1 (MQ=255)
ggAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACACTTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCAGCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTgtgt                                                                                                                                                      <  2:1098230/149‑1 (MQ=255)
        cGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGt                                                                                                                                              >  1:725249/1‑149 (MQ=255)
           tGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGATCAATTGTGTTTGTGCGTAag                                                                                                                                           >  2:492821/1‑149 (MQ=255)
             gAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCGATTGTGTTTGTGCGTAagcg                                                                                                                                         >  1:1453001/1‑149 (MQ=255)
             gAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcg                                                                                                                                         >  2:1397190/1‑149 (MQ=255)
             gAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcg                                                                                                                                         <  1:1233990/149‑1 (MQ=255)
              aaGTTTCTGCCTATCCCTCCAGCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAATCCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcgt                                                                                                                                        <  1:72847/149‑1 (MQ=255)
                      gCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCgtgt                                                                                                                                >  1:717852/1‑149 (MQ=255)
                       ccAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGCGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGGTTGTGCGTAAGCGTGAGCGtgtg                                                                                                                               >  2:1016369/1‑149 (MQ=255)
                       ccAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGtgtg                                                                                                                               <  2:329073/149‑1 (MQ=255)
                       ccAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGtgtg                                                                                                                               >  2:640508/1‑149 (MQ=255)
                       ccAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGGTTGTGCGTAAGCGTGAGCGtgtg                                                                                                                               >  2:1021992/1‑149 (MQ=255)
                       ccAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCATGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGtgtg                                                                                                                               <  2:1041526/149‑1 (MQ=255)
                       ccAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGc                                                                                                                              >  1:671765/1‑149 (MQ=255)
                          aTCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCAt                                                                                                                            >  2:619393/1‑149 (MQ=255)
                          aTCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCAt                                                                                                                            >  1:863813/1‑149 (MQ=255)
                             ccTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACCGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAg                                                                                                                         >  2:106208/1‑149 (MQ=255)
                              cTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCTTTGCTGGTGCATATGTTAAAACAGCCGGAAGCGACCCGATCAATTGTGTTTGTGCGTAATAGTTAGCGTGTGCAAGAGc                                                                                                                        >  2:932953/1‑149 (MQ=255)
                                   cccGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCa                                                                                                                   >  2:724926/1‑149 (MQ=255)
                                   cccGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCa                                                                                                                   >  2:731019/1‑149 (MQ=255)
                                    ccGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCaa                                                                                                                  >  1:1430340/1‑149 (MQ=255)
                                           ggcAAAAAAATTCATCAGTGGTATTACCGCCCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCt                                                                                                           <  1:128322/148‑1 (MQ=255)
                                              aaaaaacTTCATCAGGGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCTTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTgcg                                                                                                        <  2:1090329/149‑1 (MQ=255)
                                               aaaaaaTTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTgcgc                                                                                                       >  2:837719/1‑149 (MQ=255)
                                                    aTTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGCCCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGc                                                                                                  >  2:421003/1‑149 (MQ=255)
                                                        atcaGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGc                                                                                              >  1:920301/1‑149 (MQ=255)
                                                         tcaGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGGGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCa                                                                                             <  2:319037/149‑1 (MQ=255)
                                                          caGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCTTGAGCTGGCAAACTGGCTGCGCGAAGCGGGCAt                                                                                            <  1:1256740/149‑1 (MQ=255)
                                                           aGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATc                                                                                           <  2:1124738/149‑1 (MQ=255)
                                                            gTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATca                                                                                          >  2:650291/1‑149 (MQ=255)
                                                              ggTATTACCGCGCCGATGATCTTGAGCATAAACCCGCGTTGCTGGTGCATCTGTTAACACAGCCGGAAGCGACCCGCGCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCGGGCAAACTGGCTGCGCGAAGCGGGCATcaac                                                                                        <  1:321396/149‑1 (MQ=255)
                                                               gTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATcaaca                                                                                       >  2:22333/1‑149 (MQ=255)
                                                                tATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATcaacaa                                                                                      >  2:911133/1‑149 (MQ=255)
                                                                 aTTACCGCGCCGATGATCTTGACCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCaacaac                                                                                     <  1:211405/149‑1 (MQ=255)
                                                                  ttACCGCGCCGATGATCTTGAGCATAACACCGCGGTGCTTGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCACTTGTGTTTGTGCGTAAGCTTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACt                                                                                    <  1:932953/149‑1 (MQ=255)
                                                                   tACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTg                                                                                   <  2:1492017/149‑1 (MQ=255)
                                                                        cgcCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTAtc                                                                              <  1:421003/149‑1 (MQ=255)
                                                                          ccGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTAtctc                                                                            <  2:885683/149‑1 (MQ=255)
                                                                            gatgatCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGa                                                                          >  2:1525236/1‑149 (MQ=255)
                                                                                    tGAGCATAACACCGCGTTGCTGGTGCATCTGATAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTgaga                                                                  <  2:870003/149‑1 (MQ=255)
                                                                                    tGAGCATAAAACCGCGTTGCTGTTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTgaga                                                                  <  2:247961/149‑1 (MQ=255)
                                                                                     gAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGGTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGAt                                                                 <  1:1021992/149‑1 (MQ=255)
                                                                                     gAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGGTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGAt                                                                 <  1:1016369/149‑1 (MQ=255)
                                                                                        cATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGt                                                              <  2:1355305/149‑1 (MQ=255)
                                                                                          tAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCACTTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTAc                                                            >  1:762266/1‑149 (MQ=255)
                                                                                           aaaaCCGCGTTGCTGTTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACa                                                           <  1:1452491/149‑1 (MQ=255)
                                                                                            aaaCCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAg                                                          >  2:1154180/1‑149 (MQ=255)
                                                                                              aCCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAggg                                                        <  1:853835/149‑1 (MQ=255)
                                                                                                cgcgGTGATGGTCCATCTGATAAAACAGCCGGAAGCCACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGTTACATGGCa                                                      <  2:301041/149‑1 (MQ=255)
                                                                                                       cTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTaa                                               <  2:159213/149‑1 (MQ=255)
                                                                                                        tGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTTTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAAc                                              <  2:441723/149‑1 (MQ=255)
                                                                                                        tGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGAGTAAc                                              >  1:1274181/1‑149 (MQ=255)
                                                                                                         ggTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACg                                             >  2:1149659/1‑149 (MQ=255)
                                                                                                          gTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGa                                            >  2:407644/1‑149 (MQ=255)
                                                                                                                 tGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATc                                     <  1:640508/149‑1 (MQ=255)
                                                                                                                     aaaaCAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGc                                 >  1:41741/1‑149 (MQ=255)
                                                                                                                       aaCAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGt                               >  1:945661/1‑149 (MQ=255)
                                                                                                                        aCAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGtt                              <  1:1525236/149‑1 (MQ=255)
                                                                                                                         cAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCTTAAGCGTGAGCGTGTGGATGAGCTGGCAAACTGGCGGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGGGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGttt                             <  2:798304/149‑1 (MQ=255)
                                                                                                                           gCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGa                           <  1:106208/149‑1 (MQ=255)
                                                                                                                            ccGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGAc                          >  2:1110414/1‑149 (MQ=255)
                                                                                                                            ccGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGAc                          <  2:725249/149‑1 (MQ=255)
                                                                                                                            ccGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGAc                          >  1:1375461/1‑149 (MQ=255)
                                                                                                                                aaGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGACCGaa                      <  1:650291/149‑1 (MQ=255)
                                                                                                                                aaGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGACCGaa                      <  2:41741/149‑1 (MQ=255)
                                                                                                                                         gCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGACCGAAGGTCGCGTa             >  1:312869/1‑149 (MQ=255)
                                                                                                                                            cAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGACCGAAGGTCGCGTAAAc          <  1:1397190/149‑1 (MQ=255)
                                                                                                                                             aaTTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAATAGGGAATCCACAACTGCTATCTCGAAGGTCAGATGGAACAGTGCAACCGTAACCAAGCCATCAAGCGATTGACCGAATGTCTCGTactcc         >  2:995056/1‑145 (MQ=255)
                                                                                                                                             aaTTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGACCGAAGGTCGCGTAAACg         >  2:696394/1‑149 (MQ=255)
                                                                                                                                             aaTTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGGTTGACCGAAGGTCGCGTAAACg         <  2:863813/149‑1 (MQ=255)
                                                                                                                                              aTTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGACCGAAGGTCGCGTAAACGt        >  2:159556/1‑149 (MQ=255)
                                                                                                                                                 gtgtTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGAACGAAGGTCGCGTAAACGTACt     >  1:1098734/1‑149 (MQ=255)
                                                                                                                                                  tgtTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGACCGAAGGTCGCGTAAACGTACTg    >  1:1197886/1‑149 (MQ=255)
                                                                                                                                                    tttGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATGTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGACCGAAGGTCGCGTAAACGTACTGGt  <  2:597589/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCGAAGCGGGCATCAACAACTGCTATCTCGAAGGTGAGATGGTACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGACCGAAGGTCGCGTAAACGTACTGGT  >  NC_000913/2713507‑2713803

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: