Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,477,882 T→G 7.1% L12L (CTT→CTG)  dsdX → D‑serine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,477,8820TG 7.1% 290.0 / 10.9 85L12L (CTT→CTG) dsdXD‑serine transporter
Reads supporting (aligned to +/- strand):  ref base T (38/41);  new base G (4/2);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

ATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCGG  >  NC_000913/2477736‑2478028
                                                                                                                                                  |                                                                                                                                                  
aTCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGCTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTAt                                                                                                                                                  >  4:631713/1‑149 (MQ=255)
    cTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGc                                                                                                                                              >  2:589820/1‑149 (MQ=255)
     tCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCa                                                                                                                                             >  1:115704/1‑149 (MQ=255)
       cATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATc                                                                                                                                           >  2:141328/1‑149 (MQ=255)
           gATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCgtgt                                                                                                                                       >  4:545661/1‑149 (MQ=255)
           gATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCgtgt                                                                                                                                       <  2:633280/149‑1 (MQ=255)
             tGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTa                                                                                                                                     >  2:90033/1‑149 (MQ=255)
             tGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTa                                                                                                                                     >  2:71575/1‑149 (MQ=255)
                aaTTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAAtt                                                                                                                                  <  4:99561/149‑1 (MQ=255)
                aaTTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAAtt                                                                                                                                  >  1:388521/1‑149 (MQ=255)
                 aTTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTg                                                                                                                                 >  2:562428/1‑149 (MQ=255)
                  tttACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGt                                                                                                                                <  3:601866/149‑1 (MQ=255)
                     aCCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACt                                                                                                                             <  1:562428/149‑1 (MQ=255)
                       cTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGa                                                                                                                           >  3:263997/1‑149 (MQ=255)
                        tCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGAc                                                                                                                          >  4:286965/1‑149 (MQ=255)
                        tCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGCTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGAc                                                                                                                          >  1:163369/1‑149 (MQ=255)
                             tgAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCg                                                                                                                     >  4:575147/1‑149 (MQ=255)
                                   gCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGt                                                                                                               >  3:389455/1‑149 (MQ=255)
                                   gCAAATTTTATCGTTTGTAAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCATCATTGTGTTAATTGTACTGACCATCGTGAAGt                                                                                                               <  2:509160/149‑1 (MQ=255)
                                        ttttATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCaaa                                                                                                          >  4:652894/1‑149 (MQ=255)
                                        ttttATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCaaa                                                                                                          <  1:71575/149‑1 (MQ=255)
                                             tCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCa                                                                                                     >  1:113925/1‑149 (MQ=255)
                                               gTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCAcc                                                                                                   >  3:451766/1‑149 (MQ=255)
                                                       gCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTg                                                                                           <  3:631713/149‑1 (MQ=255)
                                                                 tttttttGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGc                                                                                 >  1:245141/1‑149 (MQ=255)
                                                                   tttttGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTg                                                                               <  2:388521/149‑1 (MQ=255)
                                                                         tCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGc                                                                         >  4:654367/1‑149 (MQ=255)
                                                                         tCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACACGTTTCTGGCGCTGTTGCTGGCCAGc                                                                         >  2:655183/1‑149 (MQ=255)
                                                                           cAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCt                                                                                            <  2:424935/128‑1 (MQ=255)
                                                                           cAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGctt                                                                       >  4:406527/1‑149 (MQ=255)
                                                                           cAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGctt                                                                       >  4:528014/1‑149 (MQ=255)
                                                                           cAATAGCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCt                                                                                            >  1:424935/1‑128 (MQ=255)
                                                                            aaTATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGcttc                                                                      <  2:56307/149‑1 (MQ=255)
                                                                            aaTATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGcttc                                                                      >  4:262260/1‑149 (MQ=255)
                                                                                 catcaGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGt                                                                 <  3:406527/149‑1 (MQ=255)
                                                                                   tcaGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTgg                                                               <  2:115704/149‑1 (MQ=255)
                                                                                   tcaGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTgg                                                               <  4:355672/149‑1 (MQ=255)
                                                                                       gTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCa                                                                                                                                                <  1:350022/64‑1 (MQ=255)
                                                                                       gTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCa                                                                                                                                                >  2:350022/1‑64 (MQ=255)
                                                                                       gTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAAc                                                           <  4:503828/149‑1 (MQ=255)
                                                                                       gTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAAc                                                           >  3:421913/1‑149 (MQ=255)
                                                                                       gTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAAc                                                           >  3:429760/1‑149 (MQ=255)
                                                                                          aaTCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACgat                                                        <  1:655183/149‑1 (MQ=255)
                                                                                            tCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACgatga                                                      >  3:665888/1‑149 (MQ=255)
                                                                                            tCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACgatga                                                      >  4:317054/1‑149 (MQ=255)
                                                                                            tCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACgatga                                                      <  2:194870/149‑1 (MQ=255)
                                                                                             cacaGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACgatgat                                                     <  1:334424/149‑1 (MQ=255)
                                                                                             cacaGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACgatgat                                                     <  4:665888/149‑1 (MQ=255)
                                                                                             cacaGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACgatgat                                                     <  3:262260/149‑1 (MQ=255)
                                                                                             cacaGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACgatgat                                                     <  3:412443/149‑1 (MQ=255)
                                                                                              acaGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGatgatg                                                    <  3:652894/149‑1 (MQ=255)
                                                                                              acaGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGatgatg                                                    <  3:528014/149‑1 (MQ=255)
                                                                                              acaGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGatgatg                                                    <  4:422232/149‑1 (MQ=255)
                                                                                              acaGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGatgatg                                                    <  1:141328/149‑1 (MQ=255)
                                                                                              acaGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGatgatg                                                    <  3:286965/149‑1 (MQ=255)
                                                                                              acaGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGatgatg                                                    <  2:113925/149‑1 (MQ=255)
                                                                                                        gTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATggg                                            >  2:280946/1‑147 (MQ=255)
                                                                                                        gTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATggg                                            <  1:280946/147‑1 (MQ=255)
                                                                                                           agaTATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCAc                                       <  1:430231/149‑1 (MQ=255)
                                                                                                             atatGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGACAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTg                                     <  2:163369/149‑1 (MQ=255)
                                                                                                                tGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGat                                  <  4:419601/149‑1 (MQ=255)
                                                                                                                  cACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGatat                                >  3:685021/1‑149 (MQ=255)
                                                                                                                     tctcAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGt                             >  3:167193/1‑149 (MQ=255)
                                                                                                                     tctcAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGt                             >  1:76199/1‑149 (MQ=255)
                                                                                                                     tctcAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGt                             <  2:245141/149‑1 (MQ=255)
                                                                                                                     gctcAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAACTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGt                             >  2:587539/2‑149 (MQ=255)
                                                                                                                        cAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTaaa                          <  4:263997/149‑1 (MQ=255)
                                                                                                                        cAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTaaa                          <  1:90033/149‑1 (MQ=255)
                                                                                                                             cTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTa                     <  4:167860/149‑1 (MQ=255)
                                                                                                                             cTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTa                     <  4:685021/149‑1 (MQ=255)
                                                                                                                                   tgtgAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGaaa               <  3:545661/149‑1 (MQ=255)
                                                                                                                                    gtgAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAg              <  4:421913/149‑1 (MQ=255)
                                                                                                                                      gAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTg            >  3:522484/1‑149 (MQ=255)
                                                                                                                                       aacACGATGCTGATAAGCAGCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTgg           >  1:685062/3‑149 (MQ=255)
                                                                                                                                          aCGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAAt        >  1:631487/1‑149 (MQ=255)
                                                                                                                                            gatgctGATAAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTg      >  2:490559/3‑149 (MQ=255)
                                                                                                                                            gatgctGATAAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTg      >  2:180290/3‑149 (MQ=255)
                                                                                                                                            gatgctGATAAGCAGCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTg      >  2:95077/3‑149 (MQ=255)
                                                                                                                                              tgctTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGc    <  4:389455/149‑1 (MQ=255)
                                                                                                                                              tgctTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGc    <  4:522484/149‑1 (MQ=255)
                                                                                                                                              tgctTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGc    >  1:355555/1‑149 (MQ=255)
                                                                                                                                              tgctGATAAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGc    <  1:490559/149‑1 (MQ=255)
                                                                                                                                              tgctGATAAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGc    <  1:180290/149‑1 (MQ=255)
                                                                                                                                                ctTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCgg  >  4:33966/1‑149 (MQ=255)
                                                                                                                                                ctTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCgg  <  3:575147/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCGG  >  NC_000913/2477736‑2478028

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