Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,765,565 A→G 5.9% M71T (ATG→ACG)  yfjM ← CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,765,5650AG 5.9% 206.8 / 2.5 67M71T (ATG→ACG) yfjMCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (36/27);  new base G (2/2);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

ATGGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAA  >  NC_000913/2765417‑2765712
                                                                                                                                                    |                                                                                                                                                   
atggcactatggcactatggcactatggcactatggcactatggcacCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACa                                                                                                                                                     <  2:90587/149‑1 (MQ=255)
      ctatggcactatggcactatggcactatggcactatggcacCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCgt                                                                                                                                               <  4:655895/149‑1 (MQ=255)
       tatggcactatggcactatggcactatggcactatggcacCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCgtg                                                                                                                                              >  2:647044/1‑149 (MQ=255)
          ggcactatggcactatggcactatggcactatggcacCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTc                                                                                                                                           >  3:31254/1‑149 (MQ=255)
            tactatggcactatggcactatggcactatggcacGGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATGATGTCCTGGTTGGTAGGGTCATATCCCGATGCTCTTTGAGCAATTTCAACATCGCGTGTTCAt                                                                                                                                         <  2:221921/148‑1 (MQ=255)
            cactatggcactatggcactatggcactatggcacCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCAt                                                                                                                                         <  1:597861/149‑1 (MQ=255)
             actatggcactatggcactatggcactatggcacCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACGTCGCGTGTTCATg                                                                                                                                        >  1:310220/1‑149 (MQ=255)
             actatggcactatggcactatggcactatggcacCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACGTCGCGTGTTCATg                                                                                                                                        >  3:88748/1‑149 (MQ=255)
             actatggcactatggcactatggcactatggcacCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCgt                                                                                                                                               <  1:472236/142‑1 (MQ=255)
             actatggcactatggcactatggcactatggcacCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCgt                                                                                                                                               >  2:472236/1‑142 (MQ=255)
             actatggcactatggcactatggcactatggcacCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATg                                                                                                                                        >  2:521642/1‑149 (MQ=255)
             actatggcactatggcactatggcactatggcacCGCCAGCGATTGTAATCACAATTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATg                                                                                                                                        >  1:452077/1‑149 (MQ=255)
                     actatggcactatggcactatggcacCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCta                                                                                                                                >  1:676154/1‑149 (MQ=255)
                       tatggcactatggcactatggcacCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCtata                                                                                                                              <  2:30594/149‑1 (MQ=255)
                        atggcactatggcactatggcacCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTatat                                                                                                                             <  3:515342/149‑1 (MQ=255)
                         tggcactatggcactatggcacCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGAGTTCATGGCTTTCTATAtt                                                                                                                            >  4:155944/1‑149 (MQ=255)
                            cactatggcactatggcacCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATAttgtt                                                                                                                         <  3:518496/149‑1 (MQ=255)
                                  ggcactatggcacCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGACTTTCTATATTGTTGATGcc                                                                                                                   >  2:110326/1‑149 (MQ=255)
                                       tatggcacCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGcc                                                                                                              <  2:347362/149‑1 (MQ=255)
                                                        ttGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTcc                                                                                             <  1:378059/149‑1 (MQ=255)
                                                           tAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTta                                                                                          <  3:678443/149‑1 (MQ=255)
                                                             aTCACACTTGATATTATGAAACACAGTTGCACGCATTATTTCCTGGCTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCgt                                                                                                                                               <  3:233602/94‑1 (MQ=255)
                                                             aTCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACGTCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTtata                                                                                        <  2:310220/149‑1 (MQ=255)
                                                             aTCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACGTCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTtata                                                                                        <  4:88748/149‑1 (MQ=255)
                                                             aTCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTtata                                                                                        >  2:318788/1‑149 (MQ=255)
                                                             aTCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTGGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCgt                                                                                                                                               >  4:233602/1‑94 (MQ=255)
                                                              tCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAg                                                                                       >  1:647535/1‑149 (MQ=255)
                                                              tCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAg                                                                                       >  1:421112/1‑149 (MQ=255)
                                                              tCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAg                                                                                       >  1:354249/1‑149 (MQ=255)
                                                              tCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAg                                                                                       >  2:326517/1‑149 (MQ=255)
                                                                 cacTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCtc                                                                                    >  2:699926/1‑149 (MQ=255)
                                                                 cacTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCtc                                                                                    >  2:296826/1‑149 (MQ=255)
                                                                    ttGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCtcatc                                                                                 >  2:434364/1‑149 (MQ=255)
                                                                       atatTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCatg                                                                              >  1:57642/1‑149 (MQ=255)
                                                                                aaCACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGt                                                                     >  4:51940/1‑149 (MQ=255)
                                                                                  cacaGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGttt                                                                   <  2:299304/149‑1 (MQ=255)
                                                                                   acaGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGtttt                                                                  >  2:578946/1‑149 (MQ=255)
                                                                                     aGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTgg                                                                <  1:647044/149‑1 (MQ=255)
                                                                                       ttGCACGCATTATTTCCTGGTTGGTAGGGTCGTATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTTTTCATGGCTTTCTATATTGTTGATGCCTTGCCAATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGc                                                              >  2:626374/1‑149 (MQ=255)
                                                                                            cGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGCTTTGGGctctct                                                         <  2:452077/149‑1 (MQ=255)
                                                                                             gCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTg                                                        >  4:186105/1‑149 (MQ=255)
                                                                                                   tttCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGtt                                                  <  1:521642/149‑1 (MQ=255)
                                                                                                   tttCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGtt                                                  <  3:638897/149‑1 (MQ=255)
                                                                                                   tttCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGtt                                                  <  3:597142/149‑1 (MQ=255)
                                                                                                       cTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGt                                              >  3:178756/1‑149 (MQ=255)
                                                                                                        tggttggtAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATcc                                                                                                         <  1:583605/89‑1 (MQ=255)
                                                                                                        tggttggtAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATcc                                                                                                         >  2:583605/1‑89 (MQ=255)
                                                                                                        tggttggtAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGtt                                             >  2:416137/1‑149 (MQ=255)
                                                                                                                 gggTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAg                                    >  1:506525/1‑149 (MQ=255)
                                                                                                                  ggTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGa                                   >  1:132974/1‑149 (MQ=255)
                                                                                                                      atatCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAAttt                               >  4:370734/1‑149 (MQ=255)
                                                                                                                       tatCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAActtt                              >  2:609531/1‑149 (MQ=255)
                                                                                                                          ctcGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTc                           >  1:280488/1‑149 (MQ=255)
                                                                                                                           tcGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCt                          >  4:84717/1‑149 (MQ=255)
                                                                                                                            cGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCtt                         <  1:434364/149‑1 (MQ=255)
                                                                                                                             gATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTg                        <  3:330033/149‑1 (MQ=255)
                                                                                                                             gATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTg                        >  1:588476/1‑149 (MQ=255)
                                                                                                                              aTGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGt                       <  3:51940/149‑1 (MQ=255)
                                                                                                                              aTGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGt                       <  2:421112/149‑1 (MQ=255)
                                                                                                                                gCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGACTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTct                     >  2:483495/1‑149 (MQ=255)
                                                                                                                                    tttGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTctgctg                 >  1:652110/1‑149 (MQ=255)
                                                                                                                                         gCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGaa            >  1:106731/1‑149 (MQ=255)
                                                                                                                                             tGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGaataat        >  4:486009/1‑149 (MQ=255)
                                                                                                                                                  aCATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTa   <  3:155944/149‑1 (MQ=255)
                                                                                                                                                   cATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTaa  <  2:588476/149‑1 (MQ=255)
                                                                                                                                                   cATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTaa  <  1:318788/149‑1 (MQ=255)
                                                                                                                                                   cATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAACTTTTTCTTGTCTGCTGGCGAATAATTACTaa  <  1:609531/149‑1 (MQ=255)
                                                                                                                                                   cATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTCCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTaa  <  1:451806/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATGGCACTATGGCACTATGGCACTATGGCACTATGGCACTATGGCACCGCCAGCGATTGTAATCACACTTGATATTATAAAACACAGTTGCACGCATTATTTCCTGGTTGGTAGGGTCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAA  >  NC_000913/2765417‑2765712

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑