Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,224,242 A→T 6.2% intergenic (+335/+37) ycgI → / ← minE protein YcgI/Z‑ring positioning protein MinE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,224,2420AT 6.2% 162.1 / 2.1 61intergenic (+335/+37)ycgI/minEprotein YcgI/Z‑ring positioning protein MinE
Reads supporting (aligned to +/- strand):  ref base A (27/28);  new base T (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

ATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTG  >  NC_000913/1224102‑1224383
                                                                                                                                            |                                                                                                                                             
aTACACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAAct                                                                                                                                             <  1:905481/143‑1 (MQ=255)
     cATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGcc                                                                                                                                        <  1:383231/143‑1 (MQ=255)
          tgttttgtGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGag                                                                                                                                   <  1:414352/143‑2 (MQ=255)
            ttttgtGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTATGCGTTGAAGa                                                                                                                                 >  2:263489/1‑143 (MQ=255)
            atttgtGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGAGGAAATAGCAAAGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGCAGa                                                                                                                                 >  1:496277/2‑143 (MQ=255)
                 tgggTGGCAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                                                                                            <  2:683675/140‑1 (MQ=255)
                     tGGGAACAGCCTTAAGGGGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGGATTCTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgcgc                                                                                                                        <  2:665376/143‑1 (MQ=255)
                        gAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCttt                                                                                                                     <  2:502966/143‑1 (MQ=255)
                          aCAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACCACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                                                                                                   <  2:834914/143‑1 (MQ=255)
                              ccTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGAGAAATGCGCTTTTACAGc                                                                                                               >  2:742315/1‑143 (MQ=255)
                               cTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCg                                                                                                              >  2:7744/1‑143 (MQ=255)
                               cTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCg                                                                                                              >  1:641874/1‑143 (MQ=255)
                                ttAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                                                                                             >  1:894260/1‑143 (MQ=255)
                                 tAAGGTGTACAGGGGGAGGTGGAAATAGCGATGAGGAGTATCAGCACGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCggg                                                                                                            <  2:225701/143‑1 (MQ=255)
                                 tAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCggg                                                                                                            <  1:73363/143‑1 (MQ=255)
                                 tAAGGTGTAAAGGGGGACGTGGCAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCggg                                                                                                            <  1:342992/143‑1 (MQ=255)
                                          aacagcgAGGTGGTCATAGCAATGAGGAGTATCAGCAAGAATACGCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATACCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                                                                   <  2:589191/137‑1 (MQ=255)
                                          aaGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                                                                   <  1:823906/143‑1 (MQ=255)
                                                 aGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttc                                                                                            <  2:923579/143‑1 (MQ=255)
                                                  ggTGGAAATAGCGATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttct                                                                                           <  1:475807/143‑1 (MQ=255)
                                                   gttgAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTTGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTg                                                                                          <  2:927166/140‑1 (MQ=255)
                                                    tGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGc                                                                                         >  1:242979/1‑143 (MQ=255)
                                                       aaaTAGCAATTAGGAGTATCAGCAAGATTACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTc                                                                                      <  1:382109/143‑1 (MQ=255)
                                                         atatCAAGGAGGAGTATGAGCAAGCATTCTCCCGACTTGCGCACAACTTGGATGAGAGGGATGAAAACCTCAAGGCAGAGATAACTCTGCCTTGAAGGTAAATGCGCTCTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCg                                                                                    <  2:394887/143‑1 (MQ=255)
                                                           aGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGCGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTAGAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGt                                                                                  >  1:37213/1‑143 (MQ=255)
                                                             ccaggaggagTATCGGCAAGAATACTCGCCGCTTCACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTaa                                                                                <  1:662168/139‑1 (MQ=255)
                                                             cAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCGTGAAGATGAATGCGCTTTTACAGCGGGCTTATTTCAGTTCTTCTGCttacagaaa                                                                                >  1:38455/1‑136 (MQ=255)
                                                                tGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGt                                                                             >  1:122881/1‑143 (MQ=255)
                                                                 gaggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTc                                                                            <  1:199269/143‑1 (MQ=255)
                                                                         ccAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTACCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAg                                                                    <  1:989159/142‑1 (MQ=255)
                                                                           aGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTATGCCTTGTAGATAAATGCGCTTTTACAGCGGGCTTCTTTCAGCTCTTCTGCTTCCGGTAAGGTCCCGTTCAGCt                                                                  >  1:695774/1‑143 (MQ=255)
                                                                            gCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCGTCCGGTAAGGTCACGTTCGGCTc                                                                 >  1:285059/1‑143 (MQ=255)
                                                                            gCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTACTTCCGGGAAGGTAGCGGTCAGCTc                                                                 >  1:104721/1‑143 (MQ=255)
                                                                                gAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGa                                                             >  2:138787/1‑143 (MQ=255)
                                                                                  aTACTCGCCGCTTTACCACAACGTGGATGGGAGGGATGAAAGACTCAAGGCAGAGATAACTCTGCCTTGAAGAGAAATGCGCTTTTACAGCGGGCTCATTTCAGCTCTTCTGCTGCCGTTCAGGTCACGTGCAGCTCGAGaac                                                           >  2:26076/1‑142 (MQ=255)
                                                                                   tACTCGCCGCTTTACCACAACGTGTATGAGAGGGATGAAAAACTCAAGGCAGAGATAGCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATa                                                          <  2:909702/143‑1 (MQ=255)
                                                                                   tACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATa                                                          <  2:595054/143‑1 (MQ=255)
                                                                                        gccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCATCTCAAGCATAGAAat                                                     >  1:438217/1‑143 (MQ=255)
                                                                                        gccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCACGGATAGaaag                                                     >  1:730963/1‑142 (MQ=255)
                                                                                           gcTTTACCGAAGCGTCGATGCGAGGGATGAAAAACTCAGGGCAGAGATAAGTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGTTTAGTTCAGCTCTTCTGCTTCCTGTAAGGTCACGTTTAGCTTAAGAATAGAAATATc                                                  <  2:932318/143‑1 (MQ=255)
                                                                                            cTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCg                                                 >  2:397747/1‑143 (MQ=255)
                                                                                              ttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGcc                                               <  1:946907/143‑1 (MQ=255)
                                                                                                 ccACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATc                                            >  1:742550/1‑143 (MQ=255)
                                                                                                    caACGTGGATGAGAGGGATGAAAACCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCttt                                         <  2:407986/143‑1 (MQ=255)
                                                                                                       cGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGATCTTCTGCTCACGGTAAGCTCACGTTCAGCACAAGAATAGAACAATCGCCctctacctt                                      >  2:837108/1‑138 (MQ=255)
                                                                                                          ggATGAGAGGGATGAACACCTCAAGGCAGAGATAACTCTGCTTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTc                                   <  2:888898/143‑1 (MQ=255)
                                                                                                             tGAGAGGGATGAAGCACGCAAGGCAGAGATAACTCTGCCTTGAAGATAAACGCGCTTTTACAGCATCTTTGCTTCTGTTTGTGTGGGCCCCACAAGGTCACGTTCAGCTCAAGATTAGGATTATCGCCATCTTTTTGCTCAAg                                <  2:1077129/143‑1 (MQ=255)
                                                                                                             tGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAg                                <  1:93358/143‑1 (MQ=255)
                                                                                                              gagaGGGATGAAAAACTCAAGGAAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGGCACGTTCCGCTCAAGAATAGAAATATTGCCATCTTTTTGCTCAAGc                               >  1:366599/1‑143 (MQ=255)
                                                                                                                 agGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGCAATATCGCCATCTTTTTGCTCAAGCTGt                            >  1:37418/1‑143 (MQ=255)
                                                                                                                 agGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGt                            >  2:594942/1‑143 (MQ=255)
                                                                                                                       gAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAACTATCGCCATCTTTTTGCTCAAGCTGTTCGGtt                      >  1:914092/1‑143 (MQ=255)
                                                                                                                        aaaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCACGCATAGAAATATCGCCATCTTTTTGCTCAAGCTGTAGGGCTa                     >  2:358181/1‑143 (MQ=255)
                                                                                                                         acaaCTCAAGGCAGTGAGAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTAc                    <  1:821884/141‑1 (MQ=255)
                                                                                                                         aaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCCGCTCTTCTGCTTCCGGTAAGGTCACTTTAAGCTCAAAAATAAAAATATCACCATATTTTTATTCAAGTTGTACGGTTAc                    >  1:92014/1‑143 (MQ=255)
                                                                                                                         aaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTAc                    <  2:574988/143‑1 (MQ=255)
                                                                                                                             ttCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGACATATCGCCATCTTTTTGCTCAAGCTGTAGGGTGACCAtc                >  2:818548/2‑143 (MQ=255)
                                                                                                                                aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTCTTTGCTCAAGCTGTACGGTTACCATCTCa             >  2:260150/1‑143 (MQ=255)
                                                                                                                                aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCACCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCa             <  2:200983/143‑1 (MQ=255)
                                                                                                                                 aGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAg            >  1:44647/1‑143 (MQ=255)
                                                                                                                                         aTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAAtt    <  1:7744/143‑1 (MQ=255)
                                                                                                                                           aaCTCTGCCTTGTAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTCTTTGCTCAAGCTGTACGGTTCCCATCTCAGGATCAATTTg  >  2:90050/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                             
ATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTG  >  NC_000913/1224102‑1224383

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: