Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,345,152 A→G V136A (GTC→GCC)  yjdF ← conserved inner membrane protein YjdF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,345,1520AG98.6% 255.9 / ‑4.1 70V136A (GTC→GCC) yjdFconserved inner membrane protein YjdF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (39/30);  minor base T (1/0);  total (40/30)
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.87e-01

AACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGACCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTC  >  NC_000913/4345020‑4345288
                                                                                                                                    |                                                                                                                                        
aaCATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCACTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGaa                                                                                                                                   >  1:880605/1‑140 (MQ=255)
 aCATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAAc                                                                                                                                  >  2:1059517/1‑140 (MQ=255)
   atatCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGc                                                                                                                                <  2:463266/140‑1 (MQ=255)
      tCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGccac                                                                                                                             <  2:401265/140‑1 (MQ=255)
          aTTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCAtt                                                                                                                         >  1:94318/1‑140 (MQ=255)
          aTTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCAtt                                                                                                                         >  1:843524/1‑140 (MQ=255)
              cGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTAc                                                                                                                     >  2:265801/1‑140 (MQ=255)
              cGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTAc                                                                                                                     >  2:953289/1‑140 (MQ=255)
              cGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTAc                                                                                                                     >  2:432274/1‑140 (MQ=255)
              cGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTAc                                                                                                                     >  2:765505/1‑140 (MQ=255)
                   cccACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCcg                                                                                                                <  1:265801/140‑1 (MQ=255)
                   cccACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCcg                                                                                                                <  1:267167/140‑1 (MQ=255)
                   cccACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCcg                                                                                                                >  1:427480/1‑140 (MQ=255)
                    ccACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCcgc                                                                                                               <  1:432274/140‑1 (MQ=255)
                     cACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCcgcg                                                                                                              <  1:641439/140‑1 (MQ=255)
                     cACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCcgcg                                                                                                              <  2:937955/140‑1 (MQ=255)
                        tGGTCGCCCTGGGTGCCGAGAAAATCAACCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCAc                                                                                                           >  2:341813/1‑140 (MQ=255)
                                cTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCcg                                                                                                                <  2:877935/127‑1 (MQ=255)
                                cTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCcg                                                                                                                >  1:877935/1‑127 (MQ=255)
                                      gCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGaa                                                                                             <  2:94318/140‑1 (MQ=255)
                                         gagaAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACga                                                                                          >  1:426281/1‑140 (MQ=255)
                                         gagaAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACga                                                                                          >  1:502026/1‑140 (MQ=255)
                                             aaaTCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAAt                                                                                      >  2:1242028/1‑140 (MQ=255)
                                                        tCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGccac                                                                          <  2:992816/140‑1 (MQ=255)
                                                         cccTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGccac                                                                          <  1:622985/140‑1 (MQ=255)
                                                              aCCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTg                                                                     >  2:554811/1‑140 (MQ=255)
                                                                  aTTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAgg                                                                 >  2:201816/1‑140 (MQ=255)
                                                                  aTTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAgg                                                                 >  2:989558/1‑140 (MQ=255)
                                                                   ttGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCa                                                                   <  1:287193/137‑1 (MQ=255)
                                                                   ttGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCa                                                                   >  2:287193/1‑137 (MQ=255)
                                                                        aGCGCTGCCCACCACTCGATTAATTCATACATCGCACTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAg                                                           <  2:467242/140‑1 (MQ=255)
                                                                          cgcTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGcc                                                         >  2:297193/1‑140 (MQ=255)
                                                                          cgcTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGcc                                                         >  1:313113/1‑140 (MQ=255)
                                                                              gCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTg                                                     >  1:1272443/1‑140 (MQ=255)
                                                                               cccACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTgg                                                    >  2:307224/1‑140 (MQ=255)
                                                                                   ccacTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCTCCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGaaaa                                                <  2:847643/140‑1 (MQ=255)
                                                                                       tCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGTCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTg                                            >  2:901970/1‑140 (MQ=255)
                                                                                         gATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGcc                                          <  1:297193/140‑1 (MQ=255)
                                                                                         gATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGcc                                          >  1:103704/1‑140 (MQ=255)
                                                                                          aTTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGccc                                         <  2:103704/140‑1 (MQ=255)
                                                                                              aTTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGc                                     <  1:989558/140‑1 (MQ=255)
                                                                                              aTTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGc                                     >  1:743604/1‑140 (MQ=255)
                                                                                                   tACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTc                                >  1:787523/1‑140 (MQ=255)
                                                                                                   tACATCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTc                                >  2:527952/1‑140 (MQ=255)
                                                                                                      aTCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATa                             <  1:1059517/140‑1 (MQ=255)
                                                                                                      aTCGCGCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATa                             <  2:732767/140‑1 (MQ=255)
                                                                                                        cgcgCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGccc                                         <  2:859441/126‑1 (MQ=255)
                                                                                                        cgcgCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGccc                                         >  1:859441/1‑126 (MQ=255)
                                                                                                        cgcgCTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACg                           >  2:232516/1‑140 (MQ=255)
                                                                                                         gcgcTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTg                                  <  1:622361/132‑1 (MQ=255)
                                                                                                         gcgcTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTg                                  >  2:622361/1‑132 (MQ=255)
                                                                                                         gcgcTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACgg                          >  2:849167/1‑140 (MQ=255)
                                                                                                          cgcTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGa                         >  1:1243443/1‑140 (MQ=255)
                                                                                                          cgcTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGa                         >  1:783720/1‑140 (MQ=255)
                                                                                                          cgcTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGa                         >  2:9164/1‑140 (MQ=255)
                                                                                                          cgcTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGa                         >  1:487126/1‑140 (MQ=255)
                                                                                                          cgcTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGa                         >  2:1068143/1‑140 (MQ=255)
                                                                                                          cgcTTATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACAAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGa                         >  2:791268/1‑140 (MQ=255)
                                                                                                             ttATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTa                      <  1:232516/140‑1 (MQ=255)
                                                                                                              tATCGCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTAc                     <  1:307224/140‑1 (MQ=255)
                                                                                                                  gCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACt                 <  2:487126/140‑1 (MQ=255)
                                                                                                                  gCCAGCGCCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACt                 <  1:765505/140‑1 (MQ=255)
                                                                                                                       cgcCACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAAcc            >  1:406810/1‑140 (MQ=255)
                                                                                                                          cACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGAtt                       <  2:548104/126‑1 (MQ=255)
                                                                                                                          cACGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGAtt                       >  1:548104/1‑126 (MQ=255)
                                                                                                                            cGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGc       >  1:403426/1‑140 (MQ=255)
                                                                                                                            cGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGc       <  1:953289/140‑1 (MQ=255)
                                                                                                                            cGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGc       <  1:527952/140‑1 (MQ=255)
                                                                                                                            cGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGc       <  1:201816/140‑1 (MQ=255)
                                                                                                                            cGCAGCAGGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGc       >  2:1106775/1‑140 (MQ=255)
                                                                                                                                 cagGCCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTc  <  1:341813/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                        
AACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGACCAGGAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTC  >  NC_000913/4345020‑4345288

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