Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,391,460 G→T 7.4% R132S (CGC→AGC)  ycjY ← S9 homolog non‑peptidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,391,4600GT 7.4% 206.5 / 5.7 68R132S (CGC→AGC) ycjYS9 homolog non‑peptidase family protein
Reads supporting (aligned to +/- strand):  ref base G (29/34);  new base T (2/3);  total (31/37)
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.33e-01

TTCTCCGCTGCTAATATCGCTGGTTCTGGCGTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAAC  >  NC_000913/1391314‑1391608
                                                                                                                                                  |                                                                                                                                                    
ttCTCCGCTGCTAATATCGCTGGTTCTGGCGTTTGAACCCGCTTCACCGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATgcgg                                                                                                                                                     <  6:238049/148‑1 (MQ=255)
   tcCGCTGCTAATATCGCTGGTTCTGGCGTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGAt                                                                                                                                                 >  4:276125/1‑149 (MQ=255)
     cGCTGCTAATATCGCTTGTTCTGGCGTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGAt                                                                                                                                                 <  7:245091/147‑1 (MQ=255)
          ctAATATCGCTGGTTCTGGCGTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCTGCAAttcaggat                                                                                                                                          >  1:174771/1‑142 (MQ=255)
            aaTATCGCTGGTTCTGGCGTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTg                                                                                                                                        >  1:312222/1‑149 (MQ=255)
             atatCGCTGGTTCTGGCGTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAAtt                                                                                                                                         <  5:4127/147‑1 (MQ=255)
              tttCGCTGGTTCTGGCGTTTGATCCCGCTTCAACGTATGTGAGGGCTTCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTgca                                                                                                                                      <  6:216803/147‑1 (MQ=255)
              tatCGCTGGTTCTGGCGTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTgca                                                                                                                                      <  5:48436/149‑1 (MQ=255)
                   cTGGTTCTGGCGTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCAtt                                                                                                                                 <  7:113264/149‑1 (MQ=255)
                      gTTCTGGCGTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCTGCAATTCAGGATCGCCGCAtcaag                                                                                                                              <  5:252002/149‑5 (MQ=255)
                      gTTCTGGCGTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCTGCAATTCAGGATCGCCGCAtcaa                                                                                                                               <  2:174771/148‑4 (MQ=255)
                       ttCTGGCGTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCTATCCTGAATTGCAGCATTGGCt                                                                                                                             >  5:371105/1‑149 (MQ=255)
                       gtCTGGCGGTTGAACCCGCTTCCACGTATGGCAGGGCATCAATAGATTTCGCATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCt                                                                                                                             <  6:228002/148‑1 (MQ=255)
                            gCGTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCTCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGt                                                                                                                           <  2:268365/146‑1 (MQ=255)
                              gTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGaa                                                                                                                                           >  3:200928/1‑128 (MQ=255)
                              gTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGaa                                                                                                                                           <  4:200928/128‑1 (MQ=255)
                              gTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGCATTGCAGCATTGGCTGTATATc                                                                                                                      >  1:346437/1‑149 (MQ=255)
                              gTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATc                                                                                                                      >  4:125479/1‑149 (MQ=255)
                               tttGAACCCTCTTCAACGTATGGCAGGGCATCAATTGATTTCACATTGTTTTGCCATCCATTACGGAATATCGAACCAATATTGACCGCAGTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATcc                                                                                                                     <  8:324151/149‑1 (MQ=255)
                               tttGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCTGCAATTCAGGATCGCCGCAtcaag                                                                                                                              >  3:8409/1‑136 (MQ=255)
                               tttGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCTGCAATTCAGGATCGCCGCAtcaag                                                                                                                              <  4:8409/140‑5 (MQ=255)
                                     cccGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCt                                                                                                                             >  3:40329/1‑135 (MQ=255)
                                     cccGCTTAAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCt                                                                                                                             <  4:40329/135‑1 (MQ=255)
                                      ccGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTgcagc                                                                                                                                    <  4:52689/127‑1 (MQ=255)
                                      ccGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTgcagc                                                                                                                                    >  3:52689/1‑127 (MQ=255)
                                          ttCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTCCCGATGGCATTGATGAGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGca                                                                                                           >  6:319516/1‑148 (MQ=255)
                                               cgaaTGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAAACATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCACCATTGGCTGTATATCCTGCACCGGCACAAAtt                                                                                                     <  6:106186/146‑1 (MQ=255)
                                                gTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTc                                                                                                    <  3:125479/149‑1 (MQ=255)
                                                   tGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCa                                                                                                 <  3:315774/149‑1 (MQ=255)
                                                       aGGGCATCAATAGATTTCACATTGTTTTCCCAACAATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTgca                                                                                                                                      >  2:359692/1‑108 (MQ=255)
                                                            aTCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCt                                                                                                                    <  6:44879/121‑1 (MQ=255)
                                                            aTCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCt                                                                                                                    >  5:44879/1‑121 (MQ=255)
                                                               aatcgcattCACATTGTGTTCCCCACCCTTCCGGACTATCGACCCAATCTGGACCGCACCAACGGTACAGAGGGCCTTGATGCGGCGATCCGGAATGGCAGCAGTGGCGGTATATCATGCACCGGCACAAATTCCCCTCGCACCAAt                                                                                       <  2:79611/140‑1 (MQ=255)
                                                               aaTAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTc                                                                                     >  7:178062/1‑149 (MQ=255)
                                                                 tAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGTCCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTggg                                                                                   <  2:40546/149‑1 (MQ=255)
                                                                 tAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCTGCACAAATTCCCATCGCACCAATTAgg                                                                                   <  6:371105/149‑1 (MQ=255)
                                                                 tAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCgg                                                                                   >  5:174405/1‑149 (MQ=255)
                                                                   gATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCggg                                                                                  >  1:234950/1‑148 (MQ=255)
                                                                   gATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCGTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGt                                                                                 <  3:276125/149‑1 (MQ=255)
                                                                       tCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTAtt                                                                              <  2:312222/148‑1 (MQ=255)
                                                                           aTTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGCCCGCACTAACGCTACCGATGGCCATGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGa                                                                         >  4:270758/1‑149 (MQ=255)
                                                                           aTTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGa                                                                         <  5:107542/149‑1 (MQ=255)
                                                                           aTTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGa                                                                         >  3:192891/1‑149 (MQ=255)
                                                                            ttGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGAc                                                                        <  8:258528/149‑1 (MQ=255)
                                                                            ttGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGAc                                                                        <  2:309666/149‑1 (MQ=255)
                                                                               ttttCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACAt                                                                      >  1:11354/1‑148 (MQ=255)
                                                                               ttttCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACAt                                                                      >  5:24448/1‑148 (MQ=255)
                                                                                    ccAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGaaa                                                                <  1:141379/149‑1 (MQ=255)
                                                                                    ccAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGa                                                                  >  6:81241/1‑147 (MQ=255)
                                                                                     cAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTc                                                                                                    >  4:63659/1‑112 (MQ=255)
                                                                                     cAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTc                                                                                                    <  3:63659/112‑1 (MQ=255)
                                                                                      aaCCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGc                                                              >  2:137787/1‑149 (MQ=255)
                                                                                         cATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGtt                                                           >  2:25647/1‑149 (MQ=255)
                                                                                         cATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGtt                                                           <  8:178062/149‑1 (MQ=255)
                                                                                         cATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGtt                                                           <  4:192891/149‑1 (MQ=255)
                                                                                                     tCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTaatca                                               <  1:137787/149‑1 (MQ=255)
                                                                                                       gAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTaatcaa                                              <  1:25647/148‑1 (MQ=255)
                                                                                                         aCCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAatcaatca                                           <  5:81241/149‑1 (MQ=255)
                                                                                                           cAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACg                                         >  4:343840/1‑149 (MQ=255)
                                                                                                           cAATAATGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACg                                         >  2:211836/1‑149 (MQ=255)
                                                                                                                       ccctcaCGGTACCGATGGCCTTGCTGCGGCGATCCTGAATTGCAGCATGGGAGGTATATCCTGCAGCGGCACAACTTCCCAGCGCACCAATTCGGGTAGTGTCGACATAAGCAAGCGTTGTCAGGTAATCAATCACGGCAATGATGTCt                             >  7:269913/6‑149 (MQ=255)
                                                                                                                       cACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCt                             >  6:298568/1‑149 (MQ=255)
                                                                                                                        aCTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCAc                                                                                           >  4:86613/1‑86 (MQ=255)
                                                                                                                        aCTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCAc                                                                                           <  3:86613/86‑1 (MQ=255)
                                                                                                                        aCTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCtt                            <  2:346437/149‑1 (MQ=255)
                                                                                                                            aCGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCTTCa                          <  7:241493/147‑1 (MQ=255)
                                                                                                                                    gATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCTTCAGTACGAATg                 >  2:345978/1‑148 (MQ=255)
                                                                                                                                        gCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCTTCAGTACGAATGTATgg            <  1:211836/149‑1 (MQ=255)
                                                                                                                                               tGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTc      >  6:106487/1‑148 (MQ=255)
                                                                                                                                                gcggcgATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCt     >  6:160460/1‑148 (MQ=255)
                                                                                                                                                 cggcgATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTaa   <  2:11354/149‑1 (MQ=255)
                                                                                                                                                 cggcgATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTaa   >  7:3466/1‑149 (MQ=255)
                                                                                                                                                 cggcgATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTaa   <  6:24448/149‑1 (MQ=255)
                                                                                                                                                  ggcgATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAAc  >  3:260330/1‑149 (MQ=255)
                                                                                                                                                  ggcgATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAAc  <  2:234950/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTCTCCGCTGCTAATATCGCTGGTTCTGGCGTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAAC  >  NC_000913/1391314‑1391608

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: