Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,391,463 G→A 7.4% R131C (CGC→TGC)  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,4630GA 7.4% 204.7 / 5.2 68R131C (CGC→TGC) ycjYS9 homolog non‑peptidase family protein
Reads supporting (aligned to +/- strand):  ref base G (28/34);  new base A (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 = 2.80e-01

TCCGCTGCTAATATCGCTGGTTCTGGCGTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAAC  >  NC_000913/1391317‑1391608
                                                                                                                                                  |                                                                                                                                                 
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                                                                                                                                           <  4:200928/128‑1 (MQ=255)
                           gTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGaa                                                                                                                                           >  3:200928/1‑128 (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                                                                                                                    >  5:44879/1‑121 (MQ=255)
                                                         aTCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCt                                                                                                                    <  6:44879/121‑1 (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   <  6:24448/149‑1 (MQ=255)
                                                                                                                                              cggcgATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTaa   >  7:3466/1‑149 (MQ=255)
                                                                                                                                               ggcgATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAAc  >  3:260330/1‑149 (MQ=255)
                                                                                                                                               ggcgATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAAc  <  2:234950/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TCCGCTGCTAATATCGCTGGTTCTGGCGTTTGAACCCGCTTCAACGTATGGCAGGGCATCAATAGATTTCACATTGTTTTCCCAACCATTACGGAATATCGAACCAATATTGACCGCACTAACGGTACCGATGGCCTTGATGCGGCGATCCTGAATTGCAGCATTGGCTGTATATCCTGCACCGGCACAAATTCCCATCGCACCAATTCGGGTATTGTCGACATAAGAAAGCGTTGTCAGGTAATCAATCACGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAAC  >  NC_000913/1391317‑1391608

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: